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Chapter 1 



Motivation: Gravity in 2+1 dimensions 



Einsteinian gravity in 2+1 dimensions is much simpler than in 3+1 dimensions because 
it has no local gravitational degrees of freedom. This can be seen by considering the 
Riemartn tensor. In any dimension the Riemann tensor satisfies the following (anti)- 
symmetric relations under exchange of its indices 

■^ixvpa ^pcTfiv 

Consequently we may write the Riemann tensor in d dimensions as 

Rnupa = i^^u^faQAB, (1-2) 

where the u^^, with A = 1,.. . , ^d{d — 1) form a complete basis of 2-forms and Qab is a 
symmetric 2-tensor. That is, we can write the Riemann tensor as a symmetric tensor on 
the bundle of 2-forms. For a 3-dimensional manifold M, Poincare duality tells us that 
the bimdle of 2-forms, (TM), is isomorphic to the tangent bundle and that a complete 
basis is given by the Levi-Civita tensors, = e^^^. The Ricci tensor is therefore given 
by 



= re^p^^.J'Q 



(1.3) 



AB 

Consequently, we see that the Einstein tensor is equal to G^^ = —Qp.v The Riemann 
tensor in 3 dimensions simply is the Einstein tensor. Einstein's equation therefore com- 
pletely fixes the curvature in terms of the energy-momentum tensor. 
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In particular, if the energy-momentum tensor vanishes so does the Riemann tensor. 
Vacuum solutions of the Einstein equations in 2+1 dimensions are flat and do not have 
any local structure. This means that there are no local gravitational degrees of freedom 
or gravitational waves. 

If spacetime is simply connected, this implies that spacetime is diffeomorphic to (the 
covering space of an open subset of) Minkowski space. Non-local gravitational degrees 



of freedom do exist if the spacetime has a non-trivial topology. Witten 1 102 1 showed 
that, if spacetime has a topology of the form R x Sg, where Sg is a Riemann surface of 
genus g, the gravitational field equations can be solved exactly, and there is a (12(7—12)- 
dimensional space of solutions. The free parameters were identified as the holonomies 
of a flat ISO(2, 1) connection aroimd the non-contractable loops of the spacetime. 

By writing gravity in 2+1 dimensions as a Chem-Simon theory of ISO(2, 1)|^ Witten 
was able to provide a finite quantization of this system. Thereby he disproved earlier 
belief that gravity in 2+1 dimensions was non-renormalizible. 

The relatively simple setting of pure gravity in 2+1 dimensions (i.e. without any 
kind of matter content) has proven to be a valuable testing ground for approaches to 
quantizing gravity. Many different approaches have been applied. When compared, 
the results do not always appear to be equivalent. (See p3) for a review.) 

The idea behind Witten's approach was to first use the classical constraints to reduce 
the infinite number of gravitational degrees of freedom to a finite number, and then 
quantize the resulting system. This is the general idea behind 'reduced phase space' 
quantizations. Since in 2+1 dimensions the gravitational phase space can be reduced 
to a finite dimension, quantization becomes a matter of traditional quantum mechanics 
making such approaches particularly powerful. 

Different approaches to describing 2+1 dimensional gravity and reducing its phase 
space have been employed. Some, like Witten, have used a 'first order' description,|[97) 
while others have used 'second order' ADM formalism. fW, 57 1 Nelson and Regge have 
made an extensive study of the classical algebra of observables to use as the basis of a 
canonical quantization. | [63}]70) A yet other approach is to use geometrical techniques to 
make a complete classification of Lorentzian manifolds with constant curvature. [62, 95| 
These each lead to different approaches to quantizing the pure gravity theory, ||10[|21j|35j 



49 57 77 98 103 J which do not necessarily seem to agree.] 22] 



One thing to learn from this is that understanding the space of classical solutions 
is important for the quantization. In particular, quantizing the pure gravity theory 
divorced from its coupling to the matter content may not lead to the same results as 



1 A few years earlier Achucarro and Townsend |2 1 had already demonstrated this connection in the context 
of three dimensional supergravity. 
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quantizing the full theory. It is therefore essential to understand the classical space of 
solutions of gravity coupled to matter. 

1.1 Point particles 

Arguably the simplest form of matter one could consider is a collection of point par- 
ticles. The problem of describing moving point particles in (2+l)-dimensional gravity 
was first considered by Staruszkiewicz in 1963,|79| providing a solution for the two 
body problem. The general A^-body problem — assuming a topologically trivial back- 
ground — was solved by Deser, Jackiw, and 't Hooft in 1983.p4| 

In 2+1 dimensions the Schwarzschild metric describing the gravitational field of a 
point particle is 

dcr^ = -A{r)dt^ + A{r)-^dr^ + r'^d(j?. (1.4) 

The vacuum Einstein equation implies that the function A{r) is constant. If we set the 
constant value to (1 — a/(27r))^, then we can rescale t and r by a factor (1 — a/ (27r)) to 
obtain 

da^ = -df + dr^ + (1 - —fr^d^\ (1.5) 

This is the metric of a cone. The tip of the cone (at r = 0) is not a traditional curvature 
singularity since the Riemann curvature (and hence any curvature invariant) vanishes in 
any neighbourhood of the point. Its singular nature becomes apparent if one calculates 
the holonomy of a loop around the origin. The holonomy of a loop that encloses the 
origin once does not vanish as the loop is contracted to a point — as it should if the 
metric was regular at the origin. Instead it stays constant at a rotation of a degrees. The 
angle a is called the deficit angle of the conical singularity. 

As shown by Regge | [74) such a singularity should be associated with a finite amount 
of curvature, i.e. the curvature is a Dirac delta peak. As a consequence, the Energy- 
momentum tensor will also contain a delta peak at this point. In units where c = 1 
and G = g^, we find that the conical singularity is associated with a stationary point 
particle of mass a. This conclusion can be confirmed by smearing out the singularity to 
a smooth manifold and integrating the associated energy-momentum tensor. 

A cone can be constructed from the Euclidean plane by removing a wedge at the 
location of the particle and identifying the opposing edges of the wedge (see figure [rT| . 
The deficit angle is the angle of the removed wedge. This construction shows, in a 
geometrically explicit manner, that the space around the defect is locally flat. 

The conical defects constructed by removing a wedge are associated with positive 
curvature. Defects with negative curvature can be constructed from the Euclidean plane 



3 



Chapter 1: Motivation: Gravity in 2+1 dimensions 




by — instead of removing a wedge — inserting a wedge along a line (see figure 1.2 1. This 
creates a defect with a surplus angle, which may be interpreted as a point particle with 
negative mass generating negative curvature. 

By removing multiple wedges at different points one can create stationary configu- 
rations of multiple particles. If there are only positive mass particles, then the imiverse 
must be closed if the mass in any part exceeds 27r.||33| Application of the Gauss-Bonnet 
theorem then immediately implies that the total mass in the universe must be equal to 

47r. 

Any long range interaction between particles must be accompanied by a field, which 
has non-zero energy-momentum, and therefore generates its own curvature. Therefore, 
if we assume that point particles are the only source of curvature, the particles can only 
move along straight lines. 

The geometry generated by a moving particle can be obtained by Lorentz boosting 
the conical singularity of a stationary defect. If we start with a stationary particle at the 
origin, then the path of the particle through 3-dimensional spacetime after a boost B, is 
a straight timelike line that goes through the origin ait = 0. The wedge of spacetime 
that was removed at the particle also gets tilted (see figure 1.3 1. Before the boost the 
map identifying the two sides of the wedge was a rotation R, the holonomy of a loop 
around the origin. The map Q identifying the opposing sides of the wedge after the 
boost is a more general Lorentz transformation obtained by conjugation of R by B, i.e. 
Q — BRB^^. This is the holonomy of a loop around the moving particle. 
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Figure 1.2: Adding a wedge creates a conical defect with a surplus angle. 



The geometry corresponding to a particle located at another position than the origin 
at i = can be obtained from the previous case by applying an appropriate shift S. The 
map identifying opposing sides of the wedge no longer is a Lorentz transformation, but 
an element P of the Porncare group ISO(2, 1) obtained by conjugation of Q with 5", i.e. 
P = SQS^^. This is the so called Poincare holonomy of a loop around the defect]^ 

For a stationary particle, R represented the mass of the particle. After a boost B, 
the holonomy Q encodes both the energy and the velocity of the moving particle. In 
fact, as shown by Matschull and Welling,|[58) the holonomy of the defect may be inter- 
preted as the covariant momentum of the particle. The Poincare holonomy P captures 
the complete data of the trajectory of the particle. In particular, the trajectory of the 
particle is given by the set of fixed points of the Poincare holonomy. In this sense the 
Poincare group may be viewed as the covariant phase space of a point particle in (2+1)- 
dimensional gravity. 

1.1.1 Timelike loops 

In equation | |1.4^ we considered a Schwarzschild-like solution in 2+1 dimensions; the 
most general static and rotationally symmetric solution one can write down. More 

^Because the spacetime around the defects is completely flat, it is possible to parallel transport not just 
the tangent space but the entire coordinate frame along a curve. Transporting the coordinate frame along a 
closed loop will in general result in a coordinate frame that differs from the one you started with by a Poincare 
transformation. 
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Figure 1.3: A moving defect is obtained by boosting a stationary defect. 

generally, one can relax the condition of the metric being static to being stationary and 
allow a Kerr-like solution (see |32| for a derivation) 

dcr^ = -dt^ + 2ujdtd(j) + dr^ + [(1 -a/(27r))V2 _^2j^^2^ ^-^ 

This metric corresponds to a spinning point particle. The corresponding geometry can 
be obtained from Minkowski space by removing a wedge of spacetune like in the static 
case, but identifying the opposing sides with a shift in time. The Poincare holonomy 
of such a defect therefore includes a timelike shift — even in the frame in which it is at 
rest. 

This shift opens up the possibility of creating timelike loops. For example, in the 
above metric the curve given by the constant values r = < lj / {1 — and t = to is 
closed and timelike. Since closed timelike curves are irreconcilable with causality, this 
strongly implies that such sources should not be physically allowed. 

This, however, is not enough to bar the occurrence of closed timelike curves. As 
pointed out by Deser, Jackiw, and 't Hooft in 1983,||34j the geometry of a pair of particles 
moving past each other resembles that of a single pointlike source with the same angular 
momentum. In 1991, Gott |47J showed that if one constructs a pair of particles that 
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move past each other with a sufficient velocity, then it is possible to find closed timelike 
curves around the pair of particles. Deser, Jackiw, and 't Hooft |36J quickly pointed out 
that such a solution would have unacceptable boundary conditions at infinity. 

In principle, this left open the option that a Gott pair could form d5mamically from 
the decay of slow moving particles. Carroll, Farhi, Guth, and Olum f25',"26l showed that 
this would be impossible in an open universe, since no part of an open universe could 
contain enough energy to create a Gott pair. In the case of a closed universe 't Hooft 
| [85) showed by considering a Cauchy surface tesselated by polygons that although Gott 
pairs could form, the universe would collapse in a big crunch before a closed timelike 
curve could form. 

1.1.2 Quantization 

The exact solubility of a system of gravitating point particles in 2+1 dimensions invites 
a reduced phase space approach to quantizing the system. That is, one first classically 
solves the equations of motion for the gravitational field to obtain a restricted phase 
space for the point particles, and then quantizes this reduced system. This obviously 
requires a thorough understanding of the classical phase space. 

The polygon tessellation by Cauchy surfaces introduced by 't Hooft |85l allowed 
him to completely formulate the d5mamics of the system in terms of the edge lengths 
and their rapidities. The quantization of this system was studied in.|86| Because the 
Hamiltonian of the system is given by a periodic angle, it was observed that time was 
discretized after quantization. Due to Lorentz invariance one would expect space to 
be discretized as well, but this was not the case because the momenta (the rapidities 
of the edges) are not periodic. In a different parametrization using the coordinates of 
the vertices of the polygons rather than the edges length as coordinates, 't Hooft |90| 
found that the conjugate momenta lay on a sphere, and consequently foimd that in the 
canonical quantization spacetime was given by a lattice. 

Using different methods Matschull and Welling [58 1 showed that the covariant mo- 
mentum space of a gravitating point particle in 2+1 dimensions should be taken to be the 
spin-1 /2 representation of the Lorentz group rather than the x topology proposed 
by 't Hoof t|^ Besides the discretization of time, Matschull and Welling remarked that the 
curvature of momentum space leads to non-commuting operators for the coordinates in 
the quantum theory — non-commutative geometry seems to be generated naturally. 

Waelbroeck and Zapata | [99) have argued that the conclusion whether time is dis- 
cretized in the quantization of the polygon model depends on the details of the quan- 

^In fact, it is more inline witli 't Hoof t's covariant momenta used in 1 86 [ which used one angle and two 
hyperbolic angles. 
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tization procedure used. In particular, care needs to be taken in the implementation of 
gauge fixings. Further study of the phase space of the polygon model by Kadar |50^.^51J 
and Eldering | [40) have revealed further subtleties which complicate the quantization of 
the model. 

Meanwhile other approaches to quantum gravity in 2+1 dimensions have yielded 
conflicting results with respect to the presence of any fundamental discreteness. A 
loop quantum gravity based analysis |44| indicates that the length of spacelike intervals 
should be continuous, whereas the length of timelike intervals should be discrete. How- 
ever, a more recent quantization using Dirac variables associated to physical lengths and 
time intervals has found no indication of any discreteness . |20 J 

1.2 Towards 3+1 dimensions 

The broad idea of the work described in this thesis is to generalize the 2+1 dimensional 
model of gravitating point particles to 3+1 dimensions. The naive approach would be 
to study a system of gravitating point particles interacting according to the rules of 
general relativity in 3+1 dimensions. Of course, such an approach would not share the 
nice properties of the system of gravitating point particles in 2+1 dimensions. In fact, 
in 3+1 dimensions one cannot even restrict one's attention to point particles because 
their gravitational interaction would inevitably lead to the production of gravitational 
waves. 

Instead we take a more unconventional approach and attempt to construct a model 
that preserves some of the essential features of the model of gravitating point particles 
in 2+1 dimensions. In particular, we want to preserve the local finiteness of the model, 
which is key for its quantization. In essence this means that we want to preserve the 
property of gravity in 2+1 dimensions that the Einstein tensor completely determines 
the geometry. In particular, we want that regions of spacetime that are completely 
devoid of matter (i.e. where the energy-momentum tensor vanishes) are (Riemann) 
flat. 

We find that this is possible, if we study a model of propagating straight cosmic 
strings propagating in 3+1 dimensions according to the rules of general relativity. 

Although our motivation for studying this model is the possibility of formulating 
a theory of quantum gravity in 3+1 dimensions, we completely focus on the classical 
aspects of the model. As the history of quantization of 2+1 gravity shows, complete 
and correct knowledge of the classical phase space of a model is essential for its (re- 
duced phase space) quantization. Understanding the classical behaviour of the model 
is already quite a challenge. 
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The general plan of this thesis is as follows. Chapter |2] establishes the various foun- 
dational issues of the model studied, and describes a number of ways in which config- 
urations of straight cosmic strings may be parametrized. In chapter |3] we then study 
what happens in the model when two straight cosmic strings collide. It is found that 
the model may not be consistent for certain extreme collisions. Ignoring these issues we 
continue to study the continuum limit of the model using a linearized approximation 
in chapter |4j finding restrictions on the types of matter that may be approximated by 
the model. In chapter|5]we then show that the model reproduces gravitational waves as 
an emergent phenomenon in continuum limit. Finally, in chapter |6] we will discuss the 
open problems with the model and how these relate to a possible quantization. 

1.3 Conventions 

Throughout this thesis, we shall use the following conventions unless explicitly noted 
otherwise. 

• Spacetime metrics in 3+1 dimensions have a signature ( — h ++). 

• All quantities are expressed in natural units such that c — h= SttG = 1. 

• Lightcone coordinates u and v are defined as 



z + 1 



and V 



z 



- t 



(1.7) 



u = 



V2 ' 



V2 ■ 



In particular the Minkowski metric in lightcone coordinates is 



(1.8) 



• The Fourier transform of a function f{x) is defined as 




(1.9) 



Consequently, the inverse Fourier transform is 




(1.10) 
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Chapter 2 



Foundational matters 



Unlike in 2+1 dimensions, the Riemann tensor in 3+1 dimensions cannot be completely 
expressed in terms of the Einstein tensor. Consequently, the Einstein equation does not 
completely fix the geometry in terms of the matter content. As a result general relativity 
in 3+1 dimensions has geometrical local degrees of freedom, which are dynamical. 
These degrees of freedom make the (3+l)-dimensional theory much richer than the 
(2+l)-dimensional one. For example, they allow the appearance of gravitational waves 
and long range gravitational fields. 

However, these local degrees of freedom are also notoriously problematic when 
trying to quantize general relativity. There are many ways to phrase the problems that 
occur. A simple one is that from a field theory point of view these degrees of freedom 
form a spin-2 gauge field, the graviton. The coupling constant of this field, Newton's 
constant, has negative mass dimension. Consequently, a naive power coimting tells us 
that the resulting Feynman diagrams in perturbative quantum field theory will have 
non-renormalizable UV divergences. This can be confirmed by explicitly doing the loop 



integrations up to two loops. p6)|94) 

Many possible ways around this problem have been investigated over the years. 
For example, in supergravity [43] the hope is that the addition of new fields together 
with supersymmetry would help cancel the divergent terms |^ In other approaches, like 
string theory, general relativity and its metric only appear as an effective field theory 
of some more fundamental quantum theory that is renormalizable. Another possibility 
— expressed in theories like loop quantum gravity |[9)[76) — is that the renormalization 



*This hope was somewhat dampened by the discovery that the cancellations found do not generically 
persist beyond loop order. p7^.^52| More recently however, it was found that in the specific case of maximal 
(N = 8) supergravity the cancellations persist to much higher orders, leading to the suggestion that it might 
in fact be finite.|14||^[80] 
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problems disappear if the theory is quantized using a different set of variables. Yet 
others blame the use of perturbative methods and hope that the issues with perturbative 
quantum gravity may be circumvented by using a non-perturbative approach — for 
example as used in causal dynamical triangulations.Q 

Our take in this matter is that the local gravitational /geometrical degrees of freedom 
are themselves the root of the problem. The relative success of the quantization of (2+1)- 
drmensional gravity relies on the fact that a system having a finite number of particles 
will only have a finite number of degrees of freedom. Empty space is essentially feature- 
less. We would like to recreate this situation in a theory of gravity in 3+1 dimensions. 

Our approach to reach this will be somewhat heavy handed. We simply impose as a 
new principle that empty space should be featureless. 

Principle 1. No local structure. The Riemann tensor in empty spac^vanishes. 

At first sight, this is a substantial departure from conventional general relativity. 
How nevertheless general relativity is recovered in the continuum limit is the subject of 
chapters|4]and|5]of this thesis. In addition, we require that the Einstein equation contin- 
ues to hold in all of spacetime — including regions that are not empty. Combined with 
principle 1 this requirement puts very stringent limitations on the "matter" contents of 
the theory. For example, one cannot add point particles, because the Einstein equation 
would then imply that the surrounding empty space should have the Schwarzschild 
metric, which does not have a vanishing Riemann curvature. 

One can however add conical curvature defects. Conical defects have codimension 
2, so in 3+1 dimensions they are 2-dimensional planes. These are the only degrees of 
freedom that we will allow in the model. We want to interpret these 2-dimensional 
planes as the world sheet of a 1-dimensional straight line propagating through space. 
A timelike conical defect thereby corresponds to a line defect travelling at a subluminal 
speed, while a spacelike conical defect corresponds to a line defect travelling at a super- 
luminal speed. In the latter case, this means that there exists a Lorentz frame in which 
the spacelike plane lies in a single time slice. It would therefore represent a defect in- 
stantaneously appearing and disappearing in an extended region. This is unacceptable 
behaviour for a physical excitation from the point of view of local causality. This leads 
to the second pillar of our model: 

Principle 2. Local Causality. The only allowed degrees of freedom are non-spacelike conical 
curvature defects. 

The only allowed defects are therefore either (l+l)-dunensional or lightlike surfaces. 
We can view these as 1-dimensional defect lines propagating through space. Physically, 
they correspond to straight cosmic strings moving at a constant velocity. 

^We define "empty space" as a region of spacetime where the energy-momentum tensor vanishes. 
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Although our model is motivated as a step towards a renormalizable theory of quan- 
tum gravity, we are not actually going to discuss any quantization of the model. In- 
stead we are taking the lesson from (2+l)-dimensional gravity to heart, that thorough 
knowledge of the classical phase space is required for a proper quantization. The rest of 
this thesis is devoted to developing the classical description of a system of moving line 
defects, and discussing the subtleties that this involves. 

In this chapter we discuss the foundational aspects of the model, and the necessary 



formalism needed to describe the possible configurations. In section 2.1 we first discuss 



the geometry of a single stationary defect. We continue with the discussion of moving 



defects in section |2.2| Section 2.3 then discusses different ways of parametrizing a 
single defect. In particular we discuss the use of the holonomy of a loop around a 
defect to identify the defect. We discuss some different types of defects in section 

we 



2.4 



and elaborate on the possibility of massless lightlike defects. In section 2.5 



discuss geometries with multiple defects. This opens the possibility of junctions of 
defects, which are discussed in the following section | |2.6| . Instead of focussing on the 
curvature defects, the considered geometries can also be described with a focus on the 
local flatness by describing them as piecewise flat manifolds, this is discussed in section 
The final section of this chapter discusses the relation between the gravity model 



2.7 



discussed in this thesis and other piecewise flat approaches to gravity. 



2.1 Stationary defects 

We first describe a single stationary conical defect in 3+1 dimensions. Since the defect 
is stationary, all time slices are identical, and we just need to consider a conical defect in 
3-dunensional Euclidean space. We can obtain the geometry of such a defect as follows 
(see figure [ZT| . Like a conical defect in 2 dimensions, we start from a 3-dimensional 
Euclidean space, R'^. From this space we remove a single wedge, and identify the 
opposing sides — in this case planes. 

In flat polar coordinates (r, 9, z) the effect of removing a wedge is that at the location 
where the wedge was removed the azimuthal coordinate 9 will make a jump, say from 
(3 to (3 + a, where a is the angular size of the removed wedge, a; the deficit angle of the 
conical defect. The metric on the remaining region (including the time direction) is still 
the common Minkowski metric, 

ds^ ^~dt^ + dr^ + r^d9^ +Az^. (2.1) 

However, because the azimuthal angle 9 makes a jump, its period is reduced to 2tt — a. 
Alternatively, we can rescale the azimuthal coordinate Ohy a factor 1 — ^ such that its 
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Figure 2.1: A conical singularity in 3 spatial dimensions can he constructed by removing a 
3-dimensional wedge. 



period remains 2tt. As a result the metric becomes, 

ds^ = -df + dr^ + (1 - ^fr^de^ + dz^. (2.2) 

27r 

Note that this metric does not record the value of p. This angle, which gives the direction 
of the missing wedge, plays no role in the resulting geometry. In fact, one could have 
removed multiple wedges along the same axis, and obtained the same result as long as 
the total angle of the removed wedges was the same. 

Since this metric is locally isometric to the Minkowski metric, its curvature vanishes 
everywhere, except possibly at the coordinate singularity at r = 0. That the curvature 
does not vanish at r = can be inferred from calculating the holonomy of a loop around 
the r = axis. Consider a path 7'' (A) = (0, , A, 0). The holonomy Q^'^ e SO(3, 1) of 



14 



2.1: Stationary defects 



this path is 



Q^=cxp - / r;^,7'^(A)dA 



/I 0\ 

COS a sin a 
— sin a cos a 



Vo 



(2.3) 



1/ 



where F^^ is the Levi-Civita connection of the metric | |2.2) . Not coincidently, this rota- 
tion about the z-axis is the map that identifies the opposing sides of the wedge|^ Since 
this loop has a non-trivial holonomy, there must be curvature somewhere inside the 
loop. Since the space outside the r = axis is completely flat, this curvature must be 
located on the axis. 

To calculate the curvature associated to the axis it is convenient to change to a set of 
coordinates where the r coordinate is rescaled such that the azimuthal part of the metric 
takes the canonical form r^dd^, 



2tt 



'■2tt — a 



Me^ + dz^ 



(2.4) 



In these coordinates, we can "smear out" the singularity at r = by turning a into a 
function of r, 



for r < e 
otherwise. 



(2.5) 



For r > e this is the same metric as before, but as r approaches zero the metric con- 
tinuously approaches the Minkowski metric. The holonomy of a loop of radius r is a 
rotation about the z-axis of adr). The holonomy thus becomes trivial at the r = axis, 
and the space is regular 

For this non-srngular metric we can calculate the Einstein tensor. For r > e, the 
metric is unchanged, and the Einstein tensor vanishes. For r < e the Einstein tensor is. 



G flu — 







V 



2 



(2.6) 



7re ^ r 27re 



^More correctly, it is the pull-back of that map to the tangent space, but since locally we are in Minkowski 
space these two can be identified. 
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The total curvature of the smeared out singularity is therefore. 



G 



2'Kr 



27r-ao^ 






Vo 









-aj 



(2.7) 



This is independent of the smearing parameter e. Consequently, we can associate a 
delta peaked curvature to the conical singularity. If we set SttG = 1, then the Einstein 



equation identifies 1 2.7 1 as the energy-momentum tensor of the conical defect. This is 
the well-known energy-momentum tensor of a stationary straight cosmic string (see for 
example |T7J): it has a linear mass/ energy density equal to its deficit angle a, and an 
equal tension. 

Note that in this discussion there is no need for a to be a positive number. Geometri- 
cally, a negative value of a corresponds to opening the space along some hypersurface 
and inserting a wedge of angle |a|, which creates a defect with a surplus angle. Nor- 
mally, cosmic strings with a negative mass would be highly unstable, since negative 
mass implies a positive pressure along the string. Any deviation from straight would 
cause the string to buckle immediately. In our model, however, the principle that the 
space surrounding the defect is locally flat stabilizes negative mass defects because it 
forces the defects to be always straight; the destabilizing fluctuations simply cannot 
occur within the model. 

Still, one might object to the model including objects with negative mass, since these 
are in obvious violation of the various energy conditions that one generally expects 
realistic classical matter to satisfy. However, in the model that we are describing, the line 
defects are not just describing matter excitations, they also represent the gravitational 
excitations of the model. If this model is to have any hope of representing vacuum 
spacetimes at large scales, where on average the Einstein tensor vanishes, but the Weyl 
tensor has a non-zero average, then it needs to have both defects that add positively 
and negatively to the Einstein tensor. In chapter |5] we will see that having both types of 
defects is, in fact, essential in reproducing gravitational waves. 



2.2 Moving defects 

Since the spacetime around a stationary defect is locally Minkowski, we can obtain the 
geometry of a moving defect by applying a Lorentz boost to the geometry of a stationary 
defect. A (3+l)-dimensional wedge can be parametrized by four (unit) 4-vectors: two 
4-vectors uo and dp that span the leading edge of the wedge and two 4-vectors ei and 
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Figure 2.2: A (3+l)-dimensional wedge can be parametrized by four (unit) 4-vectors: two 
4-vectors vo and do that span the leading edge of the wedge and two 4-vectors ei and ei 
perpendicular to vo and do that give the directions of the sides of the wedge. 



62 perpendicular to vq and do that give the directions of the sides of the wedge. For the 
wedge removed from a stationary defect, vq = (1, 0, 0, 0) and do, ei, and 62 are chosen 
to lie in the t = hyperplane. 

The sides of the wedge are 3-dimensionaI half-hyperplanes spanned by vq, do, 
and ei, 

51 = {di (r, A, ji) = TVo + Xdo + Atei |r, A e R, e R>o} 

(2.8) 

52 = {cT2{r, A, /z) = TVo + Ado + A*e2|r, A e R, /x e R>o} • 

The conical curvature defect is created by identifying the sides Si and S2 through the 
identification 

o-i(t, A,/u) ~ cr2(r, A,/z) (2.9) 

for all T, X, and ^. We obtain the construction for a moving defect by applying a Lorentz 
boost B to this whole set-up. That is, we create a new wedge parametrized by the 
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vectors 

v'q = B - Vo 
d'o = B-do 
e[= B- ei 

e2=B- 62, 

and identify the edges, 

S'l = WiiT,>^,tJ-) = tv'q + \d'f) + iie[\T,\e R,^ e R>o} 
^2 = {'72(t, A,^) = tvq + Ado +/ie2|r, A e R,^^ e IR>o}, 



(2.10) 



(2.11) 



through the identification 

ai(T,A,M)^'T^(r,A,M) (2.12) 

for all r. A, and /i. 

The first thing we note is that, if the boost B is in the do direction, then — B-Ci — 
and this procedure is simply a reparametrization of the hjrperplanes Si and S2. This 



has no effect on the identification 1 2.9 1. We thus see that conical defects are invariant 
imder boosts in the direction of the defect line. 

This allows us to restrict our attention to boosts in the eie2-plane. Let B,, be the boost 
with velocity 77 e span(ei, 62). Then B,, leaves do invariant, while 

v'o = B^-vo^j{vQ + r]), ^^^^^ 
e- = B,, ■ a = eij_^ + 76^11,, + 7(6^ • r])vo, 

where 7 is the Lorentz factor of the boost, (1 — |?7|)^^/^ and e^j^^ and eiy,, denote the 
parts of Ci that are respectively perpendicular and parallel to 77. In general, the identifi- 
cation | |2.12) will identify points that do not lie in the same t' = constant plane. More 
specifically, the difference in t' is given by 

Ai'(r, A,/i) = ((T2(t, A,^) - a[{T,X,fi)) ■ vq 

= ^(^2 - e'l) • Vq (2.14) 

= Ai7(ei • 77 - 62 • 1]). 

Therefore, the identification | |2.12[ | will only identify points on the same time slice if 
ei • ?7 = 62 • ?/ or equivalently if 7/ c>c ei + 62- That is, if 77 is directed in the symmetry plane 
of the wedge. 

Recall, that earlier we observed that the direction of the wedge was irrelevant to the 
produced geometry. We are therefore free to always choose the wedge to be cut along 
the movement direction of the defect. This choice has the advantage that we know that 
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( 2.12^ only identifies points on the same time slice. This precludes the creation of closed 
timelike curves by the identification, and allows us to analyse the generated geometry 
on a single time slice. 

In fact, the identification | 2.12^ also produces zero shift in the cIq ~ d'^ direction, so 
we can analyse the new geometry completely in the eie2-plane. This plane is defined 
by the equations 

vo ■ X — 0, do ■ X = 0. (2.15) 

Let 

^ = {a',=fie',\fi€n>o} (2.16) 
be the intersection of and the eie2-plane. We can then find e'^ by solving 

vq ■ (t'^{t,X,^) = 0, do ■ a-ir, X, fi) ^ (2.17) 

for r and A. The result after normalizing is 

^> _ <- (e» ■ ri)v'o 
- (e, • r^Y 

The effective angle of the boosted wedge in the eie2-plane, a' , is then defined through 

a' e- • ri 
cos — = -^T—Y 

7|ry|^l- (ei-77)2 

If we observe that Ci ■ -q — \r]\ cos ^, we can apply some elementary trigonometry to 
produce a relation between the angle of the stationary wedge a, and the effective angle 
a' of the boosted wedge in the eie2-plane, 

a' a 

tan — = 7 tan — . (2.20) 

The Lorentz contraction of the boost causes the effective angle of the wedge to open up. 
As a result moving defects have larger effective deficit angles than stationary ones with 
the same mass. 



2.3 Parametrizing a defect 

To describe a general defect we need the codimension 2 subspace of the defect, and the 
deficit angle in its rest frame a. A codimension 2 subspace in d-dimensional Minkowski 



19 



Chapter 2: Foundational matters 




Figure 2.3: On the left: the eie2-plane before a boost rj. On the right:the eie2-plane after a 
boost rj. 

space is described by its normal 2-form cj^y and its displacement X^. Together these 
define the subspace through the equation 

a;^,(x^ - X^') = 0. (2.21) 

The normal 2-form is the wedge product of two perpendicular normal vectors. In d + 1 
spacetime dimensions this gives us 2(i — 2 free parameters. Equation | |2.21) is invariant 
under a shift X^' ^ X^' + X'^, with uJ^^,X'^' = 0. The displacement X^' therefore 
contributes another (d + 1) — (rf — 1) = 2 independent parameters. We therefore need a 
total of 

(2(i - 2) + 2 + 1 = 2(i + 1 (2.22) 

independent parameters to describe a codimension 2 conical defect. In our case, d = i, 
this is 7 parameters. 

The condition that the defect is timelike is captured by requiring that both normals 
to the codimension 2 subspace are spacelike. In terms of the normal 2-form this means 
that 

a;^.w^" = 1. (2.23) 

In our (3 + l)-dimensional case it is often useful to parametrize the spacetime ori- 
entation of the defect by a timelike and a spacelike 4-vector and d^^, instead of the 
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normal 2-form cj^^. If we choose these such that u'^ = 1, (f* = u^d^ = 0, then we can 
write 

= d^ = (0,d), (2.24) 

where v is the velocity of the defect with respect to the equal time slices and d is the 
direction of the defect line. Moreover, we can use the translation freedom in the dis- 
placement such that — 0, and write = (0, X) with X the position of the defect at 
t = 0. 

We can therefore describe a defect with mass a by a triple of 3-dimensional vectors 
{X,d,v) which satisfy 

V A- d, 

X _L d, and (2.25) 

lldll^l. 

2.3.1 Holonomy 

A particularly useful way of characterizing a defect is by the holonomy of a loop around 
it. The holonomy of a loop is found by parallel transporting a vector around the loop 
and comparing it to the original. The relation between the original and parallel trans- 
ported vector will be given by a Lorentz transformation, Q^^, called the holonomy. For a 
space M with an arbitrary metric the holonomy of a loop 7 : S*^ — > A/ can be calculated 
using the following formula (see [24J for a derivation), 

Q". = V exp ^dA^ , (2.26) 

where V indicates that the matrix exponential is path ordered. If the space is flat, then 
the holonomy does not depend on the details of the path taken. In particular, for any 
loop that can be contracted to a point without meeting a curvature defect the holonomy 
is the identity map. 

For our piecewise flat model this means that the holonomy of a loop around a defect 
only detects the number of times the loop wraps around the defect. This means that we 
can talk about the holonomy of a defect, meaning the holonomy of a loop wrapping the 
defect exactly once|^ 

In equation | |2.3^ we already calculated that the holonomy of a static defect was given 
by a rotation about the axis of the defect line. The angle of this rotation was found to be 
equal to the linear mass density of the defect. This holonomy coincides with the map 

^There is still an ambiguity here regarding the sign of the wrapping number, corresponding to the direction 
in which the loop wraps the defect. 
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that identifies the two sides of the removed wedge. This is true in general. For each 
wedge of spacetime there is a unique Lorentz transformation, Q, that maps Si to S2 in 
such a way that 

Q • tTi(r, A,/i) = (T2(t, A,^). (2.27) 

As this map is used to identify two points on the sides of the wedge, its pullback to 
the tangent bundle is used to identify the tangent spaces at those points. Therefore, a 
vector that is parallel transported across the identified sides of the wedge is transformed 
by this pulled back Lorentz transformation, which may be identified with the original 
identification since the tangent space of Minkowski space is isomorphic to Minkowski 
space. Since the rest of the loop is in flat space no further transformation with respect 
to the background frame occurs in parallel transport. Hence, the holonomy of a loop 
around a conical defect is equal to the Lorentz transformation associated with the re- 
moved wedge. 

So what is the holonomy of a moving defect? A defect moving at a subluminal 
speed can always be transformed into a stationary defect by a boost B. As we already 
observed, a stationary defect has a pure rotation as its holonomy. So, if a moving defect 
is constructed by identifying the sides Ei and S2 of a wedge, then applying B gives the 
sides — B ■ of a stationary wedge. The sides are identified by a rotation R such 
that 

R- a[{T,\,fj,) ^ a2{T,X,fi). (2.28) 

This implies that 

R-B-ai{T,X,n) = B-a2{T,X,fi). (2.29) 

Consequently, the Lorentz transformation Q that identifies the sides Ei and E2 is given 
by 

Q^B-^ R B. (2.30) 

The holonomy of an arbitrary moving defect therefore is a Lorentz transformation that 
is conjugate to a pure rotation. Such Lorentz transformations are called rotah'on/fA:e|^The 
angle of the pure rotation is the linear mass density of the defect in its rest frame, and 
can therefore be associated with the rest mass of the moving defect. 

Conversely, given a rotationlike Lorentz transformation Q we can construct a mov- 
ing defect with holonomy Q. Like a pure rotation, a rotationlike Lorentz transformation 
has two eigenvectors with eigenvalue 1, which span a timelike surface. This surface will 
be the leading edge of the wedge needed in the construction of the defect. We can now 
pick an arbitrary spacelike vector ei perpendicular to this surface and construct Ei as 
the span of the leading edge and ei. We define 62 as 62 = Q • ei. Since Q is an orthogonal 

'Similarly, Lorentz transformations that are conjugate to a pure boost are called boostlike. 
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transformation and Q leaves the leading edge invariant, 62 will also be perpendicular to 
the leading edge, and we can construct S2 as the span of the leading edge and 62. This 
defines a wedge whose identifying map, by construction, is Q. Therefore, if we remove 
this wedge from Minkowski space we obtain a defect with holonomy Q. 

The holonomy Q therefore completely encodes all the information about the mass 
and the movement of a line defect. Following the example of point defects in 2 + 1 
dimensions,! 58 1 we could interpret Q as the covariant 4-momentum of the line defect. 

Since we are in a locally flat background, we can go a step further than just calculat- 
ing the parallel transport of vectors. The flatness of the spacetime allows us to use the 
exponential map to expand each frame of the tangent bundle to a complete coordinate 
frame that is isometric to the Minkowski frame. Parallel transporting a frame around a 
loop, can therefore be interpreted as transporting the entire coordinate system around 
the loop. Consequently, the coordinate system obtained after parallel transport around 
a loop, will be related to the original coordinate system by an isometry of the Minkowski 
frame, a Poincare transformation. This transformation is called the Poincare holonomy of 
the loop. Like the ordinary holonomy it does not depend on the local details of the loop 
and can be interpreted as a property of the defect. 

If we use a coordinate system in which the studied defect passes through the origin 
of the coordinate system, then the Poincare holonomy will just give the ordinary holon- 
omy. A defect that passes through any other point may be constructed by applying a 
shift to the whole system. Going through the same motions as we did for determining 
the holonomy of a moving defect, we conclude that a moving defect with an arbitrary 
position produces a rotationlike Poincare holonomy, i.e. a Poincare holonomy that is 
conjugate to a pure rotation by a Poincare transformation. Conversely, given a rota- 
tionlike Poincare transformation P, we can construct a defect with Poincare holonomy 
P. 

The Poincare holonomy of a loop around a defect depends on the coordinate frame 
chosen at the initial point of the loop. When choosing a different coordinate frame, the 
holonomy of the loop is transformed by conjugating with the Poincare transformation 
associated to the change of frame. The only truly frame independent property of a 
defect, therefore is the conjugacy class of its holonomy. For rotationlike Poincare trans- 
formations, the conjugacy classes can be distinguished by the angle of the rotation in 
the rest frame of the transformation (i.e. the frame where the transformation becomes 
a pure rotation). We already saw that this angle was proportional to the linear mass 
density of the corresponding defect in its rest frame. The invariant of a defect given by 
the conjugacy class of its holonomy can therefore be associated with its rest mass. 

There is a subtlety in this association. The rotation angle of a rotationlike Poincare 
transformation can only distinguish deficit angles (and therefore mass densities) mod- 
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ulo 27r. This is not a problem if we only allow defects with a positive deficit angle, since 
a deficit angle of 27r corresponds to a spacetune with no volume. We, however, also 
wish to allow defects with a negative deficit (or surplus) angle, which in principle are 
not boimded by any value. 

The Poincare holonomy of a defect therefore contains (nearly) all the information 
needed to identify a defect. At the start of this section we saw that for complete descrip- 
tion of a conical defect we needed 7 independent parameters, one for the deficit angle 
and six to describe an arbitrary codimension 2 subspace. We should be able to reproduce 
these from the holonomy. The space of rotationlike Poincare transformations is indeed 
a 7 dimensional subspace of the 10 dimensional Poincare group. So, a rotationlike 
holonomy has the right number of parameters. 

We just saw that the deficit angle can be obtained as the conjugacy class of the 
holonomy. The codimension 2 subspace, S, forming the defect is obtained as the set 
of fixed points of the Poincare holonomy P, 

S={x€]R.^''^\P{x)=x}. (2.31) 

More specifically, the Poincare group can be interpreted as the semi-direct product of the 
Lorentz group with the abelian group of spacetime translations. As such, each Poincare 
transformation may be identified with a Lorentz transformation Q and a spacetime 
translation T. The Lorentz transformation Q associated to the Poincare holonomy P 
captures the spacetime orientation of the defect, while the displacement is captured by 
the associated spacetime translation T. 

2.4 More general defects 

In the preceding section we have seen that there is a close relation between defects and 
the Poincare holonomy of a loop around it. Thus far we have considered only defects 
with a rotationlike holonomy. One may wonder if it is possible to construct defects with 
more general holonomies. The answer is yes, but not all of these defects are physical. 

In the previous sections we constructed a conical defect by removing a wedge from 
a flat spacetime and identifying the opposing sides. This can be generalized in the 
following way. Take two hyperplanes Si and S2 intersecting on a codimension 2 surface 
E, remove one of the wedges of spacetime bounded by Si and S2, and identify Si and 
S2 using a Poincare transformation P that maps Si to S2. This requires that the leading 
edge of the wedge E is invariant under P. Moreover, if the resulting geometry is to 
be a regular topological manifold along E, then P restricted to E is to be the identity 
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map. That is, the leading edge of the wedge E must consist of fixed points of P. This 
requirement severely restricts the possible transformations P. 

Without loss of generality we may assume that E passes through the origin, other- 
wise we may apply a shift to the whole system such that it does. The requirement that 
the points in E are fixed under P implies that P has no translational part. Moreover, the 
Lorentz part of the holonomy Q, must have E as an eigenspace with eigenvalue 1. 

If E is timelike, as in the case of a timelike conical defect, then Q must be a rotation- 
like transformation. There are no other timelike defects that are the fixed point of their 
own holonomy, than the conical defects that we have already discussed. In the case that 
E is spacelike Q must be a boostlike transformation. This corresponds to a defect that is 
completely spacelike, which we do not allow since it would correspond to a non-local 
degree of freedom. 

The final possibility is that E is lightlike, in which case Q is a null rotation. The 
resulting defect can be interpreted as a line defect propagating at the speed of light. 
There is no good physical reason to disallow such objects, and in fact they will prove 
crucial in the construction of gravitational wave solutions in section |5] These defects 
will be discussed in the next section |2XT|. 



Not imposing the requirement that P is the identity when restricted to E implies that 
points in E will get identified with other points of E. The result generically is a space 
that is not a topological manifold along the points of E. For us this is enough reason not 
to consider such defects in our model. But let us briefly comment on what such defects 
would represent physically. 

Lets assume that E is timelike. Without loss of generality we may then take E to be 
the tz-plane. Even if the points of E are not fixed under P, E must still be an invariant 
subspace of P. This implies that P restricted to i? is a Poincare transformation of the 
tz-plane; i.e. it is a combination of a boost in the z-direction and shifts in the t and z 
directions. 

If P is a pure shift in the te-plane, it is possible to write down a locally flat metric for 
the resulting spacetime, 

ds^ = -{dt + (3d9f + dr^ + r^de^ + {dz + -fd<j)f, (2.32) 

with 13 the shift in the i-direction and 7 the shift in the z-direction. Smearing out 



the defect like we did for ordinary conical defects in section 2.1 reveals that /3 can 
be associated with a rotation of the defect, while 7 is associated to a constant torsion 
exerted on the defect. 

If /3 > 7, then the loops with constant t, r, and z become timelike for small r. That 
is, rotating infinitely thin defects generate closed timelike curves, just like rotating point 
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particles in 2+1 dimensions. This can actually be seen directly from the action of the shift 
on E. If the shift is tunelike, E becomes compactified in a tunelike direction. Similarly, 
if we consider a P that reduces to a pure boost on E, then this boost will identify points 
on E with a timelike separation. Consequently, there will be closed timelike curves on 
E. By continuity these can be deformed to closed timelike curves in a neighbourhood 
of in the quotient manifold. 

We conclude that allowing defects with a holonomy that includes a boost or a time- 
like shift would introduce unacceptable acausal features in the model. Even though 
they would allow the introduction of physically interesting concepts like intrinsic spin 
to the model. Defects with a spacelike shift in the holonomy seem less harmful, yet they 
do introduce a weird topology in the neighbourhood of the defect, which seems difficult 
to reconcile with the interpretation of the defect as a line defect. 



2.4.1 Massless defects 



There is no physical reason not to allow defects, where the 2-dimensional surface of the 
defect is lightlike. Since the defect surface consists of the fixed points of its holonomy, 
this means that the (Lorentz part of the) holonomy should be a Lorentz transformation 
that has a lightlike and an orthogonal spacelike eigenvector with eigenvalue one. If 

in Cartesian coordinates these are chosen to be (1, 1, 0, 0) and (0, 0, 0, 1), then the most 
general orientation preserving orthogonal transformation that leaves these vectors in- 
variant is 
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(2.33) 



where a e R is a free parameter. Lorentz transformations of this type are called null 
rotations or parabolic Lorentz transformations, and a is called the parabolic angle of the 
transformation. 

All null rotations belong to the same conjugacy class of the Lorentz group. They do 
not form a closed subgroup of the Lorentz group. In fact, the smallest subgroup of the 
Lorentz group that contains all null rotations is the Lorentz group itself. The largest 
subgroups that can be formed from just null rotations are two dimensional abelian 
groups isomorphic to (D generated by two commuting nilpotent elements of so(3, 1). 

In general a null rotation can be specified by giving a null vector, u'^, a spacelike vec- 
tor, perpendicular to and one scalar parameter a called the parabolic angle. Since 
the scaling of u/' and does not matter this gives 4 independent real parameters. The 
subset of null rotations does, however, not form a manifold due to a nodal singularity 
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at the identity element. Much like that the lightcone is not a proper submanifold of 
Minkoswki space. 

The subset of null rotations forms the boundary of the subset of rotationlike trans- 
formations in the Lorentz group. This means that any null rotation can be viewed as 
the limit of a sequence of rotationlike Lorentz transformations. For example the null 
rotation 1 2.33 1 can be viewed as the limit of a sequence of rotations about the z-axis 
boosted in the x direction. If Rz{4>) is a rotation of ((> about the z-axis and Bx{rj) is a 
boost in the x-direction with rapidity 77, then 
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(2.34) 



is a rotation about the z-axis boosted in the x-direction. For any constant value of (p 
(not equal to an integer multiple of 2tt) the components of this transformation diverge 
as rj approaches infinity. To approach a finite limit, cf) must approach zero as ij goes to 
infinity If we write 4> = ^^°[^ ^ , then using that 
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we find that 
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(2.37) 



This limiting procedure is analogous to applying an Aichelburg-Sexl ultraboost |^ to 
a line defect. The rotationlike transformation | |2.34| corresponds to a line defect oriented 
in the z-direction with rest mass density </> and moving in the x-direction with rapidity 
77. As 77 increases the kinetic energy of the line defect diverges. The choice (f) = ^^"^^ 
precisely ensures that the energy density of the line defect stays constant as 77 increases. 
We can therefore interpret lightlike defects as massless conical defects moving at the 
speed of light. 
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The Aichelburg-Sexl ultraboost of a line defect was studied by Barrabes et al in.|12| 
They found the metric corresponding to the holonomy | |2.33) to be. 



ds^ = -d<2 + dx^ + dy^ + dz^ - 2a\y\8{x - t)(dx - dif 



(2.38) 



The corresponding energy-momentum tensor can be found by taking the energy 
momentum tensor for a stationary defect as calculated in equation ||2.7^, 
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(2.39) 



and boosting that in the a;-direction. We then find that the energy-momentum tensor 
corresponding to the holonomy | |2.34 1 is 
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keeps the Too component constant. The limit as 77 goes to infinity 
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(2.41) 



Alternatively, we could have started directly from the metric 1 2.38 1, and calculated the 
Einstein tensor directly. This was done in |[12| and gives the same result. 

Note that in the massless limit the tension in the spacelike direction of the defect 
vanishes. This means that the energy-momentum tensor at a single point of the defect 
does not have any information on the direction of the defect. This property of massless 
defects will be very important when we construct gravitational waves in chapter |5] 
Because the energy-momentum tensor of massless strings is the same for all defects 
moving in the same direction, we can get the energy-momentum for perpendicular 
positive and negative energy defects to cancel each other. 



2,5 Multiple defects 

Geometries with multiple defects can be generated by removing multiple wedges from 
Minkowski spacetime. This introduces a new subtlety when describing the defects in 
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Figure 2.4: The path 71 and path 72 around the defect line A are topologically inequivalent. 
The path V shows that 71 can be deformed to the sequence of paths jb then 72 then 7^"*^. 



terms of their holonomies. In a spacetime with multiple defects, there are multiple 
topologically distinct loops around a defect line. Consider the situation depicted in 
figure [Z4] (for visual clarity the removed wedges have been suppressed). There are two 
topologically distinct loops from the point p around the defect line labelled A. The loop 
can either pass above the line labelled B, as the path labelled 71 does, or it can pass 
below the defect B like loop 72 . The loop 72 cannot be deformed into loop 71 without 
crossing defect B. Since the holonomy of a loop changes when a path is deformed across 
a defect, the holonomy of path 72 is different from the holonomy of 71 . 

We can however relate the holonomies of the paths 71 and 72, by observing that the 
loop r formed by first taking the loop 75 around defect B, then taking loop 72 around A, 
and finally taking 75 around B in the opposite direction, can be continuously deformed 
to path 71 without crossing any of the defects. Consequently, the holonomy of T is equal 
to the holonomy of 71 . Since the holonomy of F is the product of the holonomies of the 
individual loops, we conclude that if we denote the holomies of the individual loops 
Q71/ Q'y2' ^rid Q-ys, then they are related as through, 

Qji — Qr 

= Q-,^^Qj2Q7B (2.42) 

~ '37b '372 QlB ' 
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where in the last line we used that the holonomy of a loop followed backwards is the 
inverse of the holonomy of the loop. In particular, we see that Q-y^ and Q^^ belong to 
the same conjugacy class. Therefore if either one is rotationlike, so is the other. 

This relation can be expressed in formal way for a general configuration. Let X be 
a spacetime containing an arbitrary configuration of defects, let X'-^^ c X by the 2- 
dimensional subset of all defects, and let p e X \ X^^'> be a point not on any of the 
defects. Two loops in X starting at p can be continuously deformed into each other 
without crossing any of the defects if and only if they are homotopic in X\X^^K The set 
of classes of topologically equivalent loops together with the operation of concatenating 
loops therefore gives the fundamental group tti [X \X (2), p]. Since there is no curvature 
on X \ X'-^^ all homotopic loops have the same holonomy. We therefore have a map 



X\X'^^\p\ ^lSO{TpX), (2.43) 



that assigns to each homotopy class of loops 7 its holonomy Q~^, an element of the 
Poincare group at the point p; lSO{TpX). The consistency requirement of the example 
above generalizes to the requirement that the assignment Q is a group homomorphism. 
That is, for any two (equivalence classes of) loops 71 and 72, 

Qii-i2 — QfiQj2- (2.44) 

The holonomy of each simple loop (i.e. a loop enclosing a single defect once) must 
be rotationlike in order for all the defects to be timelike. Equation | |2.44| tells us that 
the holonomy of a loop enclosing multiple defects is the product of the holonomies 
of simple loops around the individual defects. Consequently, because the subspace of 
rotationlike Poincare transformations is not closed under multiplication, the holonomy 
of a general loop will not be rotationlike. In fact, rotationlike Poincare transforma- 
tions generate the whole Poincare group (in the sense that the smallest subgroup of the 
Poincare group containing all rotationlike transformations is the Poincare group itself). 
This means that a loop taken around a suitable set of rotationlike defects can have any 
holonomy. 

A homomorphism form ■ki[X \ X'^^\p] to the Poincare group can be defined by 
specifying an element of the Poincare group for each generator of tti [X \X^^\p'\. If X is 
contractible , then vri , p] is generated by a set of simple loops around the defect]^ 

Since the fundamental group of X \ X'^^' is torsion-free, the number of generators will 
be equal to the first Betty number of X \ X^^\ In the special case that X^^^ consists of a 
disconnected collection of infinite defecta^the number of generators will exactly equal 



'If the fundamental group of X is not trivial, the assignment Q also includes global "topological" degrees 
of freedom. As was the case in 2+1 dimensional pure gravity. 

^'^Configurations with finite connected defects will be considered in the next section. 
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Figure 2.5: If two defects meet at an angle, then there must exist a third defect (in white) at 
that junction. 

the number of defects. Since each of the simple loops must be assigned a rotationlike 
holonomy this gives 7n degrees of freedom for n defects. 

The assignment Q of holonomies to eqmvalence classes of loops depends on the 
choice of base point p and the choice of Poincare frame at that point. Another choice 
of point and frame will lead to an assignment Q that differs by conjugation with an 
element of the Poincare group. The choice of the base point and its frame can therefore 
be used to eliminate up to 10 degrees of freedom. For example a configuration with two 
defects has 2-7-10 = 4 degrees of freedom. Note that the choice of base point can 
never eliminate the frame independent degrees of freedom given by the conjugacy class 
of each generator. 

2.6 Finite defects and junctions 

Thus far we have only considered defects of infinite extent. More generally we can 
also consider finite defects, i.e. defects that have a finite length either in their spacelike 
direction of in their timelike direction (or both). To get a feeling for these objects let us 
first consider finite conical defects in 3-dimensional Euclidean space, where they are 1- 
dimensional lines. In a locally flat background it is impossible for a defect line to simply 
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end in mid space. So, finite defect lines have to end on other objects. In our piecewise 
flat model the only available objects to end on are other finite defects. 
If two finite defect lines end at one point there are two options: 

1. Both defects have the same deficit angles and same direction. I.e. the defects form 
a single longer line defect. 

2. There is at least one more defect ending at the same point. 



The second situation is shown in figure 2.5 The two defect lines (shown in black) meet 



at an angle. As a result the line where the deficit angles meet (shown in red) cannot be 
flat, but instead forms a third line defect. 

The same situation can occur in 3+1 dimensions. Conical defects cannot end simply 
in mid air, but instead must end on a codimension 3 submanifold called a junction, 
where at least two other bounded conical defects must end. The junction is a line that is 
shared between the defects ending there. Therefore, if Qi, . . . , Qn are the holonomies of 
defects ending at a junction, then the points of the junction must satisfy, 

Qix = Q2X = • • • = QnX = X. (2.45) 

We can distinguish three types of junction based on their metric signature: timelike 
subluminal jimctions, lightlike null junctions, and spacelike superluminal junctions. 

Of these the superluminal junctions appear troublesome. If three or more timelike 
defects end on a spacelike line, then there exists a Lorentz frame in which this line lies 
in a constant time slice. In this frame, the defects ending at this junction may be divided 
into two classes depending on whether they extend to the future or the past of this time 
slice (see figure |Z6j. If there is a mismatch between the number of defects in the future 
and the past of the junction, then the junction implies a change in the number of local 
degrees of freedom. 

In principle, this does not need to be the problem. If the number of defects decreases 
at the junction, this would be an explicit example of information loss. Although a bit 
peculiar, this might simply be a feature of the model. More troublesome, is the case 
where the number of defects increases at the junction. In that case there is a spontaneous 
non-local creation of information at the junction, which seems to violate local causality. 
Yet, if there are at least two defects in the past of the junction, then the junction describes 
the frontal collision of two (or more) collinear defects. It may yet be that the proper way 
to continue such a singular collision involves the creation of new defects. (In chapter |3] 
we will see that in general collisions are accompanied by the creation of new defects.) 

The situation becomes more dire if there is just one defect (or none at all) in the 
past of the junction. In that case the geometry to the past of the junction has no prior 
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Figure 2.6: Four different types of superluminal junction. The black line is a spacelike 
junction, the grey sheets are defects ending at the junction. From left to right:... 

a) ... two incoming and two outgoing defects, constant number of degrees of freedom. 

b) ... less outgoing than incoming defects. Information is lost. 

c) ... more outgoing than incoming defects. Information is created, but the incoming defects 
do predict the location of the junction. 

d) ... just one incoming defect. Information is spontaneously created. Not compatible with 
local causality. 



information about the junction at all. The junction simply appears instantaneously at 
all points of the junction. It is impossible to reconcile such an event with local causality. 
Since local causality was one of the primary assumptions of our model, appearance 
of such junctions would indicate a possibly fatal inconsistency of the model. Much of 
chapter |3] will deal with trying to avoid the appearance of such superluminal junctions. 

Another issue we saw in the example of line defects in 3-dimensional Euclidean 
space was that not just any set of defects can connect at a junction. In the example, 
the third line defect was determined in terms of the other two line defects. In 3+1 
dimensions, there exists a similar consistency condition for general junctions. This 
condition is a result of the requirement that the space around the junction is locally 
flat. It is most easily expressed in terms of the holonomy. Recall that the assignment 
of holonomies to each loop is required to be a group homomorphism with respect to 
the concatenation of loops. The topology of a junction is such that one can find simple 
loops around each of the "legs" of the junction, which can be concatenated to a single 



contractible loop (see figure 2.7 for an example with a 3-valent junction). This implies 
that the product of the holonomies of these loops must be the identity element. I.e. if 
the loops are labelled, 71 , . . . , 7„, then the corresponding holonomies satisfy 

Q71 • • • Qi„ = I- (2.46) 
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Figure 2.7: The loop F is liomotopic to both the concatenation of paths 71 • 72 • 73 and to the 
trivial loop. The holonomies assigned to the loops 7^ must thus satisfy Q-y-^ Qj^ Q73 = I- 



This equation is referred to as the junction condition of the holonomies. 

A second restriction comes from the requirement that all defects share a common 1- 
dimensional invariant subspace. This means that all holonomies must satisfy equation 
(Z45). This condition is not completely independent of the junction condition ||Z46). If 
the junction condition 1 2A6\ is satisfied for n defects, then if n — 1 holonomies share a 
common invariant line then so does the last one. 



2.7 Piecewise flat manifolds 

The description of a configuration of conical defects in terms of the Poincare holonomies 
of loops is very powerful. In particular, it gives a succinct formulation of the consistency 
conditions at a junction. However, this method also has a few drawbacks. We already 
remarked that the holonomy does not quite uniquely describe a defect configuration if 
we include defects with negative curvature. 

A second problem is that the path dependence of the holonomy means that the local 
geometry of a defect is encoded in a non-local way depending on the path chosen to 
reach the defect. Keeping track of the chosen paths becomes increasingly inconvenient 
for larger configurations of defects. 
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An alternative approach focuses on the geometry of the locally flat patch of space- 
time rather that the curvature of the defects. In this approach the geometry of a config- 
uration of defects is described as a piecewise flat manifold. 

Definition 1. A piecewise flat manifold is a manifold X with a CW complex structure and 
a flat (pseudo)-Riemannian metric g" on each n-cell e^. A geometry on X is defined by the 
following requirements 

1. The metric g" is smooth on the interior ofeJ^ and piecewise smooth on the boundary de^. 

2. The attachment map (j) a : SJ^^^ X'^^^^^^ of each n-cell e^^^to the {n — \)-skeleton X^^^^^^ 
is an isometry. 

3. For each {n — l)-cell e^~^ in the image of (pa, the inverse image (/>Q(e^~^) has zero 
extrinsic curvature in e^- 

These conditions ensure that each n-cell is an ?i-dimensional polytope. An n-cell can 
be subdivided by adding an (n— l)-cell dividing it in two pieces. Repeated subdivisions 
can be used to ensure that the cells satisfy some additional properties. For example, one 
could subdivide the cells until all cells are convex polytopes, or ultimately until all cells 
are simplices. 

Subdivision into a simplicial complex has the additional advantage that a flat ge- 
ometry on the interior of a n-simplex is uniquely determined by the lengths of the 1- 
simplices in its boundary. This leads to a more usual definition of a piecewise flat man- 
ifold in terms of a simplicial complex. |29 J Here we opt for the more general definition 
above since it gives more flexibility. 

If we take the manifold to be four dimensional and require that the metrics on the 
4-cells have Lorentzian signature, then X can be viewed as a piecewise flat spacetime 
consisting of polytope building blocks glued together along their edges. 

Because the attachment maps are isometries, for each (n — l)-cell e^^^ in the bound- 
ary of a n-cell there is an inverse isometry : e^~^ — >^ e^. Moreover, because X is a 
manifold, each 3-cell is included in exactly two 4-cells. 

Each 4-cell is isometric to a piece of Minkowski space We can fix a frame in a 4- cell 
by choosing a specific embedding ; — > R'^+^. Given two 4-cells and that 
share a 3-cell e^^ in their common boundary and framings La and Lp, there is a unique 

i^In the pathological cases where this is not possible, it is possible to subdivide the 4-cells into cells that are 
isometric to a subset of Minkowski space. 
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Poincare transformation P^a that makes the following diagram commute: 

R3+1 ^ ^ ^3+1 




The map Pap relates the change of frame as you move from cell a to cell /3. 

A general 2-cell is in the boundary of a finite number of 4-cells, say ai, . . . , a^. We 
can follow the change of frame as we foUow a loop 7 around the 2-ceU by composing 
the Poincare transformations PajOj - The Poincare holonomy of the loop can therefore 
be calculated as 

— PaiakPctkak-l ' ' ' Pa^ai- (2.47) 

Note that this holonomy only depends on the 2-cell, around which the loop was taken 
and the frame chosen in the 4-ceU where the loop was started. 

We can use piecewise flat manifolds to describe general configurations of defects. 
In such a description the 2-cells correspond to conical defects, the 1-ceUs correspond to 
junctions, and the 0-cells correspond to collisions of defects. 

In a generic piecewise flat manifold any 2-cell will have a non-trivial holonomy 
including spacelike 2-cells (i.e. 2-cells with a Euclidian metric signature). Such a 2-ceU 
would correspond to a spacelike conical defect, which is prohibited by the local causal- 
ity principle of our model. One could attempt to enforce this principle by disallowing 
piecewise flat manifolds with spacelike 2-cells. This however seems overly restrictive. 
Instead we opt to allow spacelike 2-cells only when they have a trivial holonomy. These 
cells are treated as "virtual defects" necessary to facilitate the piecewise flat description. 

The advantage of using a piecewise flat manifold structure to describe a configura- 
tion of defects is that it encodes the geometry in a local cell-by-cell way. Moreover, if two 
configurations share the same description as a piecewise flat manifold they are geomet- 
rically indistinguishable, and should physically be considered the same configuration. 

The downside is that there are many piecewise flat structures that describe the same 
configuration of defects. Not only is there the choice of framing of the 4-cells (which 
can be considered as a kind of local gauge fixing), the process of subdivision allows the 
creation of many equivalent piecewise flat structures. A related issue is that the 4-cells 
are compact which means that we even need an infinite number of cells to describe 
Minkowski space. This last issue can easily be resolved by a slight generalization of the 
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concept of a piecewise flat manifold by allowing half -open cells that extend to infinity. 
In the rest of this thesis we will generally use this generalized concept for simplicity. 

2.8 Other piecewise flat approaches to gravity 

The use of piecewise flat structures in gravity dates back to Tulio Regge's seminal 1961 
paper "General Relativity without Coordinates". |[74) His approach — now known as 
Regge calculus — is based on the observation that any (pseudo)Riemannian manifold 
may be approximated by an appropriately fine piecewise flat simplicial complex. The 
lengths of the 1-simplices, which determine the flat metrics on the other simplices, are 
taken to be the fundamental variables. 

Regge calculus has been successfully used both as an approximation scheme for 
numerically attacking problems in classical general relativity (see |101 1 for a review and 
bibliography), and as the basis for quantum mechanical approaches to gravity (see \75j 
for a review). In many respects the piecewise flat gravity model studied in this thesis 
is similar to Regge calculus. A key difference, however, is that in Regge calculus the 
piecewise flatness of the geometry is a discretization of the classical geometry. As such 
the conical defects on the codimension 2 simplices are not viewed as physical excitations 
of the model. Consequently, no restrictions are enforced on the spacetime signature with 
which the defects may occur. 

Another difference is that quantum gravity approaches based on Regge calculus 
such as the original Ponzano-Regge model \73j and the Turaev-Viro model [96] in 3 
dimensions and the Barrett-Crane model | [T3) in 4 dimensions typically assume an Eu- 
clidean signature. Although the model under consideration here does not currently 
address any issues of quantization, it is inherently Lorentzian in nature. 

The spin foam models considered in the context of loop quantum gravity can be 
viewed as a Lorentzian generalization of the Euclidean quantum gravity models based 
on Regge calculus. Moreover, loop quantum gravity uses the holonomies of spatial 
loops as the fundamental variables, much like the holonomy description of a configu- 
ration of conical defects. This begs the question whether there is some relation with the 
holonomy description of the model studied here. In particular, one could imagine that 
this model could appear as the classical limit of loop quantum gravity. This question 
was addressed by Eugenio Bianchi.|15| He tried to reproduce the kinematical state 
space of loop quantum gravity by standard path integral quantization techniques. His 
conclusion was that restricting the path integral to locally flat metrics was too strong of 
a restriction. Instead to obtain the kinematical state space of loop quantum gravity one 
needs to integrate over the wider class of locally -flat connections. 
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Another piecewise flat approach to quantum gravity based on simplicial complexes 
is causal dynamical triangulations.|5-7| Unlike Regge calculus, causal dynamical tri- 
angulations does not use the edge lengths of the simplices as variables, but rather 
fixes all edge lengths to a constant, and relies on the sum of different triangulations to 
produce dynamics. Even more so than Regge calculus this model views the piecewise 
flat structure as a discretization tool with the theory proper being defined only in the 
limit that the edge lengths go to zero. As such in causal dynamical triangulations 
the appearance of spacelike 2-simplices with a non-trivial holonomy is not viewed as 
problematic and is in fact unavoidable for the model to work. 

In spirit, the world crystal model proposed by Kleinert PH[55) seems to be the most 
similar to the model considered here. He too considers a model of gravity where the fun- 
damental degrees of freedom are curvature defects. In his model however the defects 
are not constrained to be flat planes in spacetime, but are allowed to have an extrinsic 
curvature. As consequence they must be accompanied by a non-zero gravitational field. 
The resulting spacetime is not actually locally flat. In that aspect that model crucially 
differs from our model, where local flatness is a fundamental principle. 
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A new issue in 3+1 dimensions is that we must consider the collision of two defects. 
Generically, two codimension 2 planes have an intersection of codimension 4. Conse- 
quently, conical defects in 2+1 dimensions generically have no intersections. That is, 
in 2+1 dimension we can safely assume that two point particles never collide. In 3+1 
dimensions, however, the world sheets of two line defects generically have one common 
point. This point indicates an event where the two line defects collidej^ 

Consequently, we must add a prescription of what happens when two line defects 
collide to our a model. Because the holonomies carried by the line defects do not com- 
mute, the line defects cannot simply pass through each other. This can be understood 



geometrically if we consider the situation in figure 3.1 which shows the collision of 
two line defects. On the left the deficit angle of the white defect is drawn away from 
the approaching black defect. On the right the same situation is drawn with the deficit 
angle in the opposite direction. The black defect now has a kink after passing through 



the white defect. In section 2.6 we saw that in our piecewise flat model a kink in a defect 
implies that there is at least one other defect meeting at that point. 

Our piecewise flat model must somehow prescribe what happens after such a colli- 
sion. For this it needs to provide an intermediate configuration of conical defects that 
continues the geometry after the collision|^ In this chapter we will discuss how these 
"continuations" can be constructed for general collisions. 



^^Note that for an arbitrary time slice considered as the present this point may lie in the future or in the 
past. 

^^For the sake of brevity of expression we will simply refer to such an intermediate configuration as a 
"continuation". 
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Figure 3.1: The black defect cannot pass through the white defect without incurring a kink. 

3.1 Collision parameters 

We wish to study a general collision of two conical defects A and B. For any collision we 
can find a locally flat neighbourhood that contains no other defects than A and B (prior 
to the collision). For the study of general collisions it is therefore sufficient to study the 
collisions of infinite defects. 

Two defects have 14 free parameters. To determine how many of these can be fixed 
by the choice of Poincare frame, we proceed as follows. Let us use the translational 
freedom to fix the spacetime point of the collision at the origin. This has the advantage 
that the translational parts of the holonomies of both defects vanish and we only have to 
worry about the Lorentz part of the holonomies. We can then use two rotations to fix the 
orientation of the defect B (say along the z-axis), and use two boosts to fix the velocity 
of B (say we make B stationary). We now have only two degrees of freedom in our 
choice of frame, a rotation and a boost in the direction of the orientation of B. We can 
use the boost to make the velocity of A perpendiciilar to the defect B. The remaining 
rotation can then be used to completely fix the direction of that velocity (say along the 
X-axis). 

Thus we find that there are 4 remaining degrees of freedom: the conjugacy classes 
of the two defects given by the size of their deficit angles in their rest frame ij^A and ipB, 
the difference in orientation of the defects given by the angle </> between the two defects 
at the instant of collision, and the relative velocity between the two defects given by the 
difference in rapidity rj. 
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Figure 3.2: The parameters of a general collision oftivo defects: tp — A, ips, v (= tanh rj), 
and (p. 

Note that there are different ways in which the frame can be fixed. In the example 
above we chose the frame where defect B was stationary. We will call this the rest 
frame of B. Similarly we could have chosen the rest frame of A. These frames are often 
convenient when doing calculations with the holonomy. For geometrical calculations, 
it is often useful to make a more symmetrical choice of frame where both defects have 
an equal velocity. We will refer to that frame as the "center of velocity frame". In this 
frame the difference in velocity will not be equal to that in the rest frames, the difference 
in rapidity however is the same in all three frames. 

3.2 Orthogonal collisions 

We first consider a specific class of collisions where the colliding defects are orthogonal 
at the point of impact, i.e. (p = These types of collisions turn out to allow an 

especially simple continuation, yet they already introduce some of the issues that we 
encoimter when constructing continuations for general collisions. 

As defects A and B collide they effectively cut each other in four semi-infinite defects 
(which we label Ai, A2, Bi, and B2 as shown in figure [33| . The task of constructing a 
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Figure 3.3: The loops relating the holonomies of the defects before the collision (on the left) 
to the holonomies of the semi-infinite defects after the collision (on the right). 



continuation for the collision consists in finding a network of finite defects that interpo- 
late between the ends of these semi-infinite defects. 



Figure 3.3 shows how the holonomies of the semi-infinite defects are related to the 
holonomies of the original defects. The loops and 7^13 around the semi-infinite 
defects Ai and A2 after the collision are homotopic to the loops with the same labels 
before the collision. The loop 7^^ before the collision is homotopic to the concatenation 
1~B ■ 7Ai • iB- Consequently, for the holonomies we find 

Q A2 ^ Q~bQ AiQ B ■ (3.1) 

The loops 7bj and 7^3 after the collision are both homotopic to the loop 75 before 
the collision, and therefore have the same holonomy. Consequently, we can rewrite 
equation | |3.1) as 

QA,QB,QXQ'k=^- (3-2) 

This equation expresses that the concatenation of loops 7^1 • 7Si ■ I'a^ ' Ib^ contractible, 
and thereby that the defects formed two separate networks (i.e. two distinct defects) 
before the collision. 

In the rest frame of the original B defect, the defect Ai is spanned by the 4-vectors 



< = (l,z;,0,0), 
d'^ = (0,0,1,0). 



(3.3) 



42 



3.2: Orthogonal collisions 




Figure 3.4: An orthogonal collision allows a continuation consisting of just a single 
intermediate defect connecting the kinks in the defects A and B. The loops in the diagram 
show how the holonomy of this new defect is related to the holonomies of the colliding defects. 



The holonomy of the B defect in this frame is a rotation about the z-axis, Rz(tpB)'^u- 
Consequently, the semi-infinite defect A2 is spanned by the 4-vectors, 

u':^ = RzH^bYu^i = (l,wcosV'B,wsini/i_B,0), 
d'^^Rzii^sTJi = (0,-sin?/>B,costAi3,0). 

Consequently, the semi-infinite defects Ai and A2 share a common junction spanned by 
the 4-vector 

/'=.(l,«,«tan^,0). (3.5) 

A similar junction is shared by the semi-infinite defects Bi and B2. 

This suggests that a simple continuation is possible with a single finite defect AB 
connecting the junctions in the A and B defects. The holonomy of this defect is given 
by the junction condition of the junction in the A defect (see figure 3.4[ p^ 



Qab — QaiQa2 — QaQbQaQb- (3.6) 
^*Or equivalently the junction condition in tlie B defect. 
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We still need to check if the resulting defect is physically acceptable, i.e. whether its 
holonomy is rotationlike. Since the junction is a common line for all defects meeting 
there, if the junction is timelike so are all incident defects. Consequently, a sufficient 
(although not necessary) condition for defect AB to be timelike is that the junction in 
the A defect is timelike. 
The norm of is 

So the junction is spacelike when 



„2 > cos^ (3.8) 
Consequently, for orthogonal collisions with high velocities and/or high deficit angles 



the junctions become superluminal. In section 2.6 we observed that such junctions were 
undesirable for our model. However, both the Ai and A2 defect lie in the past of this 
junction. It is therefore of the least problematic kind, i.e. the fact that this junction is 
tachyonic needs not be seen as a violation of local causality. 

A more pressing worry is that a superluminal junction leaves room for the possibility 
that the intermediate defect AB itself is spacelike. This can be checked by finding the 
conjugacy class of the holonomy Qab of the intermediate defect. If this is rotationlike, 
the intermediate defect is timelike, and if this is boostlike, the intermediate defect is 
spacelike. A direct way to find out to what conjugacy class Qab belongs is to calculate 
the trace of Qab in the PSL(2, (D) representation of the Lorentz group (see appendix [A|. 
In the rest frame of the B defect the PSL(2, (D) representations of Qa and Qb are 

^ / cos ^ — sin ^ sinh rj sin ^ cosh 77 

sin ^ cosh 77 cos ^ + sin %^ sinh n , 

The trace of the holonomy of the intermediate AB defects is consequently given by 

TtQab = T^^QaQsQaQb 

1 — + cos tpA + cos tpB — cos TpA COS TpB (3.10) 

with V = tanhrj. A Lorentz transformation in PSL(2, (D) is rotationlike if and only if its 
trace is real and the absolute values of the trace is smaller than 2. Consequently, the 
intermediate defect AB is rotationlike only when 

^2 ^ _ (1 -COS7/)a)(1 -COsV>b) 
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Figure 3.5: Tlie regions of the parameter space of orthogonal collisions of defects with equal 
mass ip, where the junctions are timelike (in dark gray) and where the intermediate defect 
itself is timelike (In light gray). 



Since the right hand side of this inequality is smaller than one for non-zero values of ipA 
and ijjB, there are high velocity orthogonal collisions for which the intermediate defect 
becomes spacelike. 

We can compare the region of the collision parameter space where subluminal junc- 
tions become impossible to the region where the intermediate defect becomes spacelike 
— as is done for a collision of equal mass defects in figure 3.5 As we argued, the former 
is contained within the latter; The junctions being timelike is a sufficient condition for 
the intermediate defect to be spacelike. The region where the intermediate defect is 
timelike is significantly larger, showing that requiring the junctions to be timelike is not 
a necessary condition for the intermediate defect to be timelike. The region however 
still is bounded. 

Apparently, there is a large region of the parameter space of collisions where the 
single defect continuation is not consistent within our model. In particular, collisions 
with a high velocity cannot be continued in this way. Now, as we will see in the 
next section, the single defect continuation is not even an option for non-orthogonal 
collisions. We therefore need to consider more complicated continuations involving 
more intermediate defects. These continuations may also help us deal with high velocity 
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Figure 3.6: A slanted collision. The semi-infinite defect lines Ai 
after the collision. (Compare to figure i 



A2 do not meet up 



3.1 



orthogonal collisions. 



3.3 Non-orthogonal collisions 

When the defects A and B are not orthogonal at the collision point (i.e. (p ^ 7r/2), 
then the semi-infinite defect lines A\ and A2 do not share a common junction after the 
collision (see figure [3^ . Consequently, a continuation with just one intermediate defect 
is not possible. We must therefore consider more complicated continuations involving 
multiple intermediate defects. 

Another way to see that a continuation with a single defect line is impossible is to 
calculate the holonomy, Q^^, of a loop transverse to the collision (see figure [3T| . This 
is the holonomy of a defect that would connect the colliding defects. In the previous 
section we found that it is given by 



Qt — QAxQai — QaQbQaQb- 



(3.12) 



For non-orthogonal Qa and Qb the trace of this holonomy is a complex number (we 
will spare the reader the actual expression). Consequently, it cannot be the holonomy of 
a single defect line. 

This result is, in fact, slightly stronger than the previous result. If the only problem 
was that the defects Ai and A2 did not meet up in a single junction we could have tried 
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Figure 3.7: A loop jt transverse to the collision point. 




Figure 3.8: Continuations like this with triangle resolved junctions are not allowed for 
non-orthogonal collision because the middle defect would still need to have the loxodromic 
holonomy Qt- 
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Figure 3.9: There are three different types of quadrangle configurations that provide a 
continuation of a collision. 



to resolve the issue by replacing the junctions with triangle loops such as depicted in fig- 
ure 3.8 However, for such a continuation the single defect connecting the two triangles 



would still need to have the holonomy of the transverse loop 7t- Such continuations are 
therefore also impossible. We need continuations that involve at least two intermediate 
defects running between the two original defects A and B. 



3.4 Quadrangle continuations 

The simplest possible continuation involving at least 2 defects connecting the original 
defects A and S is a square loop of four intermediate defects connecting the four ends of 
the semi-infinite defects Ai, A2, Bi, and B2- There are four topologically distinct loops 



of this type possible as shown in figure 3.9 



We choose loops around the external semi-infinite defects in such a way that equa- 
tion ( |3.2) relating their holonomies becomes 

Qb,Qa,Qb,Qa,^1- (3.13) 



This can be achieved by taking the loops as defined in fig 3.4 and reversing the directions 
of the and loops. This choice is convenient in the sense that all loops are in the 
positive direction with respect to the outwards direction of each external semi-infinite 
defect. 

A set of four intermediate defects introduces 4x7 = 28 new degrees of freedom in the 
form of four holonomies. We are however not completely free to choose the parameters 
of the internal defects. First of all, all internal defects must pass through the collision 
point at the time of collision. This fixes the positions of the internal defects. Since we 
choose the collision point as the origin of our coordinate system, this means that the 
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holonomies of the intermediate defects maybe expressed as pure Lorentz transforma- 
tions. This reduces the number of free parameters from 28 to 4 x 5 = 20. 

The internal defects must also satisfy the junction conditions at each of the four junc- 



tions. For example, for a type III configuration (see figure 3.9 1 the junction conditions 
become (for a suitably chosen set of loops), 

QaiBi = QAiQAiB2 Qa-2Bi^QbiQAiBi 

QA2B2 = QA2QA2B1 QA1B2 = Qb2QA2B2- 

Each of these equations imposes 6 algebraic conditions on the parameters of the internal 
defects. The equations are however not independent. The condition | 3.13[ implies 



that if any three of the above equations is satisfied, the remaining fourth equation is 
automatically satisfied as well. The jimction conditions therefore provide 3x6 = 18 in- 
dependent algebraic conditions on the parameters of the internal defects. Consequently, 
we expect that for each type of quadrangle continuation, there is a two-dimensional set 
of configurations that satisfies all conditions. 

In this section we will go through the steps necessary to find these configurations in 
terms of two parameters, the rapidity with which two of the junctions move along the 
external defect they connect to. As an explicit example we will take a type III configura- 
tion, because as it turns out the answer takes its simplest form in this case. Type I and II 
configurations can be solved by following the same steps, but the intermediate results 
are more involved and will not be shown here. 

To do the calculations it is convenient to represent the holonomies by PSL(2, (D) 
matrices as is explained in appendix |A] This has a number of advantages. First, a 
2x2 complex valued matrix has only 8 real parameters, as compared to the 16 real 
parameters of a 4 x 4 real valued matrix. This helps reduce the number of redundant 
equations. Second, the condition the holonomy of the defect lines needs to be rotation- 
like takes a particular easy form in terms of the trace. 

As is discussed in appendix [A] a necessary and sufficient condition for a PSL(2, (D) 
matrix, Q, to correspond to a rotationlike Lorentz transformation is 

Im Tr g = 

(3.15) 

|TrQ|<2. 

So, in addition to the junction conditions | |3.14^ , the holonomies of the intermediate 
defects are also subject to the conditions 

ImTrQAiBi = ImTrQAgSi =0 

(3.16) 

IniTrQAa-Ba ImTrQAiSa = 0' 
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and 

\TtQa,bJ<2 |TrQ^,Bj<2. 

We can parametrize the space of configurations by taking one of the intermediate 
holonomies {QaiBi in a type III configuration) to be a general GL(2, (D) matrix, 

Qmb. = (^'-'1 '^'') (3.18) 
\g + ih e + ij.J 



The junction conditions | 3.14^ can then be used to obtain the holonomies of the other 



intermediate defects. Since the holonomies of the external defects all have determinant 
equal to one, the junction conditions ensure that if we require that det{QAiBi ) = 1/ the 
determinants of the other intermediate holonomies are one as well. 



The conditions | 3.16| , provide four linear constraints on the eight real parameters. 



The inequalities | |3.17^ then constrain the resulting four-dimensional parameter space 
to a h5^ercube. The complex condition det{QAiBi) = 1 provides two further real 
constraints. These, however are quadratic in the parameters. Therefore, the question 
whether these are compatible with the linear inequalities | |3.17| is somewhat delicate. 

The question can be simplified by taking a slightly different route. First we go to the 
frame where the external defect Bi is stationary and oriented along the z-axis In this 
frame. 



Qb.-[ .-W.j- (3-19) 

The equation ImTr QaiBi = implies that / = —5. The combination of the junction 
condition Q^aSi = QbiQaiBi with the reality condition on the trace of Q^aSi then 
implies that 

= ImTr (353(3^2-82 
= Im((E*'''«/2(a + ib) + (B-''^«/2(g _ j^)) (3_20) 

= (o — e) sin -03/2, 

from which we can infer that e = a. The condition that the determinant of QaiBi is real 
then implies that, 

f a + ib ±(B-'^(c + id)\ 
\— (E^(c — id) a — lb ) 

^^The requirement that B\ is stationary does not uniquely define a frame, since there is the liberty to boost 
the frame in the z-direction. Here and onwards, when we say "the frame were defect X is stationary", 
we mean the frame that can be reached from the center of velocity frame by a boost along the propagation 
direction of the incident defects. 
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In the case that the ± is a +, we can conjugate QaiBi with a boost in the z-direction 
of rapidity /i to obtain a unitary matrix. That is, in the z-boosted frame both Qb^ and 
QaiBi are pure rotations. Therefore in this frame both Bi and AiBi are stationary 
defects, consequently the junction of these defects is also stationary. The parameter /x is 
thereby identified as the rapidity of the Bi junction along the Bi defect, which we label 

In the other case where ± is a — , it is impossible to obtain a unitary matrix from 
QaiBi by boosting it in the z-direction. It is therefore impossible to find a frame where 
the junction on the Bi defect is stationary. The junction must consequently move at a 
super luminal velocity. Note that this does not necessarily mean that the defects AiBi 
and A2B1 move at superluminal velocities. 

We thus see that requiring that the junction moves at a subluminal velocity is a 
sufficient (although not necessary) condition for the incident holonomies to satisfy | |3.16) 
and | 3.17| . By requiring that QaiBi takes the form. 



where the bar on the Bi indicates that the holonomy is observed in the frame where 
the Bi defect is stationary, we can impose the conditions 13.161 and ( 3.17| for the two 



incident holonomies QaiBi and QA2B1, and the requirement that determinant of the 
intermediate holonomies is real. If we change to a frame where the external defect 
diametrically opposite Bi {B2 in the case of a type III configuration) is stationary, we 
can require that one of the incident intermediate holomies QaiB^ or Qa^b-, is of the 
form 13.22 1 to enforce the other half of the conditions ( 3.16| and 1 3.17) . This leaves only 



one equation, the requirement that the determinant is one. For a matrix of the form 
1 3.22) , this takes the simple form. 



+ 52 + 2^^1. (3.23) 



Enforcing all conditions will leave 2 free parameters. A logical choice would be the 
rapidities of the two diametrically opposite junctions, fiBi and 1132, since we already 
needed to require that these junctions move at subluminal velocities to ensure that the 
inequalities | |3.17[ | are met. However, as it turns out these parameters are not indepen- 
dent; fixing one automatically fixes the other. As far as this author has been able to 
find, there is no simple argument that explains why this should be, other than explicitly 
doing the calculations. This is what we will now do for a configuration of type III. 

To do the calculations we must first choose frames for the external defects. What 
sets configurations of type III apart from the other two types is that the diametrically 
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opposite junctions Bi and B2 lie on the two ends of what was a single defect before the 
collision. This allows us to choose a frame in which both ends are stationary simultane- 
ously. 

The holonomies of the external defects with this choice are 



Qbi = Rziips) = 



iVB/2 







<si^B/2 
(eW'-b/2 q 



<s-'^-4'b/2 
B,irj)R,{cb)R,i^A)R.W-'B,irjy' 

f cos ^ — i cosh(77 + 10) sin ^ 
\ — i sinh(77 + 10) sin ^ 



cos 



i sinh(r7 + i0) sin ^ 
^ + i cosh(77 + 10) sin - 



(3.24) 



Qa, = R,{il^B)B.,{ri)R.{4>)Rz{i^A)RA4>)-^B,{riy^R,{^BV 

( cos %^ —i cosh(?7 + 10) sin %^ — i (e""''-^ sinh(77 + i0) sin ^ 
\ i (e'"^'-^ sinh(ry + i0) sin %^ cos ^ + i cosh(r/ + i0) sin %^ 

One can readily check that these holonomies indeed satisfy ( |3.13) . We further choose 



a + 16 



"c(e''<:) 
-16 



(3.25) 



such that conditions | 3.16[ | and 1 3.17) are automatically satisfied for Qa^Bx ^^'i Q^aSi- 
We can then obtain QaiBs by using ( 3.14| ; 



(3.26) 



We now need to check that Qa^ b^ has the form 1 3.22 1 in the frame where B2 is stationary. 
Normally this would require a conjugation with an appropriate Lorentz transformation, 
but in the special frame we choose both B\ and are stationary, so 



Q^iSa — Qa^QaxBx- 



(3.27) 



Since dot Qax = 1/ the determinant of Qa^Bo is automatically real. Therefore, to ensure 
that QyijBj has the form 1 3.22) it is enough to require that 



(3.28) 



where j] denotes the entry of the matrix Q at position (i, j). This complex-valued 
equation, can be written as two linear real-valued equations for the real and imaginary 



52 



3.4: Quadrangle continuations 



part. These equations can be solved for a and b, 

sinh/^Bi . cosh/iBi 

a — c(cosx ; Sinx ; — j 

tann rj tan cp 

sinh/XBi . cosh^Bi 

6 = c (cos X — — r h sm X — — ) ■ 

tann 1] tan (p 

Substituting this in the equation for the determinant of Qa^Bi, 
det Qa^Bi ^a^ + b^ + c^ 



(3.29) 



cosh^/isi sinh^usi (3.30) 

c (1+ ,_2^ + — n— 



tanh r; 

we find that setting the determinant to one fixes c to be 

1 



-j^ I cosh /j-Bi I sinh 



(3.31) 



We now have two free parameters fiBi and x- We would like to replace x with /xb2- We 
can extract fiB2 from QaiB^ taking. 



Qa,bJ2,1]* 
(A + /i)mBj + h{hA - 1) 
(A - h)mBi - h{hA + 1) ' 



(3.32) 



where 



ipA smhrj 

A — tan — , n ~ , and tob, = tann/isi. (3.33) 

2 sin 



The dependence on x has completely dropped out; is completely determined by 
fiBi arid the collision parameters rj, (f), and We can therefore not replace the free 
parameter x with fiB2- If we want the second free parameter to be a junction rapidity, 
we will have to use the rapidity of the Ai or A2 junctions. 

The rapidity fiAi may be obtained by transforming QaiBi to the frame where Ai is 
stationary. From the assignment of the external holonomies 1 3.24^ we see that this may 
be achieved by conjugating with Bx{'q)Rx{4'), 

QMB,^R.{<l^r^BM-'QMB,BMR.{<jy). (3.34) 

Extracting ^iai and solving for x we find, 

tanx = (3.35) 
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where niAi = tanh ^^i- The same steps may be repeated to obtain /i^2 from QA2B1' ^^id 
we find 

h - BmA2 .0 
tanx=^7— , (3.36) 



where mA2 = tanh fiA2- Combrnrng expressions 1335} and 13.36 1 we find 

TOAi'TlAa - 4("^Ai + m-Aa) - = 0. (3.37) 



If we introduce 17132 = tanh/LtB2/ then equation 1 3.32) can be rewritten to take a similar 
form, 

mBimB2 — ^("^Si + n^Ba) - = 0. (3.38) 

The space of quadrangle continuations of type 111 may be parametrized by two numbers, 
one of the pair {ttia-^ , toa, ) arid another of the pair (m^i , ). Physically, these pairs of 
numbers represent the velocities of opposing junctions on the quadrangle in the center 
of velocity frame. If one of the velocities of a pair is known, then the other is given as a 
function of the collision parameters 77, 0, tpA and ipB- 

Similar relations between the velocities of opposing junctions may be found for 
configurations of type 1 and 11. They may be obtained by following the same steps as 
outlined above for a type 111 configuration. The intermediate results are somewhat 
lengthy so we just give the final relations. For a configuration of type 1 one finds, 

/hA + B A + hB , 2 1 - hAB 

ruA, ms, + h- niA^ + h- m^, — h — = (3.39) 

h-AB h-AB - h-AB 

'hA + B A + hB ,^1-hAB ,^ 

toaoTOb, + h- m^o + h- niB, — h — = 0, (3.40) 

h-AB h-AB h-AB 

and for a t5^e 11 configuration, 

/hA-B A-hB ,^1 + hAB ,^ 

'that'itib, + h- iriA-, — h- ms, + h — = (3.41) 

h + AB h + AB h + AB 

/hA-B A-hB .A + hAB ^ ,^ 

mA2^B2 + I''' "^Ao ~ IT'' 71132 + " ~ — (3.42) 

h + AB h + AB h + AB 

where in both sets of equations we used the shorthand, h = cosh{'ri)/ cos((/)). 

Comparing the set of relations for each configuration, we notice a couple of things. 
First, all configuration types share the property that the velocities of opposing junctions 
on the quadrangle are related only to each other, and are independent of the velocities 
of the junctions on the other diagonal. This decoupling was already noticed by 't Hooft 
in [ |92J based on numerical calculations, and here is confirmed analytically. 
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Second, as we already noted, configurations of type III are a bit special. Equation 
1 3.37) relating the velocities of the Ai and A2 junctions, is completely independent of 
tpA, while equation 1 3.38 1 relating the velocities of the Bi and B2 junctions, is completely 
independent of i/'s • This property makes equations ( 3.37| and 1 3.38) take a very simple 
form, and is undoubtedly related to the fact that in a configuration of type III opposing 
junctions in the quadrangle lie on ends of what was originally one defect. 

Third, in each set of equations — 0.37\ and | |3.38| , | |3.39| l and | |3.40| |, and | |3.41) and 
< 3.42) — there is a simple substitution rule that relates one equation to the other. For 
example, if in equation | |3.41) we replace with m^^, niBi with rns^, we obtain 
equation ( |3.42| . 

Finally — and most disturbingly — for some values of the collision parameters the 
equations do not allow solutions where the velocities of the junctions are all subluminal 
(i.e. where mf < 1 for all junctions). This can be seen by looking at the dependence on 
h. In each equation the h"^ term is independent of the m/s. Since the m/s are bounded, 
the lower order h terms, which do depend on the m/s, cannot cancel the term for 
large values of h. 

Consequently, for large enough values of h (in relation to the other collision param- 
eters) there are no quadrangle solutions of any type with only subluminal junctions. 
Moreover, unlike the superluminal junctions in the one defect continuation of an or- 
thogonal collision, the superluminal junctions of the quadrangle continuation have two 
defects in their future (the new intermediate defects) and only one in its past (of the 
external defects.). In a frame where such a superluminal junction is contained in a 
single time slice, it will appear as if the colliding defect instantaneously splits in two at 
the moment of collision, even at far away points of the defect whose causal past did not 
include any hint of the impending collision. This is totally at odds with local causality 
principle at the base of our model. 

Recall that h was defined as . A large value of h is therefore associated with a 
large value of or a value of (p that is close to zero. So, as was the case for the one defect 
continuation of an orthogonal collision, it are the high velocity collisions that cannot be 
adequately continued by a quadrangle configuration of defects. In addition, we find 
that collisions where the orientations of the colliding defects are close to collinear, are 
also problematic. 

In the case of orthogonal collisions the quadrangle continuations are however an 



improvement over the single defect continuation we discussed in section 3.2 Figure 



3.10 plots the regions of the collision parameter space for orthogonal collisions of defects 



with equal mass, that can be continued by the three types of quadrangle continuation 
or the single defect continuation. We see that quadrangle continuations of t5rpe I and 
II offer no improvement over the single defect continuation for orthogonal collisions. 
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Mass density ij/ 



Figure 3.10: The region of the parameter space of orthogonal collisions of defects with equal 
mass ip, for which a continuation can be found consisting of quadrangle or single defect 
configurations with only subluniinal junctions. 



Quadrangle continuations of the III type with only subliuninal jiinctions, however, can 
be found for a much larger range of collisions. 

The appearance of superluminal junctions of the opening type is a serious challenge 
to the consistency of the principles of our piecewise flat model of gravity. It seems 
that the local causality principle should disallow this type of junction, as it disallows 
superluminal defects. We may still wonder however if allowing superluminal junctions 
would allow quadrangle continuations for all collisions. 

If we do not require all jimctions to be timelike, we are no longer guaranteed that 
the intermediate defects are timelike as weE. In the discussion above we saw that for a 

general collision we can only fix the velocities of two junctions that are adjacent on the 
quadrangle. If we fix these to be subluminal, then the three intermediate defects that are 
incident to these two junctions are timelike. For certain high velocity or nearly coUinear 
collisions, the remaining two junctions will both be superluminal. Consequently there 
is no guarantee that the intermediate defect connecting these two junctions, is timelike. 
However, a finite defect connecting two spacelike junctions may still itself be timelike. 

We want to know if there are collisions for which at least one of the defects of 
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Figure 3.11: Tlie region of the parameter space of orthogonal collisions of defects with equal 
mass 4', for which a continuation can he found consisting of a quadrangle or single defect 
configuration with only non-spacelike defects. 



any quadrangle continuation is spacelike. It is to be expected that the "best" result is 
obtained when all junctions are allowed to become spacelike. The way to check whether 
a defect is timelike is to check whether the absolute value of the trace of its holonomy 
(in the PSL(2, (D) representation) is smaller than 2. 

Consequently, to test if a continuation with only timelike defects for a certain col- 
lision is possible, we need to calculate the holonomies of the intermediate defects and 
check if there exist values in the two dimensional parameter space of quadrangle con- 
tinuations (of a fixed type) where the absolute value of each trace is smaller than two. 
This problem does not seem to be tractable analytically. It can however be checked 



numerically. Figure 3.11 shows the result for orthogonal collisions of defects with equal 
mass density. 

Compared to the single defect configuration the quadrangle configurations allow 
continuation of a significantly larger portion of the parameter space of collisions. The 
quadrangle continuations of type III seem most effective in finding continuations for 
orthogonal collisions with a high relative velocityrj However, for any non-zero value 



^Type I and II quadrangle configurations are more suited for continuing slanted collisions. For collisions 
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of the mass density of the defects, there exist velocities for which none of the types of 
quadrangle continuations has only timelike defects. 

So, even if we allow superluminal junctions, quadrangle configurations cannot be 
used to continue all collisions within our model. We do however observe that for or- 
thogonal collisions, quadrangle continuations provide a substantial improvement over 
the single defect continuation considered in section 3.2 This suggests that similar im- 
provements may be made by considering more complicated continuations. It could be 
that one can find a continuation for any collision by considering an arbitrarily compli- 
cated intermediate configurations. To test this we consider a more complicated type of 
continuation in the next section. 



3.5 Tetrahedral continuations 

The quadrangle configurations considered in the previous section basically consisted of 
four vertices moving away from the collision along the external semi-infinite defects, 
which were connected by four internal finite defects. This meant we needed to choose 
which vertices should be connected by finite defects with the three possible distinct 
choices leading to the three possible types of quadrangle continuation. We will now 
consider the more general case where each vertex is connected to each other vertex. The 



result is a tetrahedral configuration of finite defects as shown in figure 3.12 

The quadrangle continuations of type I, II, and III can be considered to be special 
cases of this configuration where the deficit angles of two of the internal finite defects 
become zero. For example, a quadrangle continuation of type III is obtained by setting 
the deficit angles of the defects labelled A1A2 and B1B2 in figure 3.12 to zero. 

The six new internal defects originating from the collision point give 6-5 — 30 
degrees of freedomj^ The algebraic junction conditions | |2.46^ provide six conditions 
for each vertex. These are however not all independent. Since the holonomies of the 
external semi-infinite defects are related through | |3.2| the junction conditions for the 
four junctions give 18 independent conditions on the internal defects. The geometric 
condition that all defects incident to junction must share that junction as a common 
line, provides a further 2 conditions per 4-valent vertex. Consequently, we expect that a 
tetrahedral configuration of internal defects has four free parameters. 

This agrees with our above observation that we can obtain the quadrangle configu- 
rations (each having two free parameters) by setting two deficit angles to zero. In the 

that are close to coUinear these types will be most effective. 

^^Since the defects originate from collision their positional degrees of freedom are fixed by the fact that they 
must pass through the collision point. 
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Figure 3.12: Connecting all the vertices yields a tetrahedral configuration of internal 
strings. 

case of the quadrangle configurations we took two of the velocities of the vertices along 
the external defects as our free parameters. Since we have four vertices moving along 
four external defects and four free parameters it is tempting to take these velocities (or 
their corresponding rapidities) as our free parameters. This choice makes it rather easy 
to explicitly impose that all vertices are subluminal. 



3.5.1 Newtonian limit 

Solving the system of 26 conditions on the 30 degrees of freedom of the internal defects 
is rather complex. It is therefore instructive to first solve these conditions in the Newto- 
nian limit, i.e. the limit that the velocities of all defects are much smaller than the speed 
of light and all deficit angles (linear mass densities) are much smaller than one. 

For low velocities a Lorentz transformation can be approximated by a Galilean trans- 
formation. That is, represented as a 4 x 4 matrix, a Lorentz transformation Q may be 
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written, 

where iT s R'^ is the boost velocity and i?(0) G 50(3) is a rotation of |0| degrees around 
the (j) axis. So the holonomy Q{(j),v) of a line defect with orientation and deficit angle 
11011 moving with velocity v is given by 

0(*..^)»r' "Ml J!.) (I 1) (3.44) 



(3.45) 
(3.46) 



-V IJ \Q R{(j>) I \v I 
' 1 ^ 

' ' 'A 

-cj)xv R{(j)) / ' 



where in the last line we used that for small \\<f)\\, R{(j)) ■ x ^1 + <j) x x. 
The junction conditions 1 2.46) therefore become 



I = Qi • • • Q„ (3.47) 

y-(f)i X vi i?(0i) J y-(f>n X ■;;„ i?(0„) y 

\ -01 X wi -...-(/)„ X H h 0,1) / 



(3.48) 
(3.49) 



Consequently, in the Newtonian limit the junction conditions become 

01 H h 0„ = 0, and (3.50) 



X ui + . . . + 0„ X w„ = 0. (3.51) 



Now consider a tetrahedral configuration of defects as in figure 3.12 and let m; be 
the velocity of the vertex i along external defect i with orientation 0^ and velocity Vi. 
Since the velocities are non-relativistic the total velocity of the junction i, Wi, is then 
given by 

Wi = Vi + rrii . (3.52) 
\Wi\\ 
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From the velocities of the junctions we can obtain the orientation <j>i j and velocities Vij 
of the defect connecting junction i and j as follows. 



Wi — Wj 
\\Wi - Wj 
Wi + Wj 



and (3.53) 

(3.54) 



We have now solved 24 of the 30 degrees of freedom of six internal defects in terms of 
the velocities of the junctions along the external defects m,. The remaining unknown 



parameters are the deficit angles of the internal defects. Equation 3.50 tells us 
that we can find the aij by decomposing the vector — 0j on the base formed by the 
orientations ipij of the three internal defects connected to the external defect i. Equation 
( 3.51) is automatically satisfied as well. 



Since we independently obtain values for j and aj,i it may appear that the answer 
is overdetermrned. However, up till now we have ignored the relation ( 3.13) for the 



external defects. These give two relations among the external defect parameters 

4>Ai + (l>A2 + 4>Bi + 4>B2 = 0, and (3.55) 

(l>Ai X + 'f)A2 X ^A2 + $Bi X Vb^ + $82 X VB2 = 0. (3.56) 

It is a straightforward exercise in vector algebra to show these conditions imply that 
Q!i j = "Qfj.i- The deficit angle of each internal defect is uniquely determined. 

We find that — in the Newtonian limit — fixing the velocities of the junctions along 
the external defects indeed fixes all the internal parameters of a tetrahedral configu- 
ration. We can thus use these velocities as the four free parameters of the tetrahedral 
continuation. It is expected that this is true as well in the exact case, at least for small 
values of if) a, ''Pb, and 77. For larger values of these parameters one may expect that not 
all choices for the velocities of the junctions yield valid configurations. Ultimately, it 
may even happen that for some values of the collision parameters there are no valid 
configurations for sublumrnal choices of these velocities, as we will see. 



3.5.2 Exact piecewise flat continuations 

In the preceding sections we have studied the single defect and quadrangle contin- 
uations using the algebraic relations given by the holonomies of the defects. In the 
case of a tetrahedral continuation the algebraic relations between the 30 (or 36 if you 
start out with general Lorentz transformations) free parameters become very complex. 
Already in the case of the quadrangle continuations, solving the relations required 
significant assistance from computer algebra software like Mathematica. Despite some 
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serious attempts to directly solve these algebraic relations, that approach has proven 
impractical. 

The source of the complications here (in part) lies in the fact that the junction condi- 
tions at one jimction influence the conditions at other junctions in a non-trivial way. The 
non-local nature of the algebraic description of the defects in terms of their holonomy 
really starts to become an hindrance in solving the relations. In section [Z7| we formu- 
lated an alternative way of describing a configuration using the piecewise flat geometry 
that it generated. The advantage of that approach was that it encoded the geometry 
in an essentially local way. In this section we will use this approach to find a general 
tetrahedral continuation of a collision of defects. 



The general idea will be to split the geometry in five 4-cells, as shown in figure 3.13 
One 4-cell outside of each face of the tetrahedron, and a fifth 4-cell as the inside of the 
tetrahedron. 

Far away from the collision (well outside its lightcone) the geometry of the spacetime 
will be unaffected by the collision and will simply evolve linearly from the geometry 
before the collision. If we cut out the inner region containing the collision and its 
continuation, then this outer region contains four semi-infinite defects, each a half of 
one of the colliding defects. We will use these defects to divide the outer region in four 
(half open) 4-cells. 

The semi-infinite external defects will become (half open) 2-cells in the boundary 
of these outer 4-cells. For this we first need to choose an inner boundary for these 2- 
cells. For each external defect we choose a line contained in the defect passing through 
the collision point. This line will become the 1-cell boundary of the half-open 2-cell 
formed by the external defect. It will also serve as one of the junctions of the tetrahedral 
configuration. 

Each choice of line corresponds to one free parameter given by the velocity of the 
junction moving outwards along the external defect. The lines need to be chosen to 
move away from the junction fast enough such that the geometry in their vicinity can 
be obtained as a linear continuation from the geometry before the collision]^ 

Now, to construct one of the outer 4-cells we choose three external defects. Since 
the junctions on these defects all pass through the collision point, they lie in a common 
hyperplane. The convex hull of the junctions in that hjrperplane will form the 3-cell at 
the inner boundary of the 4-cell. The boundary of this 3-cell is formed by three 2-cells, 
each spanned by two of the junctions. These 2-cells will become the new intermediate 
defects in the tetrahedral configuration. 

i^We will later discuss what restrictions this puts on the construction. 



62 



3.5: Tetrahedral continuations 




Figure 3.13: The tetrahedral configuration can be divided in five A-cells. Here shown on 

a time slice after the collision, such that each A-cell is represented by a 3-dimensional cell. 
There are four external cells extending to infinity, and one internal tetrahedral cell. Surfaces 
(representing 3-cells) with matching colour are mapped into each other through a Lorentz 
transformation. The black lines on the external surfaces are virtual defects that subdivide the 
external surfaces and carry no holonomy. These are necessary because the external defects 
generically do not lie in the same plane. 
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A pair of external semi-infinite defects does not, in general, lie in one hyperplane. 
We can therefore not connect each pair of external 2-cells by a single 3-cell. Instead 
we need to choose a number of extra 2-cells starting at each junction, such that we can 
define 3-cells that form the boundary between the external 4-cells, as shown in figure 
I3l3l 

The spacetime geometry of the 4-cell and its boundary can be found by linearly 
continuing the geometry from before the collision. Doing this for each set of three 
external defects, we obtain the four external 4-cells complete with a flat metric on each 
boundary cell. Consequently, we have the geometry of the entire 3-skeleton of the 
tetrahedral configuration. 

The missing piece to be added is a 4-cell with a piecewise flat metric, which is to be 
attached to the boundary formed by the four interior boundary 3-cells of the external 4- 
cells. The easiest way to proceed is to observe that if we choose a common constant time 
slicing on the exterior cells — this can be done, for example, by continuing a constant 
time slicing from the center of velocity frame — then the geometry of the boundaries 
of the internal 4-cell on any two time slices are congruent to each other up to a scaling 
factor depending linearly on time. It is therefore enough to find the geometry of one 
time slice of the internal 4-cell. 

The problem of filling a tetrahedron with an internal flat metric given a flat metric on 
the boundary is well-known in Euclidean geometry. It is known that this is possible if 
and only if the boundary metrics satisfy the (generalized) triangle inequalities.] 78 1 That 
is, if and only if the sum of the areas of any three of the boundary triangles is larger than 
the area of the remaining triangle. 

A priori, it is not clear that these triangle inequalities will be satisfied for any choice 
of the collision parameters and junction velocities. However, even if it is not satisfied, it 
is still rather easy to construct a piecewise flat metric on the interior of the tetrahedron, 
as we will show presently. 

To construct an internal piecewise flat metric in the tetrahedron, we will subdivide 
it in four pieces (see figure 3.14| . We build a cell complex (in fact a simplicial complex) 
on the interior in the following way. We add a single vertex, which we connect to each 
of the four vertices on the boundary by 1-cells. Each of the triangles formed by two of 
the new internal 1-cells and one of the 1-cells on the boundary is filled by a new 2-cell, 
and each of the tertahedra formed by three of the new 2-cells and one of the boundary 
2-cells is filled by a new 3-cell. 

For each of the new cells we will have to specify a flat metric that is compatible 
with the metrics on its boundary. Any metric on a 1-cell is flat and is specified by a 
single parameter; its length. Consequently, the four new internal 1-cells give us four 
free parameters. To construct a flat metric on a 2-cell we need that the metrics on its 
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Figure 3.14: A tetrahedron can he subdivided into four smaller tetrahedra hy adding a single 
vertex, four 1-cells, six 2-ceUs, and four 3-cells. 

boundary 1-cells satisfy the (2D) triangle inequalities. This can easily be satisfied for all 
the new 2-cells. One particular way is to choose all four lengths equal and larger than 
the length of the longest 1-cell on the boundar)^^ 




To construct a flat metric on a 3-cell, we again need the metrics on its boundary 2- 
cells to satisfy the {3D) triangle inequalities. The choice above guarantees that three of 
the four inequalities are satisfied. 

To see this, call the length of the new internal 1-cells a and the length of the three 
1-cells on the outer boundary of any particular internal 3-cell bi, 62, and 63. Because 
they are on the boundary of a given triangle they satisfy 



with i,j,k e {1,2, 3}, and without loss of generality we can take the labelling such that 
bi < b2 < bj. The area Ai of the internal 2-cell incident to the 1-cell with length bi is thus 
equal to 



I'This obviously is not tlie most conservative ctioice one could make, but it is convenient for the rest of the 
argument. 




bi + bj > bk, 



(3.57) 




(3.58) 
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Since a > bi for all i, 61 < ^2 < ^3 implies that Ai < A2 < A^. Furthermore we have 

Ai+A2^ \b,^a^ - \bl + \b2^a^ - \bl (3.59) 



>\{b, + h2)Ja^-\hl (3.60) 



> ^3- (3.62) 



As a result any sum involving two of the areas of the internal 2-cells will be larger than 
the area of the remaining 2-cell. 

The remaining triangle inequality is that the sum of the areas of the internal 2-cells 
is larger than the area of the boundary 2-cell. Since the areas of the internal 2-cells can 
be made arbitrarily large by increasing a, it is possible to also satisfy this fourth triangle 
inequality for the new 3-cellsp°| 

Hence we can construct a piecewise flat metric on the interior of the tetrahedron. 
Consequently, the geometry of the missing 4-cell can be taken to be such a piecewise flat 
metric expanding linearly with time. That is, we can expand our original cell complex 
by adding four new internal 4-cells each with the geometry of a tetrahedron growing 
linearly with time (each corresponding to one of the internal 3-cells in our construction 
of the piecewise flat metric). 

The new 1-cells in the construction of the piecewise flat metric, correspond to four 
new internal 2-cells that, generically, carry non-trivial holonomy. These must therefore 
be interpreted as new internal defects. 

3.5.3 Limitations 

In the last section we have shown how to construct a tetrahedral continuation of a colli- 
sion of defects using (piecewise flat) geometrical techniques. The construction involved 
the choice of the outward trajectory of the four junctions, providing four new free 
parameters as expected The construction, however, does not automatically guarantee 
that the junctions and the new internal defects are all timelike. 

The velocities of the junctions are in principle free parameters. If we choose all 
junction to be timelike, then, since the junctions are the boundaries of the internal 



^'^A sufficient condition for this is that a is larger than the radius of the circumscribed circle of each of the 
boundary triangles. 

And possibly four additional parameters for the piecewise flat metric inside the tetrahedron. 
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defects, the internal defects are timelike as well. However, our experience with the 
quadrangle continuations has shown that such a choice may not always be available. 
We therefore need to examine what the available range of possible velocities for the 
jimctions is. 

The first issue that we need to consider is that due to the non-euclidean geometry 
induced by the curvature of the colliding defects it may occur that for a certain choice of 
the velocities of the junctions, there may not be a plane connecting a pair of the junctions. 
This typically happens when the curvature of the defects is negative. 



For example, figure 3.15 shows the 2-dimensional projection of the scattering off a 
defect of a stationary defect with surplus angle equal to tt. After the collision, it is 
impossible to connect the jimctions labelled Ai and A2 with a straight line, making it 
impossible to construct a tetrahedral configuration. In the chosen example this is true 
for any outward moving choice of the junction velocities {m > 0). Junctions that move 
inward (m < 0) can still be connected in that case. If we choose the surplus angle even 
bigger, say 2tt, then it is impossible to connect the junctions for any choice of the junction 
velocities. 

Consequently, we see that the collision parameters put limitations on the possible 
choices of the junction velocities. In extreme situations there may be no possible choices 
at all. In what cases is it possible to find junctions that move at subluminal speeds? A 
sufficient condition for such a choice to be available can be obtained in the center of 



velocity frame. In figure 3.16 the geometry of a collision of two defects with positive 
deficit angle is shown. The sphere on the right hand image indicates the points that are 
within the lightcone of the collision point. If we can choose the junctions within this 
area (as is the case in the image) then they are subluminal and all defects are timelike. 

However, if the deficit angle of defect A is large enough, then there will be no 
points on the defect B that lie within the lightcone of the collision. Consequently, 
it is impossible to choose the junction on the B defect such that it is timelike. This 
happens when the intersection of the deficit angle of A with the defect B lies outside 
the lightcone. In terms of the collision parameters this means 



+ ". > 1, (3.63) 



am 



where v is the velocity of the defects in the center of velocity frame (i.e. v = tanh?7/2). 
Consequently, the condition that all junctions can be chosen to be timelike is (in terms 
of 77) 



1 — cos^ (j) cos^ tAa/S ' 

tanh 77 < 



2 ain^ (/)COS^ 7/'b/2 



(3.64) 



1 — COS^ (f) COS^ 1 2 
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Before 



After 







1 

111 










A2 



Figure 3.15: Defect A is scattered off the stationary defect B, which has a surplus angle of 
TT. (Projected along the orientation of B.) The dashed line indicates the cut leading to the 
surplus area. On the bottom row the same process is depicted, but now as two separate pieces 
of space that are glued together along boundary of the white area. (The lines marked with a 
triangle are identified with each other, idem for the lines marked with a square.) After the 
collision there exist no straight lines going from A\ to A2. 
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Figure 3.16: Two defects colliding in the center of velocity frame. On the left, the geometry 
before the collision. On the right the continued geometry after the collision. The transparent 
sphere indicates the points that are within the light cone of the collision. 



Notice that for orthogonal {(j) = 7r/2) collisions this reduces to the limit where the 
junctions in the single defect continuation become lightlike, which we found in equation 

If the defect angles become even bigger (either because the defects have higher mass, 
or because the collision velocity is larger), than it can happen that the entire future 
lightcone of the collision point becomes contained in the wedges that are removed from 



the spacetime by the defects, as happens in the example shown in figure 3.17 



The limit where this happens can be obtained by finding the intersection of the edges 
of the two defects, and calculating the velocities of its edges. If all these velocities are 
greater than one, the wedges completely "eat up" the future lightcone of the collision. 
For an orthogonal collision of defects, this gives the limit 



,2 3 + cos tpA + COS -05 — COS IpA COS "0^ 

tann -q > 



(3.65) 



This is the same limit we found in equation | |3.11[ l for the velocity where the single defect 
in the single defect continuation became spacelike. 

It is no coincidence that we find the same limits as for the single defect continuation. 
For an orthogonal collision, choosing the junctions to lie at the intersection of the defect 
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Figure 3.17: If the deficit angles are large enough the sphere can envelope the entire future 
lightcone of the collision point. 



with the other defects deficit angle implies that the junctions for both ends of a defect 
coincide. In that case the intermediate defect connecting these two junctions disappears 
and the other four intermediate defects coincide to form a single defect. 

This result does however seem to be at odds with our earlier results for quadran- 



gle continuations as discussed in section 3.4 The algebraic solutions that we found 
there seemed to produce valid results in a much larger part of the orthogonal collision 
parameter space than the single defect continuations, both for finding continuations 
with only subluminal junctions (see figure 3.10| and for finding continuations with only 
timelike defects (see figure 3.11| . At the same time, we observed that the quadrangle 
continuations were special cases of the tetrahedral continuation. Consequently, one 
would expect to find valid tetrahedral continuations for at least the same part of the 



70 



3.6: Conclusions 



orthogonal collision parameter space as the quadrangle continuations. Something is 
not quite right. 

The problem we see here, is the result of the limitations of the algebraic description of 
a configuration of defects. Although a configuration is (almost) uniquely described by 
the Poincare holonomies of its non-contractible loops, an assignment of holonomies to 
the equivalence classes of loops does not always need to have a geometric realization. 
This happens to be the case for the supposed quadrangle continuation of orthogonal 
collisions in parts of the parameter space where there are no valid tetrahedral contin- 
uations. On closer inspection, these tetrahedral continuations correspond to configu- 
rations with junctions moving inward. Inward moving junctions are not, a priori, a 
problem, but in the case of an orthogonal collision they point to a geometric impos- 
sibility. For an orthogonal collision the two semi-infinite ends of a defect meet after 
the collision (if continued far enough). Inward moving jvinctions on these semi-infinite 
defects would imply that these defects meet before they end on their jimction — this is 
geometrically impossible. 

The geometrical argument used here to find the bound where the tetrahedral con- 
tinuation becomes incompatible with local causality can be applied to general, more 
complex continuations. If the deficit angles of the colliding defects completely "eat 
up" the future lightcone of the collision point, any continuation of that collision must 
start outside the lightcone of the collision. Since all the elements of a continuation 
originate at the collision point, this means that the continuation must contain at least 
some superliiminal elements. If some of the defects in the continuation have negative 
curvature, these can restore parts of the lightcone of the collision, but the continuation 
must still start outside the lightcone. 

3.6 Conclusions 

In this chapter we have studied the collisions of pairs of moving conical defects. Because 
the defects each carry an holonomy they cannot simply pass through each other. Instead 
an appropriate configuration of intermediate defects must be fovmd to continue the 
piecewise flat geometry after the collision. Each collision therefore creates new defects. 
Since in arbitrary configuration collisions are inevitable and the likelihood of the inverse 
process is negligible, the number of defects in a configuration steadily increases over 
time, ultimately approaching a fractal like structure. 

We have seen that for a large part of the parameter space of collisions, we can indeed 
find such continuations. When these continuations exist, they involve a number of 
free parameters. The continuations are therefore not unique. This indicates that our 
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piecewise flat model of gravity needs some additional physical input to select what 
continuation should occur for a collision. This could be taken to be an analogy of the 
fact that, in a continuous gravity theory, the evolution of a system is not only determined 
by the geometry of curved spacetime as dictated by Einstein's equation, but also by the 
matter lagrangian; the latter can still be chosen in many ways. 

More alarmingly, we find that for certain collisions it is impossible to find a continu- 
ation that contains only timelike junctions and defects. This occurs when the colliding 
defects satisfy a combination of the following conditions: 

1. At least one of the defects is very massive. 

2. The relative velocity of the defects is very high. 

3. The colliding defects come close to being coUinear. 

Since the local causality principle of our model forbids the appearance of spacelike de- 
fects (and to a lesser extent spacelike junctions), it appears that these collisions indicate 
that the model is possibly inconsistent. 

One might consider the possibility that one coiild restrict the model to configiira- 
tions where the problematic regions of the collision parameter space are not reached. 
However, although this may be possible for the first two conditions, it seems extremely 
imlikely that one could arrange a configiiration in such a way that no collisions would 
come near to being collinear. 
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At first sight, the piecewise flat model of gravity we have proposed seems very limit- 
ing. The model effectively eliminates all local structure from general relativity. Con- 
sequently, we expect — a priori — no long range gravitational fields or gravitational 
waves. This certainly seems to be at odds with the everyday world that we observe. 
Moreover, matter only seems to appear as extended objects in the form of straight 
cosmic strings. It is therefore a valid question, whether the model can serve as a model 
of real world physics at all. 

However, the conical defects in this piecewise flat model of gravity represent both 
matter and gravitational degrees of freedom (or rather a combination of the two). More- 
over, since the model allows both defects with positive and negative energy density (i.e. 
defects with a positive or negative deficit angle), one can easily imagine large scale 
configurations in which, on average, the energy-momentum (or equivalently the Ricci 
curvature) vanishes, while the Weyl curvature does not. Such configurations would 
represent effective large scale vacuum structure in the piecewise flat model. 

An indication that one can completely reproduce general relativity in this way comes 
from Regge calculus. f/S lOlj Regge calculus seeks to approximate general relativity by 



using arbitrary piecewise flat (more specifically simplicial) geometries. In that context, 
the continuum limit of piecewise flat geometries has been studied extensively. [ 29] [30| |4T} 
[45) It is known that the space of Regge manifolds is dense in the space of (pseudo)Rie- 
mannian manifolds (in some suitably defined way), and the Regge action may be ob- 
tained as the limit of the Einstein-Hilbert action. In particular, any (vacuum) solution of 
general relativity may be approximated to arbitrary precision by a suitably fine piece- 
wise flat geometryp^ Regge calculus, however, does not interpret the curvature defects 

^^Note, that although this is true on the level of geometries, the equations of motion derived from the Regge 
action may not converge due to persistent short wavelength oscillations.|18||19| 
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in the piecewise flat geometries as propagating fundamental degrees of freedom, and 
consequently does not restrict the type of defects that may occur. In the piecewise flat 
model studied in this thesis we, however, allow only non-spacelike curvature defects to 
occur. It is not immediately clear if and how this restriction affects the conclusion that 
the space of piecewise flat configurations is dense in the space of solutions of general 
relativity. This will be investigated in this chapter. 

Another issue is that the model only allows for a very particular source of energy- 
momentum in the form of straight cosmic strings. It would not be surprising if this in 
some way limits the types of matter sources that may be described by the model. This 
will be a second line of inquiry in this chapter 

The description of a configuration of defects simplifies dramatically in the limit 
where the energies of all defects are considered to be infinitesimal. This allows us to 
study continuous distributions of infinitesimal defects. Section |4.1| explains how to 
describe such a configuration in this limit using a density function on the configuration 
space of an individual defect and introduces the notation used in the rest of the chapter. 

Section 4.2 then constructs the energy-momentum tensor produced by an arbitrary 
configuration of physical defects. It finds the conditions that the energy-momentum 
will obey (and consequently must be obeyed by any theory that is to be approximated 
by this model in the limit of weak fields). 

The metric perturbation produced by a general configuration of physical defects is 
obtained in section 4.3 After this, in section 4.4 the results of the preceding sections are 
combined to find the metric perturbations that can be produced by a configuration with 
vanishing energy-momentum. We find that it is possible to reproduce the complete 
spectrum of gravitational waves found in linearized Einstein gravity. 



4.1 Defect configurations in weak field 

In chapter |2] two approaches to describing a general configuration of propagating de- 
fects were described. Both of these share that they rapidly increase in complexity as the 
number of defects increases, because the description of the state of a defect involves the 
states of other defects as well. This complication is this model's manifestation of the 
non-linear nature of general relativity. 

However, in this chapter we want to look at geometries generated by configurations 
with large numbers of defects, which makes the previously employed methods pro- 
hibitively complex. Fortimately, the description of a configuration of defects drastically 
simplifies in the limit that all defect angles are very small. In this limit each defect can 
be treated as a linear perturbation to a Minkowski background. Moreover, any new 
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intermediate defects that would be created by the collision of two defects are higher 
order in the defect angles of the colliding defects and can be neglected. 

Consequently, in the linear limit the state of a defect line can be described while 
ignoring the presence of other defects. A configuration of defects can therefore be 
completely described by giving the number of defects in any given state. Since the effect 
of two defects with the same state is simply that of a single defect with the combined 
energy of the two defects, we can completely describe the configuration by giving the 
energy in each possible state. That is, if Ai is the state space of a single defect with unit 
energy density, a configuration of defects can be described by a distribution p : A4 ^ Wi 
giving the energy density in each state. 



To parametrize the state space Ai, recall that in 2.3 we showed that the state of a 
single defect line with energy density p could be described by the following data: a 
vector X that gives the position of the defect, a vector d that gives its direction, and a 
vector V that gives its velocity. 

The triple (X, d, t7) is enough to uniquely identify the state of the defect. However, 
this characterization is not unique, since the triple 

{X + ad,pd,v + jd) (4.1) 

describes exactly the same state. In fact, the triples describing the same state as (X, d, v) 
can be completely parameterized by the numbers a, (3, and 7. The state space of a 
line defect with imit density M is thus given as the quotient of R"^ x R'^ x R^ by the 



equivalence relation defined in 1 4.1 1 



A set of unique representatives for each equivalence class in M can be formed by 
taking a triple {X, d, v) that satisfies the following conditions, 

V J- d, 

X _L d, and (4.2) 

Such a representative triple is unique up to a sign of d, i.e. to get a unique representative 
d should be viewed as an element of the real projective plane, RP'^. 

In this chapter it will be convenient to decompose density functions p on yVt in a 
set of canonical functions that we call laminar plane waves. A laminar plane wave is 
a configuration where all defects have the same direction dp ^rid velocity wg and the 
density p is a plane wave function with wave vector fco with respect to the position X of 
the defects. That is, the density function [fco, ^o, ^^o] corresponding to a laminar plane 
wave with wave vector fco, direction do/ and velocity uo can be written as a function of 
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(X, d, w) representing an equivalence class in Ai as follows, 

p'''"[ko,do,vo]{XJ,v) = e^^''^°-^S{d-do)Siv-vo)Q (4.3) 

Since X en w are always perpendicular to d it is sufficient to consider only laminar plane 
waves with kg and Uq perpendicular to do- The set Ai of laminar plane waves with 
parameters (fco, do, uo) is complete in the sense that any density function p{X, d, v) on 
Ai can be written as 

p{X,d,v)= I dkoddodvop^'''^[kQjn,VQ]{X,d,v)p{ko,do,vo) (4.4) 

JM 

dke'^'''''-^p{kj,v), (4.5) 

k±d 

where p{X,d,v) is a density fxmction on J\/l. The function p{k,d,v) can therefore be 
viewed as a partial Fourier transform of p{X, d, v)^^ 

It will also be convenient to split the parameter w in a component collinear with k 
and a component perpendicular to both k and d, 

^ k k X d 

v = uj^ — \-v^^ — , (4.6) 
k^ k"^ 

where the ambiguous direction of d is chosen such that v is non-negative. The usefulness 
of this split becomes apparent when we do a Lorentz transform. The Lorentz transform 
of a laminar plane wave with parameters (fc, d, uj, v) is again a laminar plane wave, but 
with different parameters. It can be shown that the combination k^ = {oj, k) transforms 
as a 4-vector under the Lorentz transformation. In section l4!2l we will see that this is the 
wave vector of the corresponding energy-momentum tensor. In the remainder of this 
chapter we shall denote the parameters of a laminar plane wave as (fc^, d, u) . 

The physicality condition that we impose on the defects implies that the velocity of 
each defect must be smaller than or equal to c = 1. By squaring equation | |4.6^ we see 
that this implies that 



w2 



fc2 



< 1. (4.7) 



Consequently, we see that for a laminar plane wave of physical defects, the wave vector 
k^ must be spacelike or lightlike and v must be smaller than or equal to ^Jk^Ji^^. 



^^The density function of course needs to be real. The expansion in complex exponentials is for tlie sake of 
convenience only. To obtain physical results we should consider only the real part. 
^*Conversely, p{fc, d, v) may be obtained from p{X, d, v) by an inverse partial Fourier transform. 
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4.2 Energy-momentum 

In this section we will derive the energy-momentum tensor generated by an arbitrary 
configuration of physical defects in the limit that the distribution is continuous and all 
defect angles are small. This will tell us what conditions are imposed on the energy- 
momentum tensor by the requirement that the defects are physical (i.e. non-tachyonic). 
This puts limits on the kinds of models which can be found as the continuum limit of 
the piecewise linear model considered here. 

In the linear weak field limit the energy-momentum tensor of a configuration of de- 
fects can be found by adding together the energy momentum tensors generated by the 
individual defects. Therefore, if 3),^ [X, d, v] (x^) is the energy-momentum generated by 
a single defect with position X, direction d and velocity v with unit energy density, then 
the total energy-momentum of a configuration of defects given by a density function 
p{X, d, v) on M is 

T[p]^,{x,)= fdXdddvp{Xj,v)f^,[xJ,v]{x,). (4.8) 
Jm 

Consequently, if we have an explicit expression for T^^ [X, d, v] (x^), we can compute 
the energy-momentum tensor for any configuration of defects. Alternatively, since 
the laminar plane waves form a complete basis for all configurations, the total energy- 
momentum of a configuration can also be obtained from T^^^ [fc^, d, u] (x^), the energy- 
momentum of a laminar plane wave with wave vector A:^, direction d and perpendicular 
velocity v through, 

T[p]^u{x^) = / dkdddv pik^,d,v)f^;:^[k^J,v]{x^). (4.9) 

JM 

We will now derive the energy-momentum tensor T^™ [fc^, d, v] (xk) generated by a 



single laminar plane wave. In section 2.1 we derived the energy-momentum tensor of 
a single stationary defect through the origin and directed along the z-axis. With unit 
energy density the result is. 



f^,[0,z,6\{x^) = 



/I \ 





Vo -i/ 



S{x)Siy). (4.10) 



The dependence on the position X = (X^;, Xj^, Orj can be obtained by performing 



^^Remember that X should be perpendicular to ( 
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appropriate shifts, which yields 








\ 




















Vo 





-1/ 



6ix-X.,)5iy-Xy). 



(4.11) 



By combining this result with equations 1 4.3 1 and | |4.8) we obtain the energy-momentum 
generated by a stationary laminar plane wave with wave vector — (0, k), direction z 
and zero velocity, 

fjf^[kx,z,0]{x,) = [ dXdddvp'''^[kx,z,0]{Xj,v)f^,[X,d,v]{x^) 



2iTik-X 



(4.12) 



M 

dX%, [X,z,Q\{x^)e 

/I \ 



Vo -i/ 

Any other laminar plane wave can be obtained by applying appropriate Lorentz 
transformation A^. To write the result we first notice that the right hand side of equation 
( 4.12) can be written as 

(m^m, + u^d^d,) e2"''=-% (4.13) 
if we introduce the 4-vectors = (1, 0, 0, 0) and = (0, 0, 0, 1). 
After applying and a linear redefinitiorPj 



d!^ = aAi;d^ 



such that 



d'^d'^ = 1, 



(5A>^ 



0, 



the direction d and velocity v of the new laminar plane wave can be found as 

u'^ — ( 1 , , and 

d'={Q,d). 



(4.14) 



(4.15) 



(4.16) 



^^The 4-vectors and span the plane of a single defect in the laminar plane wave. Applying a Lorentz 
transformation yields vectors spanning the plane of a defect in the transformed plane wave. However, these 
vectors will not correspond to the (non-covariant) parameters we introduced to describe the velocity and 
direction of the defect. To find a pair of 4-vectors u'^ and d'^ that span the same plane but correspond to the 
velocity v and direction d of the defect we need to do a linear transformation. 
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Applying the same Lorentz transformation and linear redefinition to equation < 4.13| 
yields the energy-momentum of the new laminar plane wave 

f';^[k^, d, v]{x,) = --^ {u'^ul + u'^d'^d'J) e^-^^-. (4.17) 

This gives us the energy-momentum for all laminar plane waves with fc^fc^ > 0. Notice 
that a laminar plane wave of defects with wave vector fc^ only contributes to the Fourier 
mode of the energy-momentum tensor with wave vector fc^,. As a result the Fourier 
transform of the total energy-momentum of a configuration specified by the density 
function p on has the especially simple form, 

T^Mikx)^ / dd[dv^kxJ,v)fPl[kx,d,v]. (4.18) 



Jd±k Jo 

The energy-momentum tensor for laminar plane wave with fc^fc'^ — 0, can be found as 



a limit of equation (4.171. There are two possibilities. 



1. fc^ — >■ corresponding to a constant "wave" of defects with the same direction and 
velocity. 

2. UpU^ — > corresponding to a laminar plane wave of lightlike defects. 



The first case is easy enough to compute since equation ( 4.17| is regular in this limit. 



However, in the second limit the 1 / u'^ factor blows up. This can be fixed by noting that 
the normalization of is arbitrary, and we are therefore free to rescale it by a factor 



-tt'^ changing equation | 4.17| to 



flTikxJ, v]ix,) = [u'^u', + u'^d'^d',) e^-fe-, (4.19) 

which is regular in the limit that u'"^ goes to zero. This fixes Tqq'" to be 1, which has 
the additional advantage of giving p the physical interpretation of the energy density 
present in a particular mode of laminar plane wave. 

We now have all the ingredients we need to calculate the energy-momentum tensor 
of a general configuration of physical defects and can explore what conditions this will 
satisfy. In particular, we are interested in what conditions are imposed on the energy- 
momentum tensor by the restriction that all defects must be physical. 

One immediate condition that we observe from equation | |4.17^ is that fc^T^™ — 
because fc^d^ — k^^u^ = 0. Consequently, equation 1 4.18) implies that the total energy- 
momentum must satisfy 

fc^T^,[p] = 0, (4.20) 
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for any distribution p on M. That is, energy is conserved, as one expects from any 
reasonable physical theory. 

Other conditions can be obtained by examining T^^^ [jj] {k\ ) mode by mode. We have 
already observed that laminar plane waves with wave vector /c^ only contribute to 
modes of the energy-momentum tensor with the same wave vector. Therefore, since 
laminar plane waves of physical defects must have fc^fc'^ > 0, it is impossible for a 
configuration of defects to generate an energy-momentum tensor with Fourier modes 
with fc^fc^ < 0. 

For modes with fc^fc^ > 0, by Lorentz invariance, we can restrict ourselves to the 
special case that fc^ = (0, k, 0, 0) without loss of generality. The direction d, which 
must be perpendicular to k, can then be parameterized by a single angle (j> with = 
corresponding to the direction of the z-axis. The Fourier transform of equation 1 4.19) 
then becomes. 



( 1 



I COS( 

Vf sin^ 



COS( 



sm 



fc2 

COS 6 sin ( 



sin( 



cos ffl sm ( 



(4.21) 



fe2 



COS^ / 



If we expand p{k^, d, v) as 



,sinm0)), (4.22) 



where the Pn are Legendre pol5momials and the coefficients rm„ are understood to be 
functions of fc^, then the total contribution of the fc^-mode of p to the energy-momentum 
(OH) becomes. 



Tf,M{kf,) 



/ 2roo 



rii+rio 

\?'ii+7'io 



rii + rio 


^''00 + ?'01 + |»'02+r20 
^20 



rii + rio 


f20 



. (4.23) 



'^roQ+roi + ^ ro2 - r2o/ 



If we compare this to the most general form of a mode energy-momentum tensor with 
wave vector fc^ — (0, fc, 0, 0) that satisfies fc^T^,y — 0, 



Tfj,^{kf^) 



/Too Toi Tq2 



Tqi T22 T23 

\To2 T23 733/ 



(4.24) 



80 



4.2: Energy-momentum 



then we find that any such energy-momentum tensor can be generated by appropriate 
choices of the coefficients r„j„. We therefore find the physicality condition on the de- 
fects puts no further restrictions on the modes of the energy-momentum tensor with 



k^kf > 0. 



In the special case that fc^, is lightlike, we can assume, without loss of generality, that 



(w, uj, 0, 0). In this limit 14.19 1 becomes. 





1 





'\ 


1 


1 




















{o 








0/ 



(4.25) 



and the total contribution of a distribution of defects p to the Fourier mode of the 
energy-momentum with /c^ = {ui, uj, 0, 0) becomes 



/rao roa 

''oo roQ 



V 








(4.26) 



Apparently, the lightlike modes of the energy-momentum tensor generated by a con- 
figuration of physical defects are subject to additional restrictions. Not only do these 
modes have to be transverse, but they also cannot have any pressure or momentum 
perpendicular to their direction of propagation. This restriction can be formalized in 
the following way: For any lightlike mode of the energy-momentum tensor T^^{kf^) 
and any lightlike vector the contraction l'^T^^{kfi) is a non-spacelike vector 

The other special case, k^ = is somewhat different since any direction d will be 
perpendicular to k ~ (0, 0, 0). The total contribution to the zero mode of the energy- 
momentum tensor will thus be found by integrating equation 1 4.19) over all mutually 
perpendicular v and d. That is. 



dvAd p(0, d, v) {u^Ui, 



u^dndp 



(4.27) 



Since tensors of the form w^Wi, + v?d^dy span the space of symmetric 2-tensors, we can 
conclude that any zero mode of the energy-momentum tensor may be produced by a 
configuration of defects. The condition that it contains only physical defects poses no 
further restrictions. 

We therefore obtain the following restrictions that the energy-momentum tensor of 
a configuration of physical defects will satisfy in the continuum weak field limit 



81 



Chapter 4: Continuum limit 



1. fc^T^,(fcA) = Oforallfc^. 

2. T^u{k\) = for all fc^ with fc'^fc^ < 0. 

3. l^T^^{kx) is a non-spacelike vector for all lightlike fc^ and Z^. 

The first and third condition hold for most physically reasonable theories. The first 
expresses conservation of energy-momentum, while the last is normally imposed as 
part of the null dominant energy condition which is employed in cosmology to ensure 
vacuum stability while allowing negative vacuum energy. |27j 

The second condition is satisfied by various simple matter models used in general 
relativity, such as dusts. However, it is typically violated in classical wave like sys- 
tems. For example, consider a standing wave solution of the Klein-Gordon equation, 
(j) = cos(aji) cos(fc • x). Even if fc^ > w^, the energy-momentum tensor — which be- 
haves like the square of — will have terms which are proportional to cos(2a;i) and 
consequently will violate the second condition. 

This indicates that the model cannot represent all types of matter at linear order. 
At this level all interactions are neglected, and we end up with a system that is very 
similar to a dust of non-interacting point particles. Beyond linear order defect lines will 
collide in a non-trivial manner, as was discussed in chapter |3] The energy-momentum 
tensor corresponding to the continuation of a collision will t5rpically violate the second 
condition. 

At non-lrnear level we need to choose a prescription, which tells us how each colli- 
sion is continued, to complete the dynamics of the model. Different prescriptions will 
lead to different types of interacting matter. At this point it is unclear what types of 
limitation on the energy-momentum tensor — if any — will persist once all possible 
types of interaction are included. 

4.3 Metric perturbations 

We now turn to the effect that a configuration of defects p(fc^, d, v) has on the metric. 
In the limit of weak fields the metric g^^^ can be separated in a static Minkowski back- 
ground r/^,y and a small perturbation 

gp^u{x\) = + hf^^{xx). (4.28) 

If we consider only the linear perturbations caused by the presence of a defect, then the 
total perturbation caused by a continuous distribution of defects can be found as the in- 
tegral of the perturbations of individual components. Consequently, if /iJj™[/ca, d, v]{Xf^) 
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is the perturbation of the metric caused by a laminar plane wave with wave vector 
kf^, direction d, and perpendicular velocity v, then the total perturbation caused by a 
configuration p is given by. 



•-i^ApK^k) ^ I dkdddv p{kx,d,v)h^^[kx,d,v]{xi^). 
M 



(4.29) 



Therefore, if we know h}^^[kx,d,v] for any combination of the parameters {kx,d,v), 
than hf^i, [p] can be computed for any configuration p. 

To calculate /ij^™ [kx, d, v] {x^) we first need to fix its linear gauge freedom, which we 
do by setting the gauge condition d'^h^y — 0|^ With this choice the linearized Einstein 
equation becomes (in its Fourier transformed form). 



(4.30) 



where h is the trace of h 

When k'^ki^ ^ 0, equation 1 4.30) can be inverted to obtain the linear metric pertur- 
bation as a function of the energy-momentum tensor In particular, if Tj^^ [kx , d, v\ is 
the Fourier mode of the energy-momentum tensor generated by a laminar plane wave, 
then the metric perturbation generated by that laminar plane wave is given by a single 
Fourier mode. 



?,lam 
''fj.i/ 



[kxJ,v] = ^(-^M^^ - - '^)v''^)f':^[kxj,v] 



(4.31) 



We can therefore study the effects of a distribution of defects cj), v) on a per mode 
basis. By plugging in the fjf^[kx,d, v] from equation < 4.19) , we find that. 



27r2fc- 



kj_iku V. 



(4.32) 



where — (1, u) and d^ — (0, d). 

When /c^fcp > 0, we can assume by Lorentz invariance that k^ = (0, fc, 0,0). Pa- 
rameterizing d as (0, — sin (/), cos (/)) and v as (0, w cos t; sin(/)), we obtain the explicit 
expression. 



I fe 



Up[kx,d,v] = 



1 



27r2fc2 



fcos0 



COS( 



COS 



COS d) sin < 



sine 



cos <p sm ( 

„:„2 / 



(4.33) 



^''There is some residual gauge freedom for tlie liglitlike modes of ttie metric perturbation as we will discuss 
later on. 
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The total metric perturbation caused by a configuration of defects given by a distri- 
bution i5 on M. can be obtained from h}i^^[k\, d, v] through the integral. 



hf,Ap\ik\) = / dddv p{kx,d,v)h^^^[kx,d,v 



(4.34) 



Performing this integral in the case that fc^ — (0, k, 0, 0) and applying the expansion of 
p as given in | |4.22) yields. 



hf_,„[p]{ki^) = 



V 





rio + rii 
rio + rii 







(4.35) 



When fc^fc^ = the linearized Einstein equation ( 4.30) cannot be inverted and equa- 
tion 1 4.32) cannot be applied directly. However, we may obtain the metric perturbation 
for these modes as a limiting case of the modes with fc^fc^ > 0. 

In the case that fc^ becomes lightlike we can assume due to Lorentz invariance that 
it goes to ~ (w,aj,0, 0). Such a laminar plane wave is the limit of waves with 
momentum fc^ — (w, k, 0, 0) as k lj. Since physicality requires that < w < y/ — uP- , 
V must simultaneously go to zero. Simply setting v = and using fc^ — (w, k, 0, 0), 
Up — (Ijw/k, 0, 0), and c?^ ~ (0, 0, — sm(j), cos(/)) in equation ( 4.32} yields. 
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(4.36) 



Since only the pre-factor depends on k the limit as k goes to oj is straightforward. 

Letting v go to zero by another route will lead to a result that (in the k lo limit) 
differs from the above by a term of the following form 



a 6 6 

6 6 
V6 6 



6\ 



0/ 



(4.37) 



where the are arbitrary (possibly infinite) parameters. 

Such a contribution can be gauged away. Under an infinitesimal coordinate trans- 
formation given by a vector field the /i^j, transforms as. 



(4.38) 
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The gauge condition = implies that 



(4.39) 



This completely fixes £.f^(kx) for fc^fc,, ^ 0. However, when fc^fc^ = it only implies 
that k'^£^^{k\) = 0. The residual gauge transformation subject to that condition for 
= {ijj,uj, 0,0) takes the form of equation ( |4.37| . Contributions of that form can 
therefore be gauged away, and we are free to adopt ( |4.36| as the gauge fixed form of 

h'^^[{uj,Lo,0,0),^,0]. 



Inserting | |4.36[ l in the integral | |4.34[ | and using the expansion 14.22 1 for p yields the 
lightlike modes of the metric perturbation caused by a configuration of defects given by 
a distribution p 



h^,y[p]iuj,uj,0,0) 
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Vo 






roo + r2o 
ha 






roa - r2oJ 



(4.40) 



In the case that fc^ — > 0, the equation 1 4.32} diverges. This signals a breakdown of 
the linear perturbation approach in the limit of constant fields. This is not unexpected, 
since constant non-zero energy densities will typically lead to non-trivial effects on the 
global level with respect to the topology or causal structure. 



4.4 Gravitational waves 

In the previous two sections we have obtained the energy-momentum and metric per- 
turbation caused by a continuous distribution of defects at the linear level. At the 
microscopic level of individual defects the energy-momentum tensor completely fixed 
the metric (up to gauge transformations) due to the ad hoc rule we imposed that the 
vacuum should be completely flat. We are now able to answer the question whether 
this property persists in the continuum limit. 

Since we are considering the linear limit, it is enough to consider what additional 
metric structure may be present for a configuration with zero energy-momentum. In 
we obtained a complete expression for the energy-momentum of a config- 



section 



4.2 



uration of defects described by a distribution p on A^. Requiring that this expression 
(equations 1 4.23) and | 4.26| ) vanishes implies the following conditions on the coeffi- 
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dents rnmik^) of p in the expansion | 4.22[ ; 

f2a{kf,)^0, rii(/c^) +rio(fc^) = 0, > for /c^fc^ > 0, and 
^■20(^^0=0, + fio(fc^) = 0, J 

roo{k^) = 0, } forfc^fc^-O. 



(4.41) 
(4.42) 



The metric perturbation caused by a vacuum configuration of defects can now be 
found by applying these conditions to the complete expressions 1 4.35) and 1 4.40) , found 
in section 4.3 for the metric perturbation caused by a configuration of defects. This 
yields that for k^kfj, > 0, h^y{k\) vanishes when 1),^ vanishes, as one would expect 



since the linearized Einstein equation is rnvertible for k^k^ ^ 0. 

However, for fc^fc^ — 0, the requirement that T^^y{k\) vanishes only fixes roo{k^) to 
be zero. The coefficients r2o{k^) and r2o(fc^) are unconstrained. Consequently, for each 
lightlike k^ there exists a two parameter family of vacuum metric structures. When 

= (w, cj, 0, 0) these are described by 



/i^^[p](w,a;,0,0) = 



1 



27r2w2 



/O \ 



r2o 7=20 

\0 f2o -r2o/ 



(4.43) 



and 



where we immediately recognize r2o and 7=20 as the coefficients of the familiar 
polarizations of gravitational waves. 

This answers the question we posed in the introduction of whether any of the vac- 
uum structure of general relativity would be recovered in the continuum limit of our 
model. The answer turns out to be affirmative. In fact, we find that all vacuum (i.e. 
Ricci flat) metrics that are available in linearized general relativity may be found as the 
continuum limit of a sequence of configurations of physical defects. 



4.5 Conclusions 

We have studied continuous distributions of conical defects. In the limit of weak fields, 
considering only linear contributions, the energy-momentum tensor of such a distribu- 
tion turns out to satisfy certain conditions. These conditions can be expressed as follows 
for the Fourier transform of the energy-momentum tensor, 

1. fc^T^^(fcA) = Oforallfc^. 
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2. Ti^t^ikx) = for all fc^ with fc'^A;^ < 0. 

3. l^Tni,{k\) is a non-spacelike vector for all lightlike and Z^. 

The first and third conditions are satisfied by most reasonable classical theories. The 
second condition, shows an inability to model wave-like phenomenon at the linear or- 
der. This property is shared by other non-interacting matter models, like dusts. Beyond 
linear order line defects collide in a non-trivial manner causing violations of this second 
condition. It would be very interesting to see if any restrictions remain once all possible 
interactions are included. 

The metric perturbation caused by a configuration of defects turns out to be mostly 
fixed by its energy-momentum. That is if two continuous distributions of defects have 
the same energy-momentum tensor, they also produce the same linear metric perturba- 
tion, at least for most modes. 

The lightlike modes of the metric perturbation form an exception, as they are only 
partially fixed by the energy-momentum tensor. The metric perturbations of two con- 
figurations with the same energy-momentum may differ by a transverse traceless light- 
like mode, i.e. by a gravitational plane wave. We thus see that the piecewise flat model 
of propagating defects — even though it does not contain any a priori vacuum structure 
— recovers the vacuum structure of general relativity in the continuum limit, at least at 
the linearized level. 
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Gravitational waves 



In the previous chapter we showed that our piecewise flat model of gravity contains 
gravitational waves as an emergent feature in the linearized continuum limit. In this 
chapter we extend this result to exact gravitational waves. In contrast to linearized gen- 
eral relativity, it is not possible to derive completely general exact vacuum solutions of 
general relativity. Deriving exact solutions requires some form of additional symmetry. 
For example, it is known how to construct completely general gravitational waves with 
cylindrical symmetry, ]100[ and how to construct exact gravitational plane waves. p6) 

Here we will focus on plane gravitational waves and show how to construct a family 
of increasingly fine discrete configurations of defects that approach an exact gravita- 
tional wave in the limit that the configuration becomes continuous. 

Our strategy will be to first consider gravitational plane waves in the linearized 
limit, and examine what linearized distributions of defects approximate these waves 



(section 5.1 1. Based on the linearized configuration, we will make an ansatz for an 
exact piecewise flat planar wavefront, which reduces to the linearized distribution in the 
linear limit. We show that this ansatz is indeed a bona-fide piecewise flat geometry, and 
that in the limit where the configuration becomes continuous it produces an impulsive 
gravitational wave solution of general relativity (section [5.2| . In section 5.3 we combine 



trains of these wavefronts to form general exact gravitational plane wave solutions. 



Finally, in section 5.4 we comment on the construction of non-planar waves and the 



additional subtleties that those entail. 
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5.1 Linear plane waves 

For simplicity we shall consider, without loss of generality, plane waves travelling in 
the positive z-direction. Expressed in lightcone coordinates (u, v, x, such a wave 
takes the following form in linearized general relativity. 



/O \ 



0/1+ hy, 

\0 hy, -h+J 



(5.1) 



where hf^^ is the perturbation of the metric on the Minkowski background. In chapter |4] 
we saw how to produce such a metric perturbation in linearized piecewise flat gravity. 
The basic building block was a configuration that we called a "laminar plane wave of 
defects". A laminar plane wave of defects is a space-filling configuration of parallel 
conical curvature defects (i.e. line defects that have the same orientation and velocity), 
whose energy density varies as the amplitude of a plane wave with a wave vector 
perpendicular to the defects. That is, given a point in our spacetime there is exactly 
one defect (spanned by 4-vectors and u^) in the laminar plane wave configuration 
that passes through this point, because all the defects in the laminar plane wave are 
parallel. The linear energy density of this defect is given by fls"'^'^^'' with some fixed 
wave vector fc^, which is perpendicular to and u^. 

In the linearized limit of our piecewise flat model an arbitrary configuration of de- 
fects may be written as a linear combination of laminar plane waves. In chapter |4] we 
foimd that a gravitational wave of the form | |5.1[ l was produced as a superposition of 
laminar plane waves of defects with the velocity — (1,0,0,1) and wave vector 

= {uj,0,0,Ld). The orientation c?^ of the defects in these laminar plane waves can 
then be parametrized by a single angle, cj), which we choose to be the angle between the 
defect lines and the x-axis. 

To be a gravitational wave the total energy density amplitude of the superposition of 
laminar plane waves must vanish. That is, if the energy density amplitude of the lam- 
inar plane wave with orientation (f) is called ^{4>)/ then this condition can be expressed 
as 

/•TT 

n{(P)d<j> = 0. (5.2) 



It was established in section 4.4 that this is the only condition on the superposition of 



laminar plane waves of defects required for it to describe a gravitational plane wave. 



^^Recall our conventions that spacetime metrics have signature ( — h ++), lightcone coordinates are defined 
as M = 2±i and v = and we use units such that c = h = SttG = 1. 

V2 V2 



90 



5.1: Linear plane waves 




Figure 5.1: A laminar plane wave consists of parallel line defects travelling in the same 
direction, whose densities vary as a plane wave in the direction of propagation. 



In equation 1 4.43[ , the polarization of the resulting gravitational wave was shown to be 
given by 



cos(20)d0, 



(5.3) 



hx = r2a = / ^i4>) sin(20)d0. 
Jo 



All higher moments of are unconstrained. That is, there is a large family of linear 
configurations of defects that produce the same linear gravitational wave. We therefore 
have quite a bit of playing room for finding a linearized configuration that can be 
interpreted as the limit of some exact configuration of defects. 

Geometrically, it would be convenient if fl{cj)) has as few discrete components as 
possible. Condition | |5.2^ implies that we need at least two components. This is, in fact, 
sufficient to create waves of any polarization. If we take 



(5.4) 



then the polarizations are given as /i+ = 2a cos 2(/)o and — 2a sin 20o- We thus have 
perpendicular components of opposite energy, with the amplitude of the defects deter- 
mining the amplitude of the wave, and the orientation determining the polarization. 
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Figure 5.2: A wavefront of perpendicular positive (white) and negative (black) energy line 
defects. 



5.2 Plane wavefronts 



This linear configuration prompts us to make an ansatz for an exact configuration of de- 
fects. Further discretizing the linear distribution, we guess that a configuration consist- 
ing of subsequent grids of perpendicular positive and negative energy massless defects 
should suffice to approximate a gravitational wave. To analyse this configuration we 
first consider a wavefront located at u = ^ = consisting of a single grid (see figure 

We first need to check if this ansatz defines a proper configuration of defects. In 
| |92| the conditions that should be met at a junction of multiple defects were discussed. 
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These can be phrased in terms of the holonomies of loops around the defect lines. They 
can be summarized as saying that the product of the holonomies of a set of loops whose 
concatenation is contractible should be the identity. For the proposed grid configuration 
the junction conditions must be satisfied at each vertex of the grid. 



Since all the horizontal line defects in figure 5.2 have the same density, the holonomy 
of each simple loop around a section of a horizontal line defect should be the same, and 
similarly for the holonomies of simple loops around vertical line defects. If we call 
the holonomy of a simple loop around a section of a horizontal grid lines Qa and the 
holonomies of the vertical grid lines Qb, then at each vertex the junction condition is 



QaQbOa^'Qr' =1. (5.5) 

The holonomy of a lightlike defect is a null rotation, i.e. a Lorentz transformation 
that leaves a lightlike and a perpendicular spacelike direction invariant. In the PSL(2, (D) 
representation of the Lorentz group, the holonomy of a general lightlike defect moving 
in the ^-direction is represented as 

where C is a complex number. The argument of C gives the direction of the defect, while 
its modulus gives the energy density of the defect, i.e. the parabolic angle of the null 
rotation. The junction condition for n lightlike defects moving in the z direction with 
complex parameters Ci , • ■ • , Cn therefore reduces to 

Ci + • • • + Cn = 0. (5.7) 

Furthermore if Qa has complex parameter Qa then Qa^^ has complex parameter —Qa- 
Hence the junction condition for the vertices of the grid | |5.5) reduces to 

Ca + Cfc - Ca - a = 0, (5.8) 



which is true for any Qa and Qb- Our ansatz configuration of a grid of perpendicular 
lightlike defects is therefore a bona-fide configuration of defects for any value of the 
amplitudes of its components. In particular, it is valid for the particular case we are 
interested in, where the components have opposite amplitudes. 

Our next step is to compare this configuration with an exact plane wave solution 
of the vacuum Einstein equations. For this we need to find a metric to describe the 
geometry of the grid of defects. 
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I 




Figure 5.3: The geometry of one "tube" of spacetime.The deficit angles of the positive energy 
defects ( white), cause the geometry to be squeezed in the y-direction as u increases. Similarly 
the surplus angles of the negative energy defects (black) cause the geometry to be elongated 
in the x-direction. 



Away from the grid the spacetime is flat. Hence, if we position the grid at u — 0, we 
can choose the metric for u < to have the familiar Minkowski form, 

ds^ = 2diidt; + dx^ + dy^ (5.9) 

For simplicity we choose to align the grid with the x and j/-axes|^such that the positive 
energy defects are aligned in the x-direction and sit at constant values of y = e{n + ^) 
for integer values of n and grid spacing e, and the negative energy defects are oriented 
in the y-direction and sit at constant values of a; = e{m + i). 

The defect and surplus angles of the grid lines can be oriented along the direction 
of propagation of the grid, i.e. the positive u direction. The cuts of the defect /surplus 
angles then divide the w > side of the grid in square "tubes" of spacetime (see figure 



5.3 1. The Minkowski metric from m < side of the grid can be continued into each tube. 
The constant u and v slices in each tube start out as a square at the base {u — 0) of the 
tube. As u increases the square is elongated in the x direction and compressed in the y 
direction. 



^'This should correspond to wavefront with "+" polarization, according to the linear analysis. 



94 



5.2: Plane wavefronts 



We now construct coordinates on each tube, which can be patched together to one 
coordinate system on the u > side of the spacetime. The space 

T™,„ = R>o X R X [e(m - i), e(m + \)] x [e(n - i), e{n + \)] (5.10) 

can be mapped to the tube behind the face of the grid centered ot\ x — em and y — en 
by the mapping, 

V 1-^ V, 

a; I— )■ em + (1 H — em), (^-H) 

y I— ^ ert + (1 u) (y — enj . 

The metric induced on T,„ „ by this mapping is, 

ds^ = 2diidt; + ((1 + f u)da; + f (x - em)du)^ + ((1 - ^u)dy - f(y - en)du)^ . (5.12) 



Together, the patches „ cover the entire u > side of the grid. The metric | |5.12[ | may 
be viewed as a metric on the imion of all „, R>o x R^, if m and n are interpreted as 
functions of x and y respectively. That isj^ 

mix) = [x/e + 1/2J , ^ 
n(y) = LyA + l/2j. 

These functions are not continuous, so we need to worry whether ( 5.12| describes a 



well-defined geometry on the whole u > half-space. The appropriate condition |31 



|56] is that the metric induced on the common boundary of two patches is the same in 
both patches. 

On the x = e(m + i) boundary the metric on Tm,n induces the 3-dimensional metric, 

ds^ = 2dudv + ((1 - ^u)dy - f (y - en)duY ■ (5.14) 

On the other side of the boundary the metric on the r„i+i,ri/ can easily be seen to 
induce the same 3-dimensional metric. Similarly, we find that the 3-dimensional metrics 
induced on the y — e{n + i) boundaries match, and hence that the metric 1 5.12) is well- 
defined on the whole u > half-space. 

Moreover, the 3-dimensional metric induced on the u — hyperplane is 

ds^=dx^+dy^, (5.15) 



"Here [x-J denotes the floor function, which returns the largest integer smaller than or equal to x. 
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which agrees with the metric induced from the Minkowski side. Consequently, the 
metric | |5.9) for u < and the metric | |5.12| for u > together describe a well-defined 
geometry on the whole spacetime, given by the metric 



ds^ = 2dudv + ((1 + fue{u))dx +f{x- em)e{u)du) 



2 



+ ((1 - f w6l(u))dy - f (y - €n)9{u)duf , (5.16) 

where 9{u) is the Heaviside stepf unction. 

This metric should describe the geometry of our piecewise flat configuration of de- 
fects. As a check, we can calculate the Riemann tensor of this metric. Since the metric 
is discontinuous along the boundaries of the patches, we need to treat the curvature 
tensors as tensor distributions as described in [56|. The upshot of this is that in a case 
like this we can and should treat the metric components as distributions rather than 
normal functions. Using the fact that in a distributional sense 

^[x\=J2H^~'i) = Mx), and (5.17) 

—ue{u)^S{u), (5.18) 

we find that the Riemann tensor is given b}]^ 

R = 4aS{u){A,{y + |)(dM A d?;)^ - A,{x + f)(du A da;)^), (5.19) 
where (x) is the Dirac comb with period e, i.e. 

A,{x) = ^S{x~ei). (5.20) 

We observe that the Riemann tensor indeed vanishes everywhere, except on the grid of 
defects where it has singular delta peaks as expected. We have therefore found a metric 
for our piecewise flat spacetime consisting of a grid of lightlike defects. 

Our next step is to find the continuum limit of this spacetime. For this we will send 
the grid size e to zero while keeping the density of defects (approximately) constant. The 
number of defects in a unit area of the grid scales inversely proportionately to e, hence 
the amplitude of the defects a should scale linearly with e to keep the total amplitude 
per unit area constant. Inserting a — eao in the above metric and taking the limit of e 
goes to zero, we find 

ds^ = 2dudv + (1 + aoue{u)f dx^ + (1 - aQue{u)f dy^ , (5.21) 



^^For notational convenience we denote the tensor product of a 2-form with itself by a square, i.e. 
{du A dy)^ = {du A dy) ® {du A dy). 
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where 6{u) is the Heaviside step function. The Riemann tensor (distribution) is given 
by 

R = Aao5{u) ((du A dyf - [Au A Axf) , (5.22) 

from which the Einstein tensor may be obtained, which is seen to vanish. Consequently, 
the metric 1 5.21) is an exact solution of the vacuum Einstein equation, a gravitational 
wave consisting of a single plane wavefront. Waves of this type are called impulsive 
gravitational plane waves, and were first considered by Penrose. |71J 



5.3 General gravitational plane waves 

We have succeeded in explicitly constructing a family of exact defect configurations that 
converges to a specific exact gravitational wave solution of general relativity in the limit 
that the configuration becomes continuous. In the construction we made some explicit 
choices for the sake of simplicity. We chose to align the grid with the coordinate axes 
and as a result described a plane wavefront with "+" polarization. Rotating the grid 
simply results in a rotation of the d.T and dy terms in | |5.21[ |, producing a wavefront with 
a different polarization. For example, rotating the grid by 7r/4 produces a wavefront 
with "x" polarization. The result always is an impulsive plane wave with constant 
polarization. 

We could also take a different geometry for the grid. Equation | |5.8) shows that the 
components of the grid do not need to be perpendicular, nor do the components have to 
have matching amplitudes. As long as the curvature density of the positive curvature 
defects is equal to the curvature density of the negative curvature defects, the Einstein 
tensor will vanish in the continuum limit. Different geometries will not produce any 
wavefronts that could not be produced by a square grid. 

Gravitational waves with longer wave packets may be constructed as a succession 
of wavefronts. Each wavefront contributes to the Riemann curvature in the following 
way, 

R = aS{u — uq) [{du A (cos0 da; + sincf) dy))^ — (du A (sin0 da; — coacj) dy))^] , (5.23) 

where (j) is the direction of the polarization of the wavefront. 

A continuous succession of wavefronts with a and varying with u therefore pro- 
duces a wave with curvature, 

R = Q;(u)[(duA(cos(/)(u) da;+sin0(u) d?/))^ — (duA(sin0(w) dx — cos (f>{u) dj/))^] . (5.24) 

This is the most general curvature produced by exact gravitational plane waves of the 
type described by Bondi et al. in |16J . 
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5.4 Other gravitational waves 

We have shown that it is possible to construct general gravitational plane waves in our 
piecewise flat model of gravity. Since a general gravitational wave will locally resemble 
a plane wave, this suggests that it should be possible to construct any gravitational 
wave. This is further supported by the observations that any gravitational wave may 
be obtained in the linear limit of the model, and that its is possible to find a piecewise 
flat approximation to any solution of general relativity] 29] Although this is all very 
suggestive, it is not enough to conclude that any gravitational wave may be constructed 
in our piecewise flat model of gravity. It is therefore instructive to explicitly examine 
the construction of a non-planar wave in our model and see what new issues arise. 

As an example, we will consider cylindrical gravitational waves. In the case of cylin- 
drical symmetry the solutions of the vacuum Einstein equation can found by solving a 



single linear differential equation.] 39 lOOj Consequently, the complete set of possible 



solutions is known. In particular, it is known that it is impossible to form isolated 
impulsive wavefronts. Instead, impulsive wavefronts leave a wake of backscattered 
gravitational waves. Explicit solutions with impulsive wavefronts ]41 show that the 
form of this backscattered component can be quite complicated. This suggests that it 
may be too much to hope for a simple description of the backscattered component in 
terms of a piecewise flat geometryp^ Instead we will focus on constructing an impulsive 
cylindrical gravitational wave front, while we allow the backscattered component to be 
any defect configuration, i.e. not necessarily one with vanishing energy-momentum. 

Let us consider a single ingoing cylindrical impulsive wavefront located at it = o j^ 
which propagates into a Minkowski background (i.e. the Riemann tensor vanishes for 
u < 0). In the limit of large radii {v oo), the wavefront will resemble an impulsive 
plane wave. Our first order guess for the piecewise flat cylindrical wavefront therefore 
consists of a grid of perpendicular massless defects with positive and negative energy 
oriented in the axial and azimuthal directions and moving in the radial direction (see 



figure 5.4 1. 



Curving a flat grid to a cylinder introduces two new complications. The first is that 
due to simple geometric considerations the velocity of the azimuthal defects must be 
smaller than the velocity of the axial defects. The azimuthal defects, therefore, are 
not quite massless. However, as the grid is made finer, this effect disappears and the 
azimuthal defects again become massless in the continuum limit. The second effect 
is that the four defects of the grid meeting at each junction do not lie in one plane. 

^^This says nothing about whether a piecewise flat description exists. We simply observe that constructing 
one wiU be very complicated. 
^^We use cylindrical lightcone coordinates defined by u = and v = . 
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Figure 5.4: To first order, a piecewise flat cylindrical wave front should consist of axial 
(black) and azimuthal (white) defects moving inwards in the radial direction. The junction 
conditions imply that any additional radial defect (gray) should he present at each junction. 



A consequence of this is that the junction condition \5.5\ cannot be satisfied by these 
four defects. To satisfy \5.5\ a fifth defect pointing in the outward radial direction must 
be added to each junction (depicted as the gray lines in figure [5^ . Just like general 
relativity, the piecewise flat model does not allow isolated cylindrical wavefronts. 



Repeating the steps of section 5.2 we can find the piecewise flat metric associated to 
this configuration and calculate the continuum limit. The result is 



ds = 2dudv 



{u + V + au6{u)y 



dcf"- + (1 - aue{u)f<iz'^ 



(5.25) 



The corresponding energy-momentum tensor is 

1 . . o _ 1 



T = a- 



-5{u)du^ - 2a 



(1 — au){v + (1 + q)u) 



6(^u)Audv. 



(5.26) 



As expected the energy-momentum vanishes in the u < region, but does not in the 
TDackscattered' u > region due to the presence of the radial defects. However, the 
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energy-momentum also has a non-vanishing delta peak at m = for v ^ 1, i.e. it fails to 
model a cylindrical impulsive gravitational wavefront. This failure is the result of the 
two singular components of the Riemann tensor, 

R ^ AaSiu) (du A dz)^ - 4-S(u)(du A -^d^)^ + 9(u) (•••), (5.27) 

V V2 

scaling differently with v. The expected behaviour — as can be deduced from linear 
cylindrical waves, or by expanding the exact cylindrical solution in Q for small values 
of u — is for both curvature components to scale inversely with y/v. 

On the level of the piecewise flat geometry, this happens because the number of axial 
defects stays constant while the radius changes, causing their density to scale inversely 
with the radius. The density of the azimuthal defects, on the other hand, stays constant. 
To change the scaling of the impulsive front, the energy densities of the defects must 
change as the defects move inward. The only way this can happen is if the defects 
periodically emit new defects into the backscattered region. 

We implement this by having the ingoing impulsive front emit an outward moving 
lightlike front of perpendicular defects. A single cell of the piecewise flat geometry then 



takes shape shown in figure 5.5 By adjusting the frequency and the energies of the 



emitted perpendicular components we can achieve any scaling with v that we want. To 
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achieve the proper scaling for the curvature of a cylindrical wave front, we choose the 
emitted components such that the metric in the continuum limit becomes, 

^ 2dudv + + - + (1 _ (5.28) 

2 y/V 

The singular part of the corresponding energy-momentum tensor is 

T = + e{u){---). (5.29) 

The impulsive wavefront has thus vanished from the energy-momentum tensor. Yet, 
the singular part of the Riemann curvature, 

R = 4^(5(u) (du A dz)2 - 4-^S(u)(du A ^d^i)^ + e{u) (•••), (5.30) 

does not vanish, and has the expected scaling for a cylindrical impulsive wavefront. 
The described configuration indeed contains a cylindrical impulsive wavefront. 

The described piecewise flat geometry is just one way an impulsive cylindrical grav- 
itational wavefront could evolve in our piecewise flat model of gravity. Recall that the 
model fundamentally describes gravity coupled to matter As discussed in chapter |3j 
the dynamics of the model is incomplete]^ It needs to be complemented by additional 
rules which relate hoe the matter content interacts, much like general relativity needs to 
be complemented by equations of motion for the matter content. Different djmamical 
completions of the model lead to different evolutions of cylindrical wavefronts. The 
specific evolution described here corresponds a matter interaction where interacting 
gravitational waves can produce other types of matter. 

One might ask if there exists a dynamical completion of the model that corresponds 
to pure gravity with no matter present at all. Of course, this is of academic interest only 
since the physical world also contains matter. In such a completion the backscattered 
component of the evolving cylindrical wavefront must also satisfy the vacuum Einstein 
equation. This requires the addition of additional defects to the configuration. These 
will break the remaining residual symmetry of the piecewise flat geometry. As a conse- 
quence it will no longer be possible to describe the geometry as an infinite collection of 
identical piecewise flat cells, substantially adding to the complexity of the description. 
This breaking of symmetry indicates that cylindrical gravitational waves may not be the 
best setting to further study the djmamical completion of the piecewise flat model. 

Nonetheless, the study of the interaction of gravitational waves in the piecewise flat 
model may provide essential guiding insights for dynamically completing the model. 



■'^This is related to the different ways in wliich a piecewise flat geometry can be continued after the colUsion 
of two defects. 
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A more promising setting for this study is the interaction of two plane wavefronts. 
This problem is well documented in the case of general relativity coupled to various 
types of matter. ijlT 48 53 STj Moreover, the intrinsic plane symmetry is much more 
suited to a piecewise flat geometry. Furthermore, it is well known that the future of 
plane wave collisions can contain regions which are locally diffeomorphic to black hole 
solutions.] 53 81J Studying plane wave collisions may therefore offer crucial insight 
in one of the major open problems of the piecewise flat model; the construction of 
non-trivial stationary vacuum solutions such as black holes. However, the study of 
dynamical completions of the model is a major project on its own and will not be 
attempted here. 



5.5 Conclusions 

In this chapter we have extended our previous linear result from chapter|4]for obtarnrng 
gravitational waves in our piecewise flat model for gravity in 3+1 dimensions to an 
explicit construction of a family of exact piecewise flat configurations that approaches 
exact plane wave solutions of general relativity. This suggests that the construction 
of general gravitational waves may be possible. As an example, we have constructed 
the evolution of a cylindrical gravitational wavefront allowing the backscattering to 
produce other t5^es of matter 

Understanding the construction of gravitational waves is an important step towards 
showing that the proposed model can indeed serve as description of gravity in 3+1 
dimensions. Further study of the interactions of gravitational waves may provide es- 
sential guiding insights for completing the incomplete dynamics of the model. 
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Conclusions and Outlook 



In this thesis we have studied the piecewise flat gravity model in 3+1 dimensions first 
suggested by 't Hooft in |92|. The primary idea behind the model is to reproduce the 
local finiteness of gravity in 2+1 dimensions by requiring that the Riemann curvature 
tensor vanishes in empty regions of spacetime (i.e. areas where the energy-momentum 
is zero). 

This is achieved by studying a system of straight cosmic strings propagating in a 
3+1 dimensional spacetime. Geometrically these strings correspond to conical defects 
of the spacetime with codimension 2. Unlike other piecewise flat approaches to gravity 
like Regge calculus, these conical defects are interpreted as physical degrees of freedom. 
This implies that they must obey the rules of causality, which means that we can allow 
only defects that are non-spacelike. 

In chapter |2] we have described two distinct ways of describing an arbitrary config- 
uration of conical defects. Each has its own advantages and disadvantages. 

In the first approach we described a configuration of defects using the Poincare 
holonomies of loops around the defects. A general configuration of defects is then given 
as a group homomorphism of the fundamental group of the spacetime minus the defects 
to the Poincare group. The requirement that all defects are non-spacelike translates 
to the condition that there must exist a complete set of generators of the fimdamental 
group whose assigned holonomy is either rotationlike or a null rotation. 

One of the advantages of this approach is that the holonomies are gauge invariant 
quantities Moreover, the consistency conditions at junctions in the network of defects 
take a clear algebraic form in terms of the holonomies. A major downside is that the 



^^Strictly speaking, there is still a global ISO(3, 1) invariance connected to the choice of frame of the base 
point of the loops. 
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holonomies entangle the geometric data of the different defects in an essentially non- 
local way This makes this method of description very complicated for configurations 
consisting of more than a handful of defects. 

A second more important drawback is that the holonomies do not properly distin- 
guish between defects with positive and negative (Ricci scalar) curvature. This means 
that there may be geometrically distinct configurations of defects that have algebraically 
indistinguishable holonomies. Possibly even more disturbingly, we found in chapter |3] 
that there exist assignments of the holonomies that appear to be consistent algebraically, 
but which do not allow a consistent geometric realization. 

The second approach emphasized the piecewise flat structure of the spacetime rather 
than the defects. It described a spacetime filled with defects as a collection of (3+1)- 
dimensional polytopes which were glued together along their faces. This gives a much 
more local description of the geometry. Moreover, it fully distinguishes between differ- 
ent configurations. 

The downside of this geometric approach is that there are many piecewise flat struc- 
tures that describe the same configuration. This means that even though this description 
has only a finite number of degrees of freedom per unit volume, it still contains many 
unphysical (gauge) degrees of freedom. 

6.1 Collisions 

In chapter |3] we studied the collision of a pair of defects. Since the defects carry non- 
commuting holonomies they cannot simply pass through each other. Instead an interpo- 
lating configuration of defects must be found that continues the piecewise flat geometry 
after the collision. Various topological configurations of intermediate defects were ex- 
amined — ranging from a simple single defect continuation for orthogonal collisions to 
more complicated 'tetrahedral' continuations involving six new intermediate defects. 

A common feature of these continuations is that they introduce new degrees of 
freedom. The current physical requirements of the model do not uniquely determine 
values of all the new parameters. The model as it currently is, is therefore incomplete. 
Additional physical input is required to uniquely determine the evolution of a configu- 
ration through a collision. 

Although this is in some sense alarming, this indeterminism may be a blessing in 
disguise. The extra physical data needed to pick a continuation of a collision may allow 
for the inclusion of different types of interaction in the model. For example, if the defects 
carried additional conserved charges, this would put new constraints on the possible 
continuations. If no such freedom were present, and the evolution of the model were 
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completely determined, than the model would allow for only one type of matter content 
with no further room for model building. 

A more pressing matter, however, is the observation that certain collisions demand 
continuations which seem to be inconsistent with the principle of local causality — a 
fundamental tenet of the model. More specifically, we found that if the parameter^^ 
of a collision exceed the bounds set by equation | |3.64^ , then some of the junctions in 
the considered continuations move faster than the speed of light. This situation occurs 
when the mass densities of the defects are large, their relative velocity is near relativistic, 
and/or the defects are close to parallel. 

If a junction moves faster that the speed of light, then there exists a Lorentz frame 
in which the junction is instantaneous in time. In that frame the junction corresponds 
to a defect spontaneously splitting into multiple defects. Since this happens at arbitrary 
distances from the collision, the split will appear to have no local cause. Superluminal 
jimctions therefore seem to be in conflict with the local causality principle of our model. 



The generic argument of section 3.5.3 shows that the appearance of superluminal 
junctions is universal for all possible continuations. It therefore appears that our pro- 
posed piecewise flat model cannot consistently find continuations for all types of colli- 
sions. Can such collisions be avoided by a fortunate choice of initial conditions? If the 
bound ( |3.64) only involved the masses and the velocities of the defects this might seem 
a possibility. However, the bound set by | |3.64^ is also violated by low energy collisions 
of almost collinear defects. The author sees no plausible physical argument that would 
disallow the creation of nearly collinear pairs of defects. It therefore seems unlikely 
that the collisions leading to continuations with superluminal junctions can be avoided 
altogether. 

Is there then a way that the superluminal junctions that appear can be reconciled 
with local causality? As we argued above, this appears to be impossible. It is, how- 
ever, not completely beyond the realm of possibility. Even though the causal past 
of a superluminal junction does not contain the actual collision event, it will contain 
parts of both colliding defects before the collision. Consequently, the junction could 
have received information of the impending collision. This does require some sort of 
signalling mechanism. One (simplistic) example of such a mechanism would be for the 
surface where the wedge of spacetime was removed to form a conical geometry, to carry 
physical informationp^ One could then always choose the junctions to occur exactly at 
the intersection of the wedge with the colliding defect. The junctions would still be 
superluminal, but it would just signify the meeting of two physical entities (the wedge 



^^The mass densities of the colliding defects and their relative angle and velocity. 

^^The natural thing to do is to treat this as a mere constructional artefact without geometrical or physical 
significance. 
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and the defect). As such there would be no problem with causality, just like there is no 
problem with the tip of a laser beam flashing across the moon moving faster than the 
speed of light. Such a solution does however seem very artificial. 

Besides superluminal junctions, we also observed that for the particular continua- 
tions studied some of the intermediate defects became superluminal for extreme values 
of the collisions parameters. This obviously is inconsistent with the base requirement 
of our model that no superluminal defects appear. However, the general argument of 
section 3.5.3 does not imply that this must be the case for any continuation. In fact, there 



appears little geometric obstruction to constructing more complicated continuations to 
remove the superluminal defects. 

In all it seems that the appearance of superluminal junctions and defects signals that 
the proposed piecewise flat model of gravity is — in some fundamental way — flawed. 
Fixing this may require some fundamental modifications to the model. 

Another issue, which we have paid little attention to in this thesis, is related to the 
fact that each collision introduces a number of new defects. Because the inverse process 
(some configuration of defects cleanly splitting in two defects) is very unlikely, this will 
cause the spacetime structure to become finer over tim^^ ultimately approaching some 
fractal-like structure. This contrasts significantly with the situation in 2+1 dimensions. 
In 2+1 dimensions, defects almost never collide head-on, so that no such fragmentation 
takes place. 

In his original paper |92| 't Hooft suggested to counteract this by restricting the 
holonomies of the defects to some discrete subset of the Poincare group. This would 
significantly increase the probability of processes where defects disappear, and the hope 
was that at some point an equilibrium in the structure would be reached. As an added 
bonus, restricting to a discrete subset of the Poincare group would introduce a minimal 
angle at which defects can collide, getting rid of the problematic nearly collinear colli- 
sion — although it would do little to prevent high energy collisions. It is, however, not 
clear that this can be implemented in a consistent way, while finding continuations for 
all collisions. 



6.2 Matter content 

In the linearized limit of the model collisions do not pose problems, because at linear 
order defects can simply move through each other. In chapter |4] we utilized this limit to 
examine the continuum limit of our piecewise flat model of gravity. 

^^Unless compensated by a comparable spacetime expansion, but this would require extreme fine tuning. 
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In particular, we examined the type of energy-momentum distribution that could 
be generated by an arbitrary configuration of defects. The result was that the model 
could reproduce the same energy-momentum distribution as a gas of non-interacting 
point particles. That is, even though the fundamental excitations of our model are 
straight strings with an infinite extent, it is still possible to approximate the matter 
configurations produced by point masses. For this conclusion it was essential that the 
model include defects with both positive and negative curvature. 

Not surprisingly, we cannot reproduce interacting forms of matter in the linearized 
limit. This is exactly, because collisions are the only form of interaction between the 
defects, and these became trivial in the linear limit. Once we include higher order terms 
collisions come back into the game and will produce some form of interaction. 

At this point, however, it is unclear what types of interaction may be reproduced. It 
seems that the extra input needed to find a continuation for a collision may in fact give 
us some freedom in describing different types of interaction. 

6.3 Gravitational waves 

One thing that we observed in our study of the linearized model, was that it was possi- 
ble to find non-trivial configurations with vanishing energy-momentum. By calculating 
the linear perturbation to the metric caused by an arbitrary linearized configuration 
of defects, we were able to show that these configurations with vanishing energy- 
momentum correspond to linear gravitational waves. 

So, even though the model removed gravitational waves as fundamental local exci- 
tations, these waves appear as an emergent feature in the continuum limit of the model. 

In chapter |5] we have explicitly extended this conclusion to exact gravitational wave 
solutions by constructing a family of defect configurations that can approximate an 
arbitrary (exact) gravitational plane wave. Obtaining general exact wave solutions is 
much more difficult, but we sketched the outlines of the construction of a cylindri- 
cal gravitational wavefront from defects. This construction highlighted the fact that 
a cylindrical wave front in general relativity cannot appear without backscatterrng. The 
backscattered component of such a front appears to be very complicated, and will likely 
contain the collisions of lightlike defects. 

Collisions of lightlike defect occur exactly on the boundary of the collision parameter 
space where superluminal effects start to occur. This boundary case was not included 
in our general study of chapter |3| this needs further study. 

Moreover, the chaos of the backscattering of a cylindrical wave front is not the best 
environment to study the collisions of piecewise flat wavefronts. A more controlled 



107 



Chapter 6: Conclusions and Outlook 



setting is the collision of impulsive wave fronts. There is extensive literature on these 
in general relativity. [ ,28^ ,48^ ,53^ ,81-83 J One thing that is known about these is that the 
resulting spacetime often locally corresponds to a black hole metric. p2) Studying such 
collisions in piecewise flat gravity may therefore also give crucial insight in the con- 
struction of static gravitational fields using configurations of conical defects. 

6.4 Quantization 

In this thesis we have focused on the classical aspects of the proposed piecewise flat 
model of gravity in 3+1 dimensions. One of the motivations behind the model is that 
it might be easier to quantize due to the similarity with gravity in 2+1 dimensions. At 
the same this similarity suggests that the quantization of this model will be subject to a 
similar set of subtleties. In this section we will comment on how some of our conclusions 
affect the possible quantization of the model. 

First notice that indeed the model has only a finite number of degrees of freedom per 
unit volume. This means that its quantization is a form of ordinary quantum mechanics 
rather than quantum field theory, bypassing any problems with renormalization. 

We have explored different approaches to describing a configuration of defects. As 
with the 2+1 dimensional model of gravitating point particles, each of these will have 
their own implications for quantization. 

The first of these was the description of a configuration of defects by the assignments 
of Poincare holonomies to equivalence classes of loops aroimd the defects. In particular, 
we found that a single defect can be described by its Poincare holonomy. 

This description is strongly reminiscent of the description of a gravitating point 
particle in 2+1 dimensions using its Poincare holonomy. This suggests that like in 
2+1 dunensions,[(58| the Lorentz holonomy may act as the covariant momentum of the 
defect. Consequently, quantization using this type of description may lead to sim- 
ilar features as observed in |58|, namely discretization of time and space and non- 
commutativity of the position variables. 

This would obviously need to be checked more carefully using a d5mamical formu- 
lation of the model starting from an action. However, there also are some reasons to 
doubt the above conjecture. For one we know that the Lorentz holonomy of a defect 
encodes not only the energy and velocity of a defect, but also its orientation. The latter, 
normally would be considered a position rather than a momentum variable. 

Moreover, we know that the description of a configuration of defects in terms of 
its holonomies is not entirely unambiguous once we allow for defects with negative 
curvature. Furthermore, we also observed that there exist consistent assignments of 
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holonomies that do not correspond to a geometrically realizable configuration of de- 
fects. Both these observations suggest that the holonomies may not be the right vari- 
ables to use for quantization. 

The second approach to describing configurations of defects that we explored was 
in terms of piecewise flat manifolds. This description certainly is unambiguous, and 
due to its geometric nature there is no doubt about the geometric realization of any 
configuration. The downside of this approach is that it has a lot of residual gauge 
freedom; there are many different descriptions of the same configuration of defects. 

Any quantization using this description therefore needs to take into account these 
gauge degrees of freedom. This could be done by choosing a gauge fixing prior to 
the quantization. A partial approach to this could be to try to find a foliation of the 
piecewise flat cell complex by Cauchy surfaces in such a way that each slice is again 
a piecewise flat manifold and that the time in each cell runs at the same rate. This (if 
possible) should yield something that is very similar to 't Hooft's polygon model in 2+1 
dimensions — which one could refer to this as a polytope model. 

Regardless of the approach to describing a configuration of defects, any quantization 
of the model will require an action. In this thesis we have attacked most problems in 
the terms of kinematics, rather than formulating it in terms of dynamics as the result of 
some action. Nonetheless the common wisdom is that an action is required to properly 
formulate the dynamics of the model and its quantization. As an alternative 't Hooft 
has suggested | [92) that this model could be quantized using his "pre-quantization" 
formalism. [|9l| 

6.5 Applications outside of quantum gravity 

The problems with consistency of the model when considering collisions may signal 
that it is not suitable as a fundamental model of gravity in 3+1 dimensions. Nonetheless, 
the studies of chapters|4]and|5]have shown that the model can at least approximate some 
systems in general relativity. 

As such, the fact that this model describes these systems in terms of discrete excita- 
tions evolving with time, may prove useful in (numerical) analysis of such system. For 
example, the ability to decompose gravitational wavefronts in terms of lightlike defects 
opens up the possibility of studying non-uniform wavefronts. 

As another example, it is relatively easy to find defect configurations that approxi- 
mate FRW solutions of general relativity with equation of state parameter w — —1/3. 
Since these configurations consist of discrete defects, they are inhomogeneous at small 
scales. Such configurations may therefore serve as interesting test cases for approaches 
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to cosmological averaging. The use here may however be limited by the fact that these 
models may belong to a somewhat pathological subset of configurations due to the 
absence of gravitational fields between the defects. 
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Some facts about PSL(2, C) 



The group of proper orthochronous Lorentz transformations in 3+1 spacetime dimen- 
sions, SO (3, 1)+ is isomorphic to the group of projective special linear transformations 
of the complex plane, PSL(2, (D), also known as the group of Mobius transformations. 
Each Mobius transformation, 

(A.1) 

cz + a 

can be represented as a complex 2x2 matrix. 

This identification is unique up to an overall scalar factor. We can use this freedom to 
ensure that the determinant of the matrix, ad — be, is equal to one. This condition fixes 
the overall scalar factor up to a sign. We can therefore represent Lorentz transformation 
by special linear 2x2 matrices modulo an overall sign, which are known as projective 
linear transformations. In terms of the groups, 

SO(3,l)+ ^SL(2,C)/K2 = PSL(2,C). (A.3) 



This representation of the Lorentz group is often very useful when dealing with the 
algebraic conditions on holonomies of defects. This is mostly because of the practical 
reason that PSL(2, (D) matrices have only four (complex) components compared to the 
sixteen (real) components of a SO(3, 1) matrix. In this appendix we review some basic 
facts about this representation and clarify some of the conventions used. See chapter 
1 of f72] for an extensive discussion of this representation albeit with slightly different 
terminology 
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A.l Lorentz transformations 

Minkowski space can be identified with the space of complex Hermitian 2x2 matrices 
as follows 

{t,x,y,z)^X=\ A.4 
\x + ly t — z J 

Under this identification the determinant of X is equal to (minus) the Minkowski norm 

det X = t^ -x^ -y^ - z^. (A.5) 

Matrices Q in SL(2, (D) act on the space of Hermitian 2x2 matrices through the adjoint 
action 

AdgX ^QXQl (A.6) 

Since the determinant of Q is one, this action preserves the determinant of X. Conse- 
quently, using the identification I A.4| , this action can be interpreted as a linear norm 



preserving action of SL(2, (D) on Minkowski space. The action | |A.6[ | and the identifica- 
tion | |A.4[ | therefore define a group homomorphism from SL(2, (D) to the Lorentz group. 
This group homomorphism is in fact a 2-1 surjection onto the identity component of the 
Lorentz group with kernel {I, —I}. Consequently, we find that 

PSL(2, C) = SL(2, €)/{!, -1} ^ SO(3, 1)+. (A.7) 



A.2 Generators 



The Lie group PSL(2, (D) is generated by the three Pauli matrices 



won 1 /o -i\ 1 (1 

^^=2U oj' '"^^^i oj' ^^=2(,0 -1 

If we denote a = {ax, cry, cfz), then any element Q of PSL(2, (D) may be written a^ 

Q = exp(<T-(6 + if)). (A.9) 
This parametrization has some useful characteristics. If & = 0, then 

= (cxp (iff • r))^ = exp {-ia-f) ^Q-\ (A.IO) 

So, Q is unitary. Conversely, if Q is unitary, then b = 0. Consequently, the elements with 
6 = form a subgroup. This subgroup corresponds to the subgroup of rotations of the 



^'Note that this a non-trivial statement since the exponential map from s[(2, (D) to SL(2, <D) is not onto. 
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Lorentz group. An element with 6 = corresponds to a rotation about the f axis with 
angle |lr1|. 

Similarly, Q is hermitian if and only if r = 0. Hermitian elements of PSL(2,(D) 
correspond to Lorentz boosts in the 6 direction with rapidity ||&||. 

Consequently, a general rotationlike Lorentz transformation wiU correspond to a 
PSL(2,(D) matrix 

R = exp(CT • b) exp(ia • r) exp(— a • b) 
= exp (i exp(CT • b){a ■ r) exp(-a • b)) 

— ' X T* — » 

= exp ( sinh||6|| + i cosh||6|| f) , 
ll^ll 

with 6 • r = 0. Conversely, we can conclude that any 

Q = exp(a- (6 + ir)), (A.12) 

with 6 • r = and ||6|| < ||r|| corresponds to a rotationlike Lorentz transformation. 
Similarly, one finds that elements with 6 • r = and ||6|| > ||f|| correspond to boostlike 
Lorentz transformations. 

In the boundary case that 6 • f = and ||5|| — \\r\\ one finds that 

(a • {b + ir)f = [h + ir) ■ (6 + ir)Ic;2 

= (ll^f -|MI'+2i6-r)Ic2 (A.13) 
= 0. 

That is, (7 • (6 + if) is nilpotent. Consequently, the corresponding PSL(2, (D) element 
is imipotent, and by extension so is the corresponding Lorentz transformation. The 
eigenvalues of the corresponding Lorentz transformation thus are all equal to one; It is 
a null rotation. 



A.3 Conjugacy classes 

In the piecewise flat model of gravity studied in this thesis, it is often necessary to deter- 
mine the conjugacy class of a Lorentz transformation. In the PSL(2, (D) representation 
of the Lorentz group, the conjugacy class of an element can be determined on basis of 
the trace. 

Two matrices in SL(2, (D) belong to the same conjugacy class if and only if they have 
the same Jordan normal form. The Jordan normal form of a general 2x2 matrix, M, is 



Co 


'1 








Mr 


or 


(o 










(A.14) 
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Lets first consider the first case. If M is in SL(2,(D), then dctM = A1A2 = 1, and 
consequently A2 = 1/Ai. The trace is thus given by 

TrA/ = Ai + -'-, (A.15) 
Ai 

and we obtain 

A,, . (A.16) 

Furthermore, if M represents an element of PSL(2, (D) we can use the overall sign free- 
dom to fix Im Tr M > 0. 

In the second case, det 7\/ = A^ = 1, implies that A = ±1, and we can use the overall 
sign in PSL(2, C) to fix A = 1. As a result, the trace of M is 2. Plugging Tr A/ = 2 in 
equation ( AT6) we see that in the diagonalizable case this implies that AI is the identity. 

We therefore obtain the following lemma: 

Lemma 1. If Mi and M2 are non-identity elements of PSL{2, (D) represented as SL(2, (D) 
matrices with ImTr A/1,2 > or IniTrAf 1,2 = and ReTr Afi,2 > 0, then Mi and M2 
lie in the same conjugacy class if and only if It Mi — Tr Af2- 

If we write an element Q of PSL(2, (D) as 

Q = exp(<T-(6 + ir)), (A.17) 

then the trace is given by 



Trg = 2cosy ||r1| - ||6|| - 2ir ■ b. (A.18) 
Consequently, a PSL(2, (D) transformation Q not equal to the identity is ... 

• ... rotationlike (||f|| > ||6|| and f-b^O), if and only if 

ImTrQ^O and |TrQ|<2. (A.19) 

• ... boostlike {\\r\\ < \\b\\ and r-b — 0), if and only if 

IniTrQ^O and |TrQ|>2. (A.20) 

• ... a null rotation = ||6|| and f-b = 0), if and only if 

ImTrQ^O and |TrQ|=2. (A.21) 

This characterization of elements of PSL(2, (D) is very useful when we need to test if a 
defect described by a certain holonomy is tunelike, spacelike, or lightlike. 
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De modeme natuvirkiinde rust op twee standaardtheorieen. Aan de ene kant geeft de 
kwantumveldentheorie aan hoe materia zich gedraagt op (sub)atomaire schaal. Aan 
de andere kant is er de algemene relativiteitstheorie van Albert Einstein, die de gra- 
vitationele interactie tussen planeten, sterren en hele melkwegstelsels beschrijft. Beide 
theorieen zijn in de afgelopen eeuw tot de kleinste details van hun voorspellingen na- 
gemeten en hebben iedere experimentele toets waaraan ze zijn onderworpen glansrijk 
doorstaan. 

Echter, als beide theorieen gelijktijdig op een systeem worden toegepast, dan lijken 
ze elkaar tegen te spreken. In experimenten leidt dit tot nu toe niet tot problemen omdat 
de inconsistenties pas zichtbaar worden in systemen die tegelijkertijd heel erg klein en 
heel erg zwaar zijn. De dichtheid waarbij beide theorieen belangrijk worden, heet de 
Planckdichtheid, ongeveer 10^^ kilogram per kubieke meter. Dit is de dichtheid die 
wordt bereikt als je de massa van het hele zichtbare heelal samendrukt tot de grootte 
van een atoomkem. 

Systemen met een dergelijke dichtheid kunnen natuurlijk niet in een lab gereali- 
seerd worden, en het is daarom niet zo gek dat de tegenstrijdigheden tussen deze twee 
theorieen niet naar voren komen in experimenten. Toch zijn er in de natuurkujn.de 
waarschijnlijk situaties waarin zulke extreme omstandigheden voorkomen. De alge- 
mene relativiteitstheorie voorspelt bijvoorbeeld dat de materie in het centrum van een 
zwart gat tot een punt wordt samengedrukt. Het omgekeerde doet zich voor in het 
vroege heelal,dat volgens dezelfde theorie vanuit een toestand met oneindige dichtheid 
is geexpandeerd tot het heelal van nu. 

Om de wisselwerking tussen deeltjes in zulke extreme omstandigheden te beschrij- 
ven is een theorie nodig die zowel de kwantumveldentheorie als de algemene relati- 
viteitstheorie omvat. Een dergelijke overkoepelende theorie wordt ook wel een theorie 
van kwantumgravitatie genoemd. De moderne theoretische natuurkunde is al meer dan 
een halve eeuw op zoek naar een dergelijke theorie. Tot nog toe zonder een definitief 
succes. 
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Een centraal probleem bij het opstellen van een overkoepelende theorie is dat als 
de conventionele technieken om een kwantumtheorie te formuleren worden toegepast 
op de algemene relativiteitstheorie, deze leiden tot onacceptabele oneindigheden in de 
voorspellingen. Deze oneindigheden vinden hun oorsprong in de geometrische aard 
van de algemene relativiteitstheorie. In dit proefschrift wordt een altematief model 
voor zwaartekracht ontwikkeld dat dit probleem probeert te vermijden. 

Zwaartekracht is kromming van de ruimtetijd 

De algemene relativiteitstheorie kan in een notendop worden samengevat als: "Zwaar- 
tekracht is kromming van de ruimtetijd." Wat wordt daarmee bedoeld? 

We beschouwen eerst het begrip kromming. Stel een wereld met maar twee ruim- 
telijke dimensies, zoals Flatland (Platland) beschreven door Edwin A. Abbott. En 
stel, dat die platlanders leven op het oppervlak van een bol (in plaats van op een plat 
vlak zoals in het werk van Abbott). Hoe zouden de twee-dimensionale inwoners dat 
kunnen zien en meten? Een mogelijkheid is het meten van de hoeken van geometrische 
figuren, een favoriete bezigheid van de inwoners van Platland. In tegenstelling tot de 
gebruikelijk constatering dat de som van de hoeken van een driehoek 180 graden is, 
zullen de platlanders constateren dat deze som afhangt van de oppervlakte van de 
driehoek; des te groter het oppervlak van de driehoek, des te groter de som van de 
hoeken. Dergelijke afwijkingen van de gebruikelijke Euclidische meetkunde worden 
door wiskundigen samengevat onder de algemene noemer 'kromming'. 

Een ander effect dat de platlanders zouden kunnen waarnemen, is dat op een bol 
parallelle lijnen elkaar in twee punten snijden. Nu liggen deze punten op tegenover- 
gestelde punten van de bol. Dus als de bol erg groot is, zullen ze dit niet snel merken. 
Maar, een oplettende platlander zal zien dat twee passanten, die hem in een rechte lijn 
parallel aan elkaar voorbij lopen, na een tijdje dichter naar elkaar toe zijn bewogen. 
Mocht de platlander in de veronderstelling verkeren in een plat vlak te leven, dan zou 
hij concluderen dat er kennelijk een soort kracht op de twee passanten werkt die ze naar 
elkaar toe trekt; een soort zwaartekracht. 

Volgens de algemene relativiteitstheorie is de zwaartekracht in onze drie-dimens- 
ionale wereld een vergelijkbaar effect. Alleen gaat het niet om kromming van alleen 
de ruimte, maar om kromming van de ruimte en de tijd samengenomen tot een vier- 
dimensionale ruimte, de ruimtetijd. De banen van planeten om de zon zijn volgens de 
relativiteitstheorie in werkelijkheid rechte lijnen in deze ruimtetijd. Maar, doordat wij 
bij onze waamemingen (impliciet) veronderstellen dat de ruimtetijd 'vlak' is, lijken deze 
banen voor ons gekromd en concluderen wij dat er kennelijk een kracht op de planeten 
werkt: de zwaartekracht. 
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De algemene relativiteitstheorie stelt verder dat er een wisselwerking bestaat tussen 
de kromming van de ruimtetijd en massa. Aan de ene kant bepaalt de kromming hoe 
massa's door de ruimtetijd bewegen, aan de andere kant bepaalt de aanwezigheid van 
massa hoe de ruimtetijd zich kromt. Dit laatste wordt vastgelegd door de zogenaamde 
Einsteinvergelijking. 

De kromming wordt niet volledig vastgelegd door de Einsteinvergelijking. Een deel 
van de kromming kan zich vrijelijk bewegen, als een soort golf door de ruimte. Deze 
golven worden zwaartekrachtsgolven genoemd. Het zijn juist deze zwaartekrachtsgol- 
ven die zorgen voor (niet renormaliseerbare) oneindigheden bij het opstellen van een 
kwantummechanische tegenhanger van de algemene relativiteitstheorie. 

Zwaartekracht in twee dimensies 

Als we een theorie niet goed begrijpen, kan het leerzaam zijn om te bekijken hoe de 
theorie zich gedraagt als we het aantal ruimtelijke dimensies veranderen. Dit geldt ook 
voor de algemene relativiteitstheorie. 

In een wereld met maar twee ruimtelijke dimensies zoals Platland — en dus met een 
drie(=2+l)-dimensionale ruimtetijd — is de algemene relativiteitstheorie veel eenvoudi- 
ger In tegenstelling tot de theorie in 3+1 dimensies, legt de Einsteinvergelijking in 2+1 
dimensies de kromming van de ruimtetijd geheel vast in termen van de massaverdeling: 
Platland kent geen zwaartekrachtsgolven. 

De kromming bevindt zich alleen op die pun ten in de ruimtetijd waar massa is. Lege 
stukken ruimtetijd zijn geheel vlak. De kromming op de plaats van een puntdeeltje 
is een zogenaamd 'conisch defect'. Dat wil zeggen: op de plek van de puntdeeltje 
ontbreekt er een taartpunt uit de ruimtetijd zoals bij de punt van een kegel (zie figuur 
1.1 op bladzijde|4}. De grootte van de ontbrekende hoek is evenredig met de massa van 
het pimtdeeltje. 

De geometric van een ruimtetijd van een verzameling puntdeeltjes kan gezien wor- 
den als allemaal blokjes vlakke ruimtetijd die langs hun vlakken aan elkaar geplakt 
zijn. De puntdeeltjes bewegen zich dan langs de ribben van deze blokjes. Een dergelijke 
constructie wordt een stuksgewijs vlakke geometric genoemd. 

Omdat ruimtetijd waar zich geen massa bevrndt helemaal vlak is, is er in 2+1 dimen- 
sies geen vrijheid voor zwaartekrachtsgolven. Dit doet vermoeden dat het formuleren 
van een kwantummechanische versie van de algemene relativiteitstheorie makkelijker 
is in 2+1 dimensies. Dit werd voor het eerst bevestigd door Edward Witten in 1988.| |102) 
Sindsdien wordt gravitatie in 2+1 dimensies veel gebruikt als voorbeeldtheorie om 
conceptuele problemen in kwantumgravitatie te onderzoeken. 
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Stuksgewijs vlakke zwaartekracht in drie dimensies 

Het idee achter het model dat is beschreven in dit proefschrift is het construeren van 
een model voor zwaartekracht in drie ruimtelijke dimensies (en dus 3+1 ruimtetijddi- 
mensies) dat de stuksgewijs vlakke structuur van algemene relativiteitstheorie in 2+1 
dimensies bezit. Omdat zo'n model geen zwaartekrachtsgolven zal bezitten, is het 
mogelijk makkelijker om een kwantummechanische tegenhanger te formuleren. In dit 
proefschrift gaan we echter niet verder op de mogelijkheid van het formuleren van 
een kwantummechanische versie van het model, maar richten we ons op de klassieke 
eigenschappen en interne consistentie van het model. 

Het model wordt bepaald door drie basis regels: 

1. De ruimtetijd is vlak op plekken waar geen materie is. 

2. De gehele ruimtetijd (zowel de regio's met als zonder materie) voldoet aan de 
Einstein vergelijking. 

3. Alle materie voldoet aan lokale causaliteit. (Elke gebeurtenis kan voorspeld wor- 
den door de gebeurtenissen in zijn verleden). Materie kan niet sneller bewegen 
dan het licht. 

De combrnatie van regels 1 en 2 legt een strikte beperkrng op aan de soort materie die we 
toe kunnen staan. De enige vorm van materie die is toegestaan, zijn rechte staven met 
een constante dichtheid en snelheid. De derde regel stelt vervolgens dat deze staven 
maximaal met de lichtsnelheid mogen reizen. 

De geometric van zo'n rechte staaf is wederom een conisch defect. De gehele geome- 
tric van een collectie staven bestaat uit vlakke (3+l)-dimensionale brokken ruimtetijd 
die langs hun ((2+l)-dimensionale) rand en aan elkaar geplakt zijn. De rechte staven 
bevinden zich vervolgens op de (l+l)-dimensionale grensvlakken. Dit is wederom een 
stuksgewijs vlakke geometric. 

Botsingen 

Een nieuw fenomeen ten opzichte van gravitatie in 2+1 dimensies is, dat in dit model 
de staven kunnen botsen. Hoofdstuk|3]richt zich op de vraag wat er gebeurt als twee 
staven met elkaar botsen. Het is geometrisch niet toegestaan dat de staven door elkaar 
heen bewegen. In plaats daarvan moet uit de botsrng een nieuwe stuksgewijs vlakke 
geometric ontstaan met nieuwe staven. Dit blijkt altijd mogelijk te zijn, echter met een 
aantal op- en aanmerkingen. 
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Ten eerste, is de stuksgewijs vlakke geometrie die ontstaat na een botsing niet uniek: 
er bestaat een enorm scala aan mogelijkheden. De regels die we aan het model hebben 
opgelegd, leggen de d5mamica van het model niet volledig vast. Dit is vergelijkbaar 
met de situatie in de algemene relativiteitstheorie waar de Einsteinvergelijku-ig aange- 
vuld moet worden met een bewegingsvergelijking voor de aanwezige materie om de 
dynamica volledig te bepalen. Om een volledig model te krijgen zal ons stuksgewijs 
vlakke model voor zwaartekracht aangevuld moet worden met regels die zeggen welke 
mogelijke vervolgtoestand gekozen moet worden bij elke botsing. 

Ten tweede blijken er extreme botsingen te bestaan, waarvoor alle mogelijk vervolg- 
toestanden knooppunten van staven bevatten die sneller bewegen dan het licht. Dit 
lijkt in strijd met onze lokale causaliteitsregel en stelt daarmee de consistentie van dit 
model ter discussie. Deze botsingen doen zich voor als de botsende staven extreem hoge 
energie hebben en/ of bijna parallel zijn. De extreemheid van deze omstandigheden 
doet vermoeden dat er mogelijk een verband is met het ontstaan van singulariteiten 
(zoals zwarte gaten) in de algemene relativiteitstheorie. 

Zwaartekrachtsgolven 

Een directe consequentie van regel 1 van het stuksgewijs vlakke model voor zwaarte- 
kracht is dat er geen zwaartekrachtsgolven zijn toegestaan. Dit was onze opzet. Het 
is echter ook een probleem als we willen dat het model de werkelijkheid beschrijft. 
Er bestaat namelijk sterk (indirect) bewijs dat er zwaartekrachtsgolven bestaan in de 
natuur. 

Het model biedt echter een uitweg. In prrncipe staat het model zowel staven met po- 
sitieve energie als staven met negatieve energi^^toe. Dit geeft de mogelijkheid tot het 
creeren van configuraties waarin (gemiddeld genomen) de energie van de staven met 
positieve energie opgeheven wordt door de energie van staven met negatieve energie. 
De fysische interpretatie van een dergelijke configuratie is een puur zwaartekrachts- 
veld. Deze mogelijkheid is onderzocht in hoofdstukken|4]en|5] 

In hoofdstuk|4]is de limietsituatie van het model onderzocht waarin alle staven heel 
erg licht zijn. In deze limiet vereenvoudigt de beschrijvuig van het model sterk doordat 
alle interactie tussen de staven verdwijnt. Het is daardoor mogelijk een beschrijving te 
geven van een algemene configuratie van staven. 

Met deze beschrijving kunnen we uitrekenen wat de energie-impulsverdeling van 
een algemene configuratie is. Deze komt overeen met de energie-impulsverdeling van 

^^De staven representeren een combinatie van materie- en gravitatievrijheidsgraden en kunnen daardoor 
effectief negatieve energie hebben. 
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een gas van puntdeeltjes zonder interactie. Dit is opmerkelijk aangezien onze funda- 
mentele bouwstenen, de staven, oneindig lang waren. Tevens blijkt uit de berekening 
dat er inderdaad een scala aan configuraties bestaat, die netto geen energie en impuls 
bezitten. Door uit te rekenen wat het effect (tot op eerste orde) op de geometrie van een 
algemene configuratie is, hebben we kunnen laten zien dat deze configuraties zonder 
energie en impuls precies overeenkomen met de vacuiimoplossingen van de gelineari- 
seerde Einsteinvergelijking. Dat wil zeggen, ze beschrijven (lineaire) zwaartekrachts- 
golven. 

In hoofdstuk|5]is vervolgens aangetoond dat het in het volledige model mogelijk is 
configuraties te construeren die een vlakke zwaartekrachtsgolf willekeurig dicht bena- 
deren. Tevens wordt een opzet gegeven voor de constructie van cilindrische zwaarte- 
krachtsgolven. Dit is een sterke aanwijzrng dat de constructie van willekeurige zwaar- 
tekrachtsgolven mogelijk is binnen het model. 

Dus hoewel het model op de fundamenteel schaal van individuele vlakke stukjes 
ruimte geen plaats biedt voor zwaartekrachtsgolven, ontstaan deze effectief toch als 
men kijkt naar grote schalen waar configuraties bestaan uit heel veel kleine vlakke 
stukjes. 

Conclusies 

De formulering van een stuksgewijs vlak model voor zwaartekracht in 3+1 dimensies 
is een gedeeltelijk succes. Het model dat we hebben beschreven heeft de eigenschap- 
pen die we wilden op korte afstand, maar lijkt op grote schalen toch zwaartekracht te 
beschrijven zoals we die kennen uit de algemene relativiteitstheorie. 

Het model is echter niet compleet en heeft extra regels nodig om de dynamica van 
een systeem volledig te beschrijven. Verder is het model mogelijk niet consistent bij 
hoge energieen. Dit laatste is mogelijk een indicatie dat dit model niet zonder verdere 
modificaties geschikt is als basis voor een fundamentele theorie voor kwantumgravita- 
tie. 
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